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Coupled-cluster  based  methods  have  become  widely  used  for  calculating  energies  and 
properties  of  molecular  ground  states.  In  this  dissertation,  I  present  a  new  method,  the 
partitioned  equation-of-motion  coupled-cluster  method  based  on  a  second  order  many- 
body  perturbation  theory  ground  state  (P-E0M-MBPT(2))  for  calculating  energies  of 
molecular  electronic  excited  states.  I  also  describe  how  to  calculate  gradients  for 
two  excited  state  methods,  P-E0M-MBPT(2)  and  the  similarity  transformed  equation- 
of-motion  coupled-cluster  singles  and  doubles  (STEOM-CCSD)  method.  These  new 
methods  are  then  applied  to  several  example  molecules. 
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CHAPTER  1 
AN  OVERVIEW  OF  CURRENT  METHODS  FOR  EXCITED  STATES 

Coupled-Cluster  Methods 

Coupled-cluster  theory'"''*  provides  a  routinely  applicable  method  for  calculating 
energies  and  properties  of  molecules  in  the  lowest  state  of  a  given  spin  and  symmetry. 
By  adding  triple,'^-  '^  quadruple,'^'  '^  and  higher  excitations,  it  is  possible  to  converge  on 
the  exact  answer  within  a  given  basis  set.  The  perturbative  inclusion  of  triple  excitation 
effects'^"^'  is  now  routine,  meaning  that  energies  and  properties  can  approach  near 
chemical  accuracy  for  many  molecules.  Unfortunately,  these  methods  do  not,  in  general, 
work  for  electronic  excited  states. 

The  primary  problem  with  describing  excited  states  with  normal  coupled-cluster  (CC) 
theory  is  its  inherent  single  reference  nature.  The  wavefunction  for  the  state  can  be 
expressed  as  a  linear  combination  of  determinants,  with  the  reference  determinant  having 
a  coefficient  of  one.  Although  in  principle  coupled-cluster  theory  can  be  made  exact,  the 
cost  of  such  a  calculation  would  increase  in  cost  at  the  same  rate  as  full  configuration 
interaction  (CI),  i.e.  as  the  size  of  the  system  factorial.^^  Therefore,  the  set  of  excitations 
used  in  the  calculation  are  truncated.  In  such  a  truncated  CC,  a  second  determinant  with  a 
coefficient  of  more  than  0.2  to  0.3  can  cause  the  accuracy  of  the  computed  results  to  drop 
significantly.  However,  for  an  open-shell  singlet  excited  state,  the  second  determinant 
must  have  a  coefficient  of  one. 

Coupled-cluster  theory  has  been  extended  in  order  to  be  able  to  describe  such 
multireference  states.    Multireference  coupled-cluster  theories  can  be  divided  into  two 
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main  categories.  They  are  the  Hilbert  space  approaches^^"^'  and  the  Foclc  spaces 
approaches. ^^"^^  A  very  simplified  view  of  Hilbert  space  multireference  coupled-cluster 
theory  is  that  each  determinant  in  the  active  space  gets  its  own  set  of  T  amplitudes, 
and  a  small  effective  Hamiltonian  is  diagonalized  over  the  space  of  these  wavefunctions. 
In  the  Fock  space  coupled-cluster  method,  a  single  exponential  operator  is  used  for  all 
determinants  in  the  active  space,  but  this  exponential  contains  operators  that  change  the 
number  of  electrons  in  the  system.  Hilbert  space  coupled-cluster  theory  is  more  suited 
to  studying  one  or  a  few  states,  while  Fock  space  coupled-cluster  theory  is  more  suited 
to  studying  energy  differences  between  many  states  and  to  calculating  spectra. 

An  alternative  method  for  calculating  excitation  energies  within  a  coupled-cluster 
framework  is  to  follow  the  response  of  the  system  to  an  external  time  dependent 
perturbation.  This  was  introduced  by  Monkhorst,"^^  and  has  been  developed  into  the 
coupled-cluster  linear  response  (CCLR)  method.'^^"^'^  The  same  equations  can  also  be 
derived  from  an  equations-of-motion^^  framework  based  on  a  coupled-cluster  ground 
state.  This  gives  the  equation-of-motion  coupled-cluster  (EOM-CC)  method.^^^^  In 
EOM-CC  (or  CCLR)  the  ground  state  coupled-cluster  operator  is  used  to  perform  a 
similarity  transformation  of  the  Hamiltonian,  which  is  then  diagonalized  over  a  space 
of  excited  determinants.  EOM-CC  and  CCLR  are  identical  for  excited  state  energies, 
but  the  approximations  made  by  the  two  methods  differ  for  oscillator  strengths  and 
polarizabilities.  Typically,  the  ground  state  operator  and  the  space  of  excited  determinants 
is  truncated  to  include  only  single  and  double  excitations.  This  gives  the  EOM-CCSD^^'  ^^ 
(or  CCSDLR)^^'  ^"^  method.  Models  have  also  been  developed  which  partially  include 
the  effects  of  triple  excitations. ^^^^ 
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Recently,  Nooijen  and  Bartlett  have  developed  a  new  method  for  calculating  excited 
state  energies  and  properties.  It  is  the  similarity  transformed  equation-of-motion  coupled- 
cluster  (STEOM-CC)  method.^^6  For  singly  excited  states,  STEOM-CC  is  closely 
related  to  the  Fock  space  coupled-cluster  method,^^  but  conceptually  they  are  very 
different.  In  STEOM-CC,  ionization  potential  EOM-CCSD^''  and  electron  attachment 
EOM-CCSD^^  calculations  are  done.  Information  about  the  states  with  one  fewer  and 
one  more  electron  is  then  used  to  perform  a  second  similarity  transformation  on  the 
Hamiltonian^^  such  that  the  terms  which  directly  couple  single  and  double  excitations 
are  set  to  zero.  A  detailed  description  of  STEOM-CC  will  be  given  in  Chapter  4. 

y  •  . 

Other  Current  Excited  State  Methods 

At. 

An  attempt  to  catalog  all  ab  initio  methods  currently  used  for  calculating  excited 
states  would  be  futile  and  will  not  be  attempted.  Instead,  several  of  the  methods  which 
are  used  most  often  will  be  reviewed.  But  first  it  will  be  useful  to  take  a  more  general 
view  of  excited  states  and  break  down  how  the  various  methods  describe  excited  states. 
In  this  view,  the  excited  state  energy  is  divided  into  several  components.  They  are  the 
zeroth  order  description  of  the  wavefunction,  the  orbital  relaxation,  and  the  dynamic 
correlation  in  the  excited  state.  • 

Properly  describing  the  zeroth  order  description  simply  means  being  able  to  include 
all  of  the  determinants  which  play  a  dominant  role  in  the  excited  state.  For  example, 
this  means  being  able  to  give  the  two  determinants  in  an  open-shell  singlet  the  same 
weight.  The  orbital  relaxation  refers  to  how  much  the  orbitals  are  allowed  to  rearrange 
themselves  in  response  to  the  excitation. 
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The  dynamic  correlation  for  the  excited  state  is  further  broken  down  into  the  dynamic 
correlation  for  the  ground  state  and  the  differential  dynamic  correlation  between  the 
ground  and  excited  states.  The  reason  for  this  division  of  the  dynamic  correlation  is  that, 
especially  for  the  single  reference  theories,  first  a  ground  state  is  calculated  and  then  the 
excited  states  are  built  upon  it.  The  justification  for  this  division  is  that  most  electrons 
do  not  take  part  in  the  excitation  and  therefore  still  interact  with  each  other  in  the  same 
way  in  the  excited  state  as  in  the  ground  state. 

These  definitions  are  not  unique,  especially  when  considering  orbital  relaxation  versus 
differential  dynamic  correlation.  The  coupled-cluster  viewpoint  will  be  taken  here.  The 
key  to  the  coupled-cluster  viewpoint  is  Thouless'  theorem,^^  which  states  that  any  two 
non-orthogonal  determinants  can  be  written  as 

|$,;)  =  e^'|$j),  (1-1) 

where  T\  is  a  single  excitation  operator.  Therefore,  any  set  of  single  excitations  or  any 
product  of  single  excitations  will  be  considered  orbital  relaxation,  and  dynamic  correlation 
will  be  reserved  for  "connected"  double  and  higher  excitations.'^ 

Table  1-1  provides  a  convenient  representation  for  how  the  selected  methods  break 
down  into  the  four  terms.  A  'V  means  that  the  method  treats  this  term  adequately. 
Two  pluses  mean  that  the  method  treats  the  term  very  well.  A  "-"  means  that  the  term 
is  included,  but  poorly.  Finally,  a  "0"  means  that  the  term  is  completely  neglected. 
The  question  mark  means  that  how  well  the  term  in  question  is  treated  varies  based  on 
specifics  of  the  calculation. 

The  first  method  to  be  considered  is  probably  the  simplest  of  all  ab  initio  excited  state 
methods.  It  is  the  Tamm-Dancoff  approximation,^'- "^^  also  known  as  the  configuration 


^'i'.W*'' 


Table  1-1:  An  overview  of  excited  state  methods. 


zeroth  order 
description 


orbital 
relaxation 


ground  state 

dynamic 

correlation 


differential 

dynamic 

correlation 


CIS 

+ 

- 

0 

0 

RPA 

+ 

+ 

0 

0 

MCSCF 

++(?) 

++ 

- 

- 

CASPT2 

++(?) 

++ 

+ 

+ 

CIS(D) 

+ 

- 

+ 

+ 

OM-CCSD 

++ 

+ 

++ 

+ 

STEOM- 
CCSD 

+ 

+ 

■++    ^ 

++  (?) 

MRCI 

++(?) 

++(?) 

++(?) 

++(?) 

interaction  singles  (CIS)^-^"^^  method.  The  names  will  be  used  interchangeably.  In 
CIS  the  excited  state  is  determined  from  diagonalizing  a  Hamiltonian  matrix  made  up 
of  single  excitations  with  respect  to  the  Hartree-Fock  reference  determinant.  Since,  by 
Brillouin's  theorem,^^  single  excitations  cannot  couple  with  the  Hartree-Fock  determinant, 
this  diagonalization  leaves  the  ground  state  unchanged  and  gives  excited  states.  Because 
all  possible  single  excitations  are  included,  any  excited  state  which  corresponds  to  a  single 
electron  being  promoted  can  be  described.  However,  the  orbital  relaxation  is  extremely 
limited  and  the  method  contains  no  dynamic  correlation. 

The  next  method  is  the  random  phase  approximation  (RPA).^^'  ^^  Although  RPA  can 
be  derived  several  ways,  the  one  most  useful  here  is  to  consider  RPA  as  the  response 
of  a  wavefunction  written  as'^- 13  |^^  ^  (i  +  7^^  _,_  ^tI)\^q).  Clearly,  then,  RPA  also 
contains  no  dynamic  correlation,  but  it  does  contain  more  orbital  relaxation  than  CIS. 

The  next  method  to  be  considered  is  multiconfigurational  SCF  (MCSCF).^^-  ^^  In 
MCSCF,  a  set  of  configurations  is  chosen.  Then  the  coefficient  for  each  configuration 
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is  variationally  optimized  while  simultaneously  optimizing  the  orbitals  from  which  the 

configurations  are  built.  Since  the  procedure  involves  variationally  choosing  the  optimal 
orbitals,  orbital  relaxation  can  be  treated  quite  well.  Because  of  the  flexibility  to  choose 
any  configurations  wanted,  any  excited  state,  no  matter  its  excitation  level,  can  be 
described.  Thus  the  zeroth  order  description  can  also  be  very  good  if  the  configurations 
are  chosen  properly.  Because  of  the  multideterminantal  nature  of  the  final  wavefunction, 
some  dynamic  correlation  is  recovered,  but  normally  the  amount  is  small.  A  special  way 
of  choosing  the  configurations  is  to  choose  a  set  of  orbitals  and  to  create  all  possible 
configurations  with  a  given  number  of  electrons  distributed  among  that  set  of  orbitals. 
This  gives  rise  to  the  complete  active  space  SCF  (CASSCF)^^  variant  of  MCSCF. 

The  CASSCF  wavefunction  can  be  used  as  the  zeroth  order  wavefunction  in  a 
multireference  perturbation  theory  expansion.  The  most  popular  of  the  methods  derived 
this  way  is  the  CASPT2  method.^''  ^^  Since  it  is  based  on  a  CASSCF  wavefunction,  the 
arguments  above  about  the  zeroth  order  description  and  the  orbital  relaxation  still  apply. 
The  perturbation  theory  then  allows  a  treatment  of  the  dynamic  correlation.  Overall, 
this  method  works  well,  but  because  of  issues  involving  properly  choosing  the  active 
space  and  a  potential  intruder  state  problem,  it  is  sometimes  difficult  to  get  consistently 
good  results. 

Starting  with  the  CIS  wavefunction  being  used  as  the  zeroth  order  wavefunction  in  a 
perturbation  expansion,  it  is  possible  to  derive  the  CIS(D)  method.^^  CIS(D)  shares  CIS's 
weakness  in  term  of  orbital  relaxation,  but  it  does  allow  for  a  second  order  perturbation 
description  of  both  the  ground  state  dynamic  correlation  and  the  differential  dynamic 
correlation. 
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The  next  method  to  be  considered  is  EOM-CCSD,^^'  ^^  discussed  above.  Since  the 
excited  state  wavefunction  consists  of  both  single  and  double  excitations,  it  can  have  a 
proper  zeroth  order  description  of  not  only  singly  excited  states,  but  also  doubly  excited 
states  and  states  which  are  mixtures  of  singles  and  doubles.  On  the  other  hand,  since  the 
reference  state  determines  the  orbitals,  they  are  fixed,  and  the  excited  state  wavefunction 
has  only  limited  flexibility  to  relax  them  through  the  other  singles  and  through  the  doubles 
acting  as  a  single  times  a  single.  This  linear  excited  state  also  limits  to  what  extent 
differential  dynamic  correlation  is  included.  More  is  probably  included  than  for  any  other 
method  considered  so  far,  but  in  EOM-CCSD,  unlike  CASPT2,  the  differential  dynamic 
correlation  needs  to  compensate  for  the  orbital  relaxation.  The  ground  state  dynamic 
correlation  is  very  well  described  through  the  ground  state  coupled-cluster  calculation, 
and  therefore  for  states  without  large  orbital  relaxation  and  differential  correlation,  EOM- 
CCSD  works  well. 

In  STEOM-CCSD,^'  ^^  the  use  of  a  singles  only  excited  state  operator  limits  the 
zeroth  order  description  to  those  kinds  of  states  describable  in  CIS.  The  orbital  relaxation 
is  also  somewhat  limited,  but  some  orbital  relaxation  comes  indirectly  from  the  con- 
sideration of  the  ionized  and  electron-attached  states.  The  ground  state  coupled-cluster 
calculation  again  describes  the  ground  state  dynamic  correlation  very  well.  The  differen- 
tial dynamic  correlation  is  handled  through  the  second  similarity  transformation,  which  is 
built  from  an  active  space  of  ionized  and  electron-attached  states.  How  well  this  works 
in  practice  for  a  large  variety  of  molecules  is  not  clear  yet,  but  the  amount  of  differential 
correlation  included  will  depend  heavily  on  the  completeness  of  the  active  space.^^ 
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Finally,  there  is  multireference  CI  (MRCI).^"^  Because  of  the  flexibility  offered  by 
a  MRCI  calculation,  it  is  possible  to  achieve  extremely  good  accuracy  by  properly 
accounting  for  all  four  parts  of  the  excited  state  energy.  But  MRCI  has  two  very 
significant  drawbacks.  It  often  takes  an  experienced  practitioner  to  properly  build  the 
trial  wavefunction,  and  the  cost  of  the  method  is  such  that  it  can  only  be  applied  to  very 
small  systems.  It  is  not  size  extensive,  either,  and  its  lack  of  size  extensivity  will  cause 
the  error  in  the  calculation  to  grow  as  the  molecule  grows. 

Gradients  for  Excited  State  Methods 

One  of  the  most  important  advances  in  quantum  chemistry  was  Pulay's  derivation 
of  gradients  for  SCF  wavefunctions.^^  This  allowed  for  the  convenient  computation 
of  the  derivative  of  the  energy  with  respect  to  a  general  perturbation.  By  taking  the 
perturbation  to  be  moving  an  atom,  his  discovery  made  practical  the  exploration  of 
ground  state  potential  energy  surfaces.^^  People  could  readily  find  minima  and  transition 
states  and  calculate  vibrational  frequencies.  Pople  et  al.^'^  extended  gradient  technology 
to  correlated  methods  with  the  development  of  gradients  for  MBPT(2).  For  a  review  of 
current  analytical  gradient  techniques,  see  Ref.  88. 

In  terms  of  excited  states,  gradients  were  first  introduced  for  the  MCSCF  method.^^ 
Although  not  in  general  use,  gradients  are  also  available  for  MRCI  wavefunctions.^^ 
Gradients  for  CIS  were  first  developed  in  1992."^^  Then  came  gradients  for  EOM-CCSD 
in  1993,90-92  followed  by  CIS(D)  gradients.^^  For  a  general  overview  of  gradients  for 
coupled-cluster  based  methods  for  excited  states,  see  Ref.  94. 


CHAPTER  2 

PARTITIONED  EQUATION-OF-MOTION  COUPLED-CLUSTER 

METHODS  FOR  EXCITED  STATE  ENERGIES 

Introduction 

The  equation-of-motion  coupled-cluster  (EOM-CC)  method^^"^^  is  a  conceptually 
single  reference,  generally  applicable,  unambiguous  approach  for  the  description  of 
excited,^^"^^  electron-attached^^  or  ionized  states.^^  All  follow  from  simple  consideration 
of  the  Schrodinger  equation  for  two  states,  a  reference  state  |'^o)  (not  necessarily  the 
ground  state),  and  an  excited  (electron-attached  or  ionized)  state  |^h-).  Considering  H  to 
be  in  second  quantization,  where  the  number  of  particles  is  irrelevant,  the  Schrodinger 
equations  for  the  two  states  are  ., 

Hl^lfo)  =  Eol^o),  (2-1) 

and 

H\^,)^E,\^,).  '  (2-2) 

By  choosing  to  represent  the  excited  state  eigenfunction  as 

|^«)  =  n^l^o),  (2-3) 

it  is  readily  obtained  that 

[H,n^]\^o)  =  Lu,n,\^o)  (2-4) 

for  uj^  =  E^  -  Eq,  from  subtracting  Eq.  (2-1)  from  Eq.  (2-2)  after  left  multiplication 
by  7^«. 
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The  choice  of  TZ^  defines  the  particular  EOM  with  /,  j,  k,  .  .  .  indicating  occupied 
orbital  indices  and  operators,  while  a,  b,  c,  .  .  .  are  unoccupied  orbitals  and  operators. 
Also,  p,  q,  r,  .   .   .   refer  to  orbitals  and  operators  of  either  occupation.   For  electronic 
excited  states, 

^P  =  ^o(ac)  +  ^rf(/.){atr}  +  ^  5]rJ(/c){atifetj|  +  . . . ,  (2_5) 


i.a  '••) 


for  electron  attachment 

-^^^  =  E-'{-){o']  +  \  E  -f  (-){-' Jf^'}  +  •  •  •  >  (2-6) 

"  n,b,j 

and  for  ionization 

T^'/  =  E^'('^)0'}  +  IT.<j(^){^^^^}  +  •  •  ■ '  (2-7) 

t  i,j,a 

Coupled-cluster  theory  is  introduced  by  choosing 

|*0>  =  e^|0)  (2-8) 

where  T  is  the  usual  excitation  operator,  and  |0)  represents  some  independent  particle 
model  reference.  Since  [7^k,  T]  =  0  for  any  of  the  above  choices,  we  can  commute  the 
operators  to  give 

[H,  n,]\0)  =  (^7^,)  JO)  =  a;.7^,|0),      -  -.  (2-9) 

where 

H  =  e-^He'^,  (2-10) 

and  [HTZ^)^  indicates  the  open,  connected  terms  that  remain  in  the  commutator.  Two 
operators  being  connected  means  that  they  share  at  least  one  index.  In  this  way,  all  of 
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the  ground  state  CC  information  is  contained  in  H,  now  generalized  to  have  three-  and 
four-body  terms.^^  Note,  H  is  also  non-Hermitian,  necessitating  that  both  its  left  £«  and 
right  T^K  eigenvectors,  which  form  a  biorthogonal  set,  Cillj  =  6ij,  be  considered  for  a 
treatment  of  properties.    In  matrix  form, 

(2-11) 
Lh-H  =  IjkUJk- 

Consequently,  EOM-CC  reduces  to  a  Cl-like  equation  that  provides  the  relevant  excitation 

energies  directly. 

Besides  the  obvious  approximations,  such  as  T  =  Ti  -f-  T2,  and  11^  being  limited  to 
single  and  double  excitations,  which  defines  EOM-CCSD,  and  various  triple  excitation 
extensions,  E0M-CCSD(T),6'  EOM-CCSDT-l,^',  EOM-CCSD(T),63  and EOM-CCSDT- 
3,^^  one  can  conceive  of  many  other  approximations  to  the  basic  EOM-CC  structure. 
Instead  of  Eq.  (2-8),  for  the  reference  state  a  perturbation  approximation  might  be 
chosen,  such  as 

*('"))  =  (^i  +  r(i)-HV)r(i)  +  T(2)  +  ...)  1 0),    ■       (2-12) 

where  one  truncates  at  a  particular  order,  m.  (See  Ref.  95  for  the  coupled-cluster,  non- 
Hartree-Fock  definition  of  the  various  orders.)  Alternatively,  rather  than  retaining  the 
full  perturbation  approximation  to  |*{,"'^),  we  could  truncate  H  itself  to  some  order.^^-  ^'^ 
It  is  possible  to  conceive  of  perturbative  approximations  to  7^K,  too.  The  latter  are, 
perhaps,  most  easily  viewed  from  the  partitioning  approach  to  perturbation  theory.^^  That 
is,  Eq.  (2-12)  can  be  partitioned  into  the  spaces  P  and  Q,  where  P  represents  the  principal 
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configuration  space  (of  dimension  p)  and  Q  (of  dimension  q)  represents  its  orthogonal 
complement.  Then  it  is  well  known  that  we  can  consider  an  effective  Hamiltonian, 

Spp(a;J  =  Hpp  +  B.pq{uj,  -  Hqq)"^Hqp,  (2-13) 

whose  eigenvectors  are  solely  defined  in  the  P  space 

Hppcp  =  cpu)^  (2-14) 

for  the  first  several  eigenvalues,  w^.  Expanding  the  inverse  in  Eq.  (2-13)  provides  a 
series  of  perturbative  approximations  to  Hpp  or  to  the  eigenvectors  cp  in  Eq.  (2-14).^^ 
With  P  chosen  to  be  the  space  of  single  excitations,  and  Q  that  of  double  excitations,  such 
partitioned  EOM-CCSD  results,  first  presented  in  1989^"^  have  been  shown  to  retain  most 
of  the  accuracy  of  the  full  EOM-CCSD  method.  Other  approximations  can  be  made.  To 
introduce  triple  and  quadruple  excitations,  selective  double  excitations  could  be  retained 
in  P  and  at  least  diagonal  approximations  could  be  made  for  the  triple  and  quadruple 
excitation  blocks  of  Hqq. 

In  this  chapter  such  partitioned  and  perturbation-based  approximations  to  the  full 
EOM-CCSD  method  will  be  reconsidered.  It  will  be  demonstrated  that  very  good 
accuracy  may  be  obtained  within  a  much  less  expensive  computational  structure.  In 
particular,  whereas  the  full  EOM-CCSD  is  proportional  to  nl^^N^-^f  oc  n^,  for  n  basis 
functions,  a  n^^^iV^.^^  ^  ^5  procedure  will  be  developed  that  shows  promise  for  large 
molecules. 
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Theory 

Partitioned  EOM-CC  for  Excited  States 

In  a  typical  EOM-CCSD  calculation,  the  excitation  energy  and  excited  state  properties 
are  calculated  via  diagonalization  of  the  non-symmetric  matrix,^^ 


(H)  = 


Hss      HsD 

Hds    Hdd 


(2-15) 


where  Hss  stands  for  the  singles-singles  block  of  the  matrix,  etc.  The  major  step  in  an 
iterative  diagonalization^^  of  the  matrix  is  multiplying  the  matrix  by  a  trial  vector  C.  In 
the  EE  case,  the  equations  for  this  multiplication  are  as  foUows:^^'  ^ 

[HSSC]^  =  ^  FaeCt  -  X;  FmiCl  ^  ^  H^am,;eC^,  (2-16) 

e  m  em 

[HsdC]?  =  Y,  FmeCZ  +  I  E  ^orncfC^i  "  \  E  ^^^-C^  ,  (2-17) 

''"'  mef  mnc 

[HdsC]  J  =P{ah)  E  W„,,,,Ct  +  P{rj)  E  W,,,,Ct 


+  ^ME     E^We^;)^-^ 

/      \  me 


(2-18) 


n      \  me  / 

[HddC]^  =P{ab)  E  F,,Ct^  -  P{rj)  E  F^^^Ct  +  ^  E  ^^^/^l/ 

+  ^  E  ^rnnijC^Jn  +  P{o^h)P{tj)  E  W,^,,Ct, 
mn  cm 

-lP(-f^)E(E^nmfeC^r:n)i' 
f      \mne  / 

n      \mcf  J 


lae 
'im 
em 


(2-19) 
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The  permutation  operator  P(qr)  is  defined  as 

Piqr)i{.  ..q...r...)=i{...q...r...)-i{...T...q...)  (2-20) 

Fpq  and  Wpqrs  are  presented  elsewhere. ^^ 

In  the  partitioned  scheme,  Hdd  is  replaced  by  Ho^d.  where  Hq  is  the  usual  M0ller- 
Plesset  unperturbed  Hamiltonian  from  many-body  perturbation  theory.'^  In  other  words, 
the  unfolded  matrix  analogous  to  Eq.  (2-13)  to  be  diagonalized  is  approximated  by 

and  Eq.  (2-19)  becomes 

[HoooC]^'  =  P{ab)  J2  heCtf  -  P{ij)  J2  fmjCfl  (2-22) 

e  m 

Here,  fpg  is  the  Fock  operator.  In  the  case  when  the  reference  function  is  composed 
of  canonical  or  semi-canonical  Hartree-Fock  orbitals,  Hqoo  becomes  diagonal  with 
differences  of  orbital  energies  on  the  diagonal.  From  examining  the  equations,  it  is 
clear  that  the  method  is  formally  an  iterative  rv'  method,  and  the  results  are  invariant  to 
rotations  among  occupied  or  unoccupied  orbitals.  This  partitioning  scheme  differs  from 
that  presented  in  Ref.  57  in  two  ways.  Geertsen  et  alr''^  included  the  full  H  elements 
on  the  diagonal  instead  of  just  including  Hq  elements  in  the  doubles-doubles  block,  thus 
destroying  the  orbital  invariance.  Also,  Geertsen  et  al.^'^  did  not  include  the  three  body 
terms  (the  last  two  terms  in  Eq.  (2-18))  in  their  work. 

After  a  matrix  in  the  form  of  Eq.  (2-13)  has  been  diagonalized,  the  eigenvectors 
consisting  of  single  excitations  and  double  excitations  are  built.  Therefore,  the  same 
techniques  used  to  calculate  the  properties  of  a  full  EOM-CCSD  wavefunction^^  can 
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be  used  to  calculate  properties  of  the  partitioned  EOM-CCSD  (i.e.     P-EOM-CCSD) 
wavefunction. 

Since  H  is  not  Hermitian,  its  left  and  right  hand  eigenvalues  are  the  same,  but  the 
eigenvectors  differ.^^  The  equations  for  the  left  hand  side  are  similar  to  those  presented 
above. 

MBPT[2]  Ground  State 

The  first  logical  perturbative  approximation  to 

=  XI  ^viKv^^}  +  XI  ^PVs{vW s'r]  +  higher  order  terms 

P5  pqrs 

is  obtained  by  keeping  only  terms  through  second  order.  Such  an  approximation  defines 
an  E0M-MBPT[2]  method  for  excited  states  (i.e.  an  EOM-CC  calculation  based  on  a 
MBPT[2]  ground  state  instead  of  a  coupled-cluster  ground  state).^^  If  we  insist  upon  a 
method  correct  through  second-order  for  the  non-Hartree-Fock  case,  the  explicit  equations 
for  ^(2)  ^g 

^i  =  0,  -^  (2-24) 

Fah  =  /afc  +  X  C/m6  +  ^  ^Lk'MW  -\Y.  ^"mn{mn\\he) ,  (2-25) 

"I  jm  emn 

Fr,  =  h,  +Y.']h.  +  Y.tl{un\\je)  -  \Y,ti,{M\ef).  (2-26) 

e  cm  eftn 

Fia  =  fia  +  ^t^m  {im\\ae) ,  (2-27) 

em 

Wabij  =  0,  (2-28) 

^ijki  =  mki) + p{ki)Y;^tnij\\ke)  +  ij^'iUm^f)^       (2-29) 

^  ef 


16 

\\r.d\-\- 


WatcJ  =  {ah\\cd)  -  P{ab)J2ti.{am\\cd)  +  lY.*tn{mn\\cd),  (2-30) 

m  mn 

Waibc  =  {ai\\bc)-J^C{mi\\hc),  (2-31) 

m 

Wijka  =  {jk\\ia)  +  J2tl{3k\\ea),  (2-32) 

e 

Wijab  =  {ij\\ah),  ■  (2-33) 

Waijt  =  {ai\\jb)  +  J2tjH\eh)  -^C{m.i\\jh)  -J2Cj{mM\eh),  (2-34) 

e  m  em, 

Wabci  ={ah\\ci)  -  P{ah)Y,tt{^m\\ce)  +  ^  ^C'n(c^1l^n> 

em  mn  ,rs    ^e\ 

+  Y.^tfmc  +  Y.^t{ah\\ce)  -  P{ab)Y,t'L{mb\\ci), 

me  m 

Wiajk  =H\jk)  -  P{^j)J2fZ,{^m\\Je)  +  ^  X]/;{{ia||e/) 

em.  p  f 

,-^  ,r^  (2-36) 

+  T.^jkfie  +  EC(^^lli^>  -  P{ij)'£i]{ia\\ek). 

Cm,  e 

In  all  of  these  equations  tf  and  t^^  refer  to  their  values  through  first  order.    That  is, 

^f^l;  =  -^,  and  i^W"''  =  ^MML_. 

Since  approximating  the  excited  state  with  a  partitioned  EOM-CC  calculation  and  ap- 
proximating the  ground  state  with  a  MBPT(2)  calculation  are  independent  approximations, 
we  can  obviously  combine  them  into  a  partitioned  E0M-MBPT(2)  (P-E0M-MBPT(2)) 
method.  When  combined,  the  H  elements  Wiju,  Wijab,  and  the  most  numerous,  Watcd, 
are  not  needed.  Therefore,  in  the  Hartree-Fock  case,  the  {ab\\cd.)  integrals  never  need  to 
be  calculated.  The  {ab\\cd)  integrals  would  contribute  to  the  Wabd  H  element  multiplied 
by  a  Tj^ ,  but  in  the  Hartree-Fock  case  all  T^'s  are  zero.  A  few  n^  terms  still  remain 
in  the  calculation  of  the  H  elements,  but  these  terms  only  need  to  be  calculated  once. 
In  most  calculations  the  cost  of  the  iterative  n^  step  in  the  excited  state  calculation  will 
dominate  over  the  n^  step  involved  in  calculating  the  H  elements. 
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Calculations 


Be  Atom 


Table  2-1  lists  energies  calculated  using  the  four  methods  described  above  (EOM- 
CCSD,  P-EOM-CCSD,  E0M-MBPT(2),  and  P-E0M-MBPT(2))  for  the  first  few  singlet 
excited  states  in  beryllium.  The  basis  set  from  Ref.  100  was  used  in  the  calculations. 
From  the  mean  absolute  errors,  it  would  appear  that  only  the  EOM-CCSD  calculation  is 
able  to  reproduce  the  full  CI  excitation  energies  (which  are  exact  within  the  basis  set  used). 
However,  if  the  1  ^D  state  is  excluded  from  the  determination  of  the  mean  absolute  error, 
the  errors  become  0.010  eV  for  EOM-CCSD,  0.151  eV  for  P-EOM-CCSD,  0.431  eV  for 
E0M-MBPT(2),  and  0.300  eV  for  P-E0M-MBPT(2).  With  an  approximate  excitation 
level  (AEL)  value  of  1.60,  the  l^D  state  is  dominated  by  double  excitation  character. 
The  AEL  is  a  measure  of  the  number  of  electrons  excited  in  an  excitation.^^  A  value  of 
1.00  is  a  pure  single  excitation,  while  2.00  corresponds  to  a  pure  double  excitation.  Since 
the  doubles-doubles  block  of  the  EOM  wavefunction  is  approximated  in  this  partitioned 
scheme,  it  is  reasonable  to  expect  any  state  with  appreciable  double  excitation  character 
will  be  poorly  described  by  this  partitioned  EOM  calculation. 

Even  with  the  I'D  excluded,  the  calculations  based  on  a  MBPT(2)  ground  state  are 
still  very  poor.  This  poor  behavior  can  be  ascribed  to  the  inadequacy  of  the  MBPT(2) 
wavefunction  in  describing  the  ground  state  of  Be.  In  a  CCSD  calculation  the  largest 
T2  amplitude  is  0.065,  while  the  largest  TJ^'  amplitude  in  the  MBPT(2)  calculation  is 
only  0.034.  Also,  the  MBPT(2)  calculation  only  obtains  73%  of  the  correlation  energy 
recovered  by  the  CCSD  calculation.  Clearly,  MBPT(2)  is  not  adequate  to  describe  the 
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ground  state,  and  any  excited  state  calculation  based  on  that  MBPT(2)  ground  state  will 
suffer  accordingly. 

In  Table  2-2  the  first  few  triplet  excited  states  of  beryllium  starting  from  the  singlet 
ground  state  are  presented.  Since  all  of  the  calculated  states  are  single  excitations,  the 
EOM-CCSD  method  does  exceedingly  well,  while  the  P-EOM-CCSD  also  does  a  good 
job  of  describing  the  states. 

Example  Molecules 

While  comparisons  with  full  CI  calculations  are  useful,  it  is  also  informative  to  look 
at  the  performance  of  the  methods  for  more  common  molecules.  These  calculations 
will  also  provide  an  opportunity  to  compare  the  current  methods  with  other  single 
reference  methods  used  today.  Tables  2-3  to  2-6  present  calculations  on  four  molecules: 
formaldehyde,  acetaldehyde,  ethylene,  and  butadiene.  The  calculations  are  performed  at 
the  MP2/6-31G*  geometries  given  in  Refs.  101-103.  A  6-31  l(2+,2+)G**  basis  set'^' 
is  used  for  formaldehyde  and  ethylene,  with  a  6-311(2+)G*  basis  set'^^  being  used  for 
acetaldehyde  and  butadiene.  All  electrons  are  correlated  except  in  butadiene,  where  the 
first  four  core  orbitals  are  dropped. 

All  of  the  methods  presented,  except  for  the  CIS-MP2,^^  can  be  viewed  as  approxi- 
mations to  the  full  EOM-CCSD  method.  TDA  (or  CIS)'°^  is  the  crudest  of  the  methods  in 
that  the  excited  state  is  given  as  a  linear  combination  of  single  excitations  out  of  a  single 
reference  ground  state.  This  method  includes  no  dynamic  correlation.  CIS(D)^^  provides 
a  non-iterative  n^  perturbation  correction  to  the  CIS  energy.  CIS-MP2  also  provides  a 
correction  to  the  CIS  energy,  but  the  method  is  not  size-consistent  and  scales  as  n'^P 
The  partitioned  methods  provide  an  iterative  n^  excited  state  method  based  on  either  a 
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n^  ground  state  (P-E0M-MBPT(2))  or  an  iterative  n*^  ground  state  (P-EOM-CCSD).  The 
full  EOM-CCSD  method  is  an  iterative  n^  method. 

Except  for  two  differences,  our  assignment  of  states  agrees  with  the  assignments 
of  Wiberg  et  al.  and  with  the  corrected  assignments  of  Head-Gordon  et  a/.,'°^.  The 
states  have  been  ordered  based  on  their  experimental  excitation  energies,  or,  where  not 
available,  their  EOM-CCSD  excitation  energies,  instead  of  their  CIS  excitation  energies. 
Also,  Head-Gordon  etalP  had  incorrectly  assigned  the  'A]  EOM-CCSD  state  at  9.27  eV 
in  formaldehyde  to  the  valence  4'Ai  state.  Based  upon  the  state's  properties,  including 
an  <r^>  value  of  54.4  au^  (compared  to  the  ground  state  <r^>  value  of  20.5  au^),  the 
'Ai  EOM-CCSD  state  has  been  reassigned  to  the  Rydberg  3'A|  state.  The  4'Ai  state 
(the  7r*<— TT  state),  with  an  <r^>  value  of  32.6  au^,  has  an  EOM-CCSD  excitation  energy 
of  10.00  eV. 

Considering  the  basis  sets  and  the  relatively  inexpensive  methods  used  here,  these 
calculations  are  not  meant  to  be  definitive.  For  previous  results  on  these  molecules 
see  the  references  in  Refs.  101-103.  Instead,  the  mean  absolute  errors  of  the  methods 
compared  to  experiment  and  compared  to  EOM-CCSD  will  be  used  to  access  the  quality 
of  the  methods.  The  mean  absolute  error  for  each  of  the  methods  with  respect  to  the 
experimental  values  given  is  0.66  eV  for  CIS,  0.41  eV  for  CIS-MP2,  0.32  eV  for  CIS(D), 
0. 17  eV  for  P-E0M-MBPT(2),  0.20  eV  for  P-EOM-CCSD,  and  0. 14  eV  for  the  full  EOM- 
CCSD  method.  Compared  to  the  more  complete  EOM-CCSD  method,  the  mean  absolute 
errors  of  the  various  methods  are  0.70  eV  for  CIS,  0.39  eV  for  CIS-MP2,  0.35  for  CIS(D), 
0.12  eV  for  P-E0M-MBPT(2),  and  0.16  eV  for  P-EOM-CCSD. 
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As  should  be  expected  from  its  simple  nature,  TDA  performs  poorly.  The  average 
errors  for  CIS-MP2  and  CIS(D)  are  similar,  although,  as  previously  noted,^^  the  CIS-MP2 
energies  are  more  erratic.  The  P-E0M-MBPT(2)  and  P-EOM-CCSD  energies  are  also 
similar,  seldom  differing  by  0.1  eV.  This  suggests  that  for  these  cases  where  MBPT(2)  is 
able  to  well  represent  the  ground  state,  P-E0M-MBPT(2)  should  be  able  to  well  describe 
singly  excited  states.  A  balance  argument  between  the  ground  and  excited  state  would 
also  tend  to  favor  the  P-E0M-MBPT(2)  method,  since  the  partitioning,  as  discussed 
above,  can  be  viewed  as  a  second-order  in  H  perturbation  expansion  for  the  excited  state. 

Looking  only  at  the  valence  states,  as  designated  by  Wiberg  et  al.,^^^~^^'^  the  mean 
absolute  errors  from  the  EOM-CCSD  energies  are  0.45  eV  for  CIS,  0.68  eV  for  CIS-MP2, 
0.33  eV  for  CIS(D),  0.23  for  P-E0M-MBPT(2),  and  0.38  eV  for  P-EOM-CCSD.  This  is 
compared  to  the  Rydberg  states,  where  the  mean  absolute  errors  from  the  EOM-CCSD 
energies  are  0.76  eV  for  CIS,  0.32  eV  for  CIS-MP2,  0.35  eV  for  CIS(D),  0.10  for  P- 
E0M-MBPT(2),  and  0.11  eV  for  P-EOM-CCSD.  It  appears  that  the  partitioned  methods 
on  average  do  not  describe  valence  states  as  well  as  they  describe  Rydberg  states.  Since 
orbital  relaxation  is  often  greater  for  valence  states,  it  is  possible  the  reduced  relaxation 
available  by  approximating  the  doubles  vector  might  cause  larger  errors  for  valence  states. 

These  calculations  also  give  a  good  chance  to  measure  the  time  savings  of  the 
partitioned  method.  For  example,  for  butadiene  the  full  EOM-CCSD  calculation  took 
25905  seconds  on  an  IBM  RS/6000  model  590  to  calculate  eighteen  excited  states  (both 
right  and  left  hand  sides).  To  calculate  the  same  states  with  the  P-EOM-CCSD  method 
took  3998  seconds.  These  times  only  include  the  time  for  the  excited  states.   P-EOM- 
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MBPT(2)  also  saves  time  in  the  calculation  of  the  ground  state  and  in  the  formation  of 
H,  as  well  as  saves  disk  space. 

Conclusions 

In  this  chapter  a  formalism  is  presented  and  results  are  given  for  partitioned  methods 
based  on  the  equation-of-motion  coupled-cluster  theory,  where  the  ground  state  can  be 
described  by  either  a  CCSD  or  an  MBPT(2)  wavefunction.  The  partitioned  methods 
provide  an  iterative  n^  method  (plus  an  n^  step  for  forming  H  elements)  for  excited  states. 
When  the  ground  state  of  the  system  is  well  described  by  an  MBPT(2)  wavefunction, 
the  P-EOM-MBPT(2)  method  provides  an  inexpensive  way  to  accurately  calculate  the 
energies  and  properties  of  singly  excited  states.  For  systems  less  well  described  by 
a  MBPT(2)  wavefunction,  the  P-EOM-CCSD  method  is  a  generally  accurate,  but  more 
economical,  alternative  to  a  full  EOM-CCSD  calculation.  For  vertical  excitation  energies, 
P-E0M-MBPT(2)  is  a  superior  n^  method  to  CIS(D). 


'I 
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Table  2-1:  Be  singlet  excitation  energies  (in  eV) 


EOM- 

P-EOM- 

EOM- 

P-EOM- 

state 

AEL^ 

CCSD 

CCSD'' 

MBPT(2) 

MBPT(2) 

full  CF 

I'P 

1.07 

5.323 

5.666 

4.869 

5.228 

5.318 

2iS 

1.06 

6.773 

6.985 

6.329 

6.576 

6.765 

l^D 

1.60 

7.139 

10.579 

6.682 

10.170 

7.089 

2ip 

1.06 

7.468 

7.653 

7.023 

7.224 

7.462 

2'D 

1.21 

8.055 

7.892 

7.618 

7.477 

8.034 

s's 

1.04 

8.084 

8.203 

7.648 

7.788 

8.076 

3ip 

1.05 

8.309 

8.432 

7.873 

8.011 

8.302 

3^0 

1.09 

8.548 

8.548 

8.115 

8.132 

8.536 

4'S 

1.04 

8.583^ 

8.694 

8.278 

8.600 

4'? 

1.04 

8.700 

8.792 

8.267 

8.374 

8.693 

mean 

0.014 

0.485 

0.428 

0.578 

abs.  error 

(0.01  Of 

(0.151f 

(0.43  If 

(0.300)^ 

^^Ref.  60.  The  AEL  is  for  the  EOM-CCSD  wavefunction. 

''^Because  of  a  different  implementation  of  the  partitioning,  these  numbers  are  slightly 
different  than  those  in  Ref.  57. 
•^^Ref.  100. 
''^Ref.  60. 
^^ Average  error  without  I'D  double  excited  state. 
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Table  2-2:  Be  triplet  excitation  energies  (in  eV) 


EOM- 

P-EOM- 

EOM- 

P-EOM- 

state 

AEL^ 

CCSD 

CCSD'' 

MBPT(2) 

MBPT(2) 

full  CF 

l3p 

1.01 

2.729 

2.819 

2.271 

2.366 

2.733 

2^8 

1.04 

6.447 

6.583 

5.997 

6.140 

6.444 

23p 

1.04 

7.301 

7.424 

6.868 

7.006 

7.295 

l^D 

1.04 

7.748 

7.866 

7.316 

7.448 

7.741 

a^s 

1.04 

7.991 

8.089 

7.557 

7.668 

7.985 

43p 

1.03 

8.278 

8.366 

7.949 

8.272 

23d 

1.03 

8.456 

8.539 

8.025 

8.121 

8.449 

43s 

8.58^^ 

8.647 

8.228 

8.560 

53p 

8.70^^ 

8.766 

8.348 

8.686 

mean 

abs.  error 

0.008 

0.104 

0.436 

0.321 

^)The 

AEL 

is  for  the  EOM-CCSD  wavefunction. 

'')See 

footnote  b. 

Table  2-1. 

^)Ref. 

100. 

''^Ref. 

59. 
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Table  2- 

-3:  Formaldehyde  (energies  in  eV) 

CIS- 

P-EOM- 

P-EOM- 

EOM- 

CIS^ 

MP2'^ 

CIS(D)b 

MBPT(2)^ 

CCSD^ 

CCSD 

Exp.^ 

l'A2(V) 

4.48 

4.58 

3.98 

4.31 

4.41 

3.95'' 

4.07 

l'B2(R) 

8.63 

6.85 

6.44 

7.10 

7.10 

7.06'' 

7.11 

2'B2(R) 

9.36 

7.66 

7.26 

7.97 

7.95 

7.89'' 

7.97 

2'Ai(R) 

9.66 

8.47 

8.12 

8.02 

8.01 

8.00'' 

8.14 

2'A2(R) 

9.78 

7.83 

7.50 

8.25 

8.23 

8.23'' 

8.37 

3'B2(R) 

10.61 

8.46 

8.21 

8.99 

8.97 

9.07'' 

8.88 

l'Bi(V) 

9.66 

9.97 

9.37 

9.61 

9.68 

9.26" 

3>Ai(R) 

10.88 

8.75 

8.52 

9.24 

9.21 

9.27'= 

4'B2(R) 

10.86 

8.94 

8.63 

9.39 

9.38 

9.40'' 

41  A,  (V) 

9.45 

9.19 

8.80 

10.08 

10.24 

10.00<= 

'')Ref. 

102. 

'')Ref. 

83. 

'^^Present  work. 
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Table  2-4 

:  Acetaldehyde  (energies  in  eV) 

» 

els'* 

CIS- 

MP2^ 

CIS(D)'' 

4.28 

P-EOM- 

MBPT(2)'= 

4.65 

P-EOM- 
CCSD'^ 

4.71 

EOM- 
CCSD 

4.26" 

Exp.^ 
4.28 

l'A"(V) 

4.89 

5.27 

2iA'(R) 

8.51 

6.71 

6.13 

6.88 

6.84 

6.78'' 

6.82 

3*A'(R) 

9.22 

7.57 

7.04 

7.57 

7.52 

7.49*^ 

7.46 

2'A"(R) 

9.37 

7.37 

6.90 

7.70 

7.64 

7.64'' 

4'A'(R) 

9.30 

8.00 

7.42 

7.77 

7.70 

7.68'' 

7.75 

.1 

< 

i 

S'A'CR) 

10.19 

8.09 

7.70 

8.42 

8.35 

8.39'' 

8.43 

\ 

6'A'(R) 

10.26 

8.08 

7.70 

8.53 

8.47 

8.51'' 

8.69 

3^A"(R) 

10.31 

8.10 

7.74 

8.58 

8.51 

8.57''  . 

4iA"(V) 

9.78 

10.34 

9.34 

9.61 

9.65 

9.23*= 

7'A'(V) 

9.73 
102. 

9.07 

8.50 

9.58 

10.07 

9.44'= 

■-f'. 

l^ 

''^Ref. 

.      ' 

''^Ref. 

83 

-■■.:-. m  ■ 

■  -a 

1 

'^^Present 

work. 

.'•',■ 
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Table  2-5:  Ethylene  (energies  in  eV) 


CIS- 

P-EOM- 

P-EOM- 

EOM- 

CIS^ 

MP2^ 

CIS(D)*' 

MBPT(2)'^ 

CCSD^ 

CCSD 

7.31" 

Exp.^ 

l'B3u(R) 

7.13 

7.52 

7.21 

7.45 

7.51 

7.11 

l'B,u(V) 

7.74 

8.39 

8.04 

8.20 

8.36 

8.14" 

7.60 

l'B,g(R) 

7.71 

8.14 

7.84 

8.08 

8.14 

7.96" 

7.80 

l'B2g(R) 

7.86 

8.12 

7.86 

8.12 

8.18 

7.99" 

8.01 

2'Ag(R) 

8.09 

8.42 

8.18 

8.43 

8.48 

8.34" 

8.29 

2'B3u(R) 

8.63 

8.92 

8.69 

8.95 

9.00 

8.86" 

8.62 

l'Au(R) 

8.77 

9.00 

8.80 

9.07 

9.12 

9.01" 

3'B3u(R) 

8.93 

9.14 

8.96 

9.23 

9.28 

9.18" 

2'B,u(R) 

9.09 

9.38 

9.18 

9.42 

9.51 

9.39*^ 

9.33 

2'Big(R) 

9.09 

9.31 

9.12 

9.42 

9.48 

9.38^^ 

9.34 

'^^Ref.  101. 

''^Ref.  83 

''^Present 

work. 
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Table  2-6:  Butadiene  (energies  in  eV) 


CIS- 

P-EOM- 

P-EOM- 

EOM- 

CIS'' 

MP23 

CIS(D)b 

MBPT(2)'^ 

CCSD'^ 

CCSD 

Exp." 

l'Bu(V) 

6.21 

7.00 

6.29 

6.52 

6.63 

6.42'' 

5.91 

l'Bg(R) 

6.11 

6.73 

6.11 

6.40 

6.39 

6.20'' 

6.22 

l^Au(R) 

6.45 

7.03 

6.44 

6.72 

6.73 

6.53'' 

2'Au(R) 

6.61 

7.11 

6.55 

6.85 

6.84 

6.67'' 

6.66 

2iBu(R) 

6.99 

7.58 

7.03 

7.29 

7.34 

7.17'' 

7.07 

2iBg(R) 

7.22 

7.66 

7.17 

7.47 

7.47 

7.31'' 

7.36 

2'Ag(R) 

7.19 

7.74 

7.19 

7.43 

7.44 

7.10'' 

7.4 

B'BgCR) 

7.25 

7.74 

7.24 

7.54 

7.53 

7.39'' 

7.62 

4'Bg(R) 

7.39 

7.87 

7.40 

7.70 

7.69 

7.55'' 

7.72 

3iAg(R) 

7.45 

7.88 

7.44 

7.74 

7.73 

7.61'' 

S'BuCR) 

8.05 

8.40 

8.01 

8.32 

8.31 

8.21'= 

8.00 

S'AuCR) 

7.78 

6.75 

7.73 

8.07 

8.07 

7.92'= 

8.18 

4'Au(R) 

7.92 

7.66 

7.86 

8.19 

8.18 

8.06^ 

8.21 

«>Ref. 

103. 

^'^Ref. 

83 

'^^Present 

work. 

CHAPTER  3 

GRADIENTS  FOR  THE  PARTITIONED  EQUATION-OF- 

MOTION  MBPT(2)  METHOD  FOR  EXCITED  STATES 

Introduction 

The  equation-of-motion  coupled-cluster  method  with  single  and  double  excitations 
for  excitation  energies  (EOM-CCSD)^^-  ^^  provides  an  accurate  method  for  calculating 
the  energy  and  properties  of  many  excited  states  of  molecules.  With  the  development  of 
gradients  for  EOM-CCSD^^^^  it  has  become  possible  to  study  the  excited  state  potential 
energy  surfaces,  just  as  the  development  of  coupled-cluster  gradients  for  ground  state 
methods'°^'°^  made  routine  the  study  of  ground  state  potential  energy  surfaces.  More 
recently,  Stanton  and  Gauss  have  developed  second  derivatives  for  E0M-CCSD'^°  and 
gradients  for  ionization  potential  EOM-CCSD  (IP-EOM-CCSD).^''  The  key  to  all  coupled- 
cluster  gradients  has  been  the  recognition  of  how  the  interchange  theorem  of  Dalgarno 
and  Stewart,'^'  initially  used  by  Adamowicz  et  aO^^  and  Bartlett'""^  for  coupled-cluster 
gradients  and  sometimes  known  in  quantum  chemistry  as  the  Z-vector  method,' '^  could 
be  used  to  avoid  computing  the  response  of  the  ground  state  T  amplitudes  to  each 
perturbation. 

One  of  the  primary  problems  with  the  EOM-CCSD  method  is  its  cost.  An  EOM- 
CCSD  gradient  calculation  involves  four  iterative  steps  which  scale  as  nl^^N^-^^,  where 
nocc  stands  for  the  number  of  orbitals  occupied  in  the  reference  determinant,  and  7Vv,>, 
stands  for  the  number  of  orbitals  unoccupied  in  the  reference  determinant  .  These 
iterative  nLc^^tX  steps  are  calculating  the  ground  state  T  amplitudes,  calculating  the 
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right  hand  eigenvector  TZ,  calculating  the  left  hand  eigenvector  C,  and  inverting  H  as 
part  of  calculating  the  Z  vector. 

In  the  philosophy  of  replacing  ground  state  cluster  amplitudes  with  second-order 
perturbation  approaches,  "•'•  ^^  Stanton  and  Gauss^^  developed  an  EOM-CC  method  based 
on  a  ^  that  was  expanded  through  second  order.  For  calculations  based  on  a  Hartree- 
Fock  reference,  expanding  H  through  second  order  and  replacing  the  CCSD  amplitudes 
with  the  MBPT(2)  amplitudes  in  the  EOM-CC  equations  are  functionally  equivalent. 
This  truncation  of  H  reduced  the  cost  of  calculating  the  T  amplitudes  to  an  non-iterative 
nocc{nocc  +  Nyirf)  step  and  reduced  the  cost  of  inverting  H  to  ri^occNlirt^''  However, 
the  iterative  n^cc^ttv^  ^^^P^  for  calculating  11  and  C  still  remained.  Since  this  did  not 
significantly  reduce  the  cost,  the  method  was  not  very  useful  in  practice. 

In  order  to  further  reduce  the  cost  and  to  create  a  more  practical  method  geared 
toward  providing  a  robust  second-order-like  treatment  of  excited  states,  in  analogy  with 
the  robust  MBPT[2]  for  ground  states,  "'*  we  proposed  partitioning  away  the  doubles- 
doubles  block  of  H.  It  is  in  the  doubles-doubles  block  of  H  that  all  of  the  iterative 
'^Icc^virt  Steps  in  calculating  7^  and  C  arise.  In  practice  this  means  replacing  H  in 
the  doubles-doubles  block  with  Hq,  which  is  diagonal  in  the  Hartree-Fock  case.  Thus, 
the  iterative  nl^^N^-^^  steps  involved  in  calculating  U  and  C  are  replaced  with  iterative 
^IccMirt  Steps.  The  overall  effect  is  that  the  cost  is  significantly  reduced,  with  little  loss 
in  accuracy  in  vertical  excitation  energies."'^ 

Now,  however,  the  non-iterative  nl^^N^^^^  step  involved  in  calculating  H  elements 
can  also  become  significant.  The  point  at  which  the  cost  of  calculating  the  one  'n?  iV*.  . 
step  in  H  becomes  dominant  over  the  cost  of  calculating  the  three  iterative  nl^^N^-^f  steps 
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in  the  excited  state  can  be  estimated.  The  two  will  have  approximately  equal  cost  when 
the  number  of  virtual  orbitals  is  three  times  the  number  of  iterations.  Since  it  typically 
takes  around  fifteen  iterations  to  solve  for  each  TZ  and  each  £,  this  would  imply  that 
the  cost  of  calculating  H  would  become  significant  when  the  number  of  virtual  orbitals 
is  more  than  about  ninety  per  excited  state  being  calculated.  In  practice,  some  of  the 
terms  not  included  in  this  estimate  also  contribute  significantly  to  the  cost,  increasing 
the  break-even  point  substantially.  Also,  when  calculating  vertical  excitation  energies, 
multiple  excited  states  can  be  calculated  while  forming  H  only  once. 

Theory 
EOM-CC  Gradients 

Before  we  beginning  with  an  overview  of  EOM-CC  theory  and  EOM-CC  gradients, 
it  will  be  helpful  to  define  some  notation.  Given  a  set  of  spin  orbitals  and  a  reference 
determinant  |0),  the  labels  /,;,  k,  .  .  .  represent  spin  orbitals  occupied  in  |0),  the  labels  a, 
b,  c,  .  .  .  represent  spin  orbitals  unoccupied  in  |0),  and  the  labels  p,  q,  r,  .  .  .  represent 
spin  orbitals  whose  occupation  is  not  specified.  The  reference  determinant  is  usually, 
but  not  necessarily,  the  SCF  determinant.  The  space  of  all  possible  determinants  with 
n  electrons  formed  from  the  spin  orbitals  is  represented  by  |h).  The  space  |h)  is  then 
divided  into  |h)  =  |p)  0  |q).  Here,  |p)  represents  the  space  spanned  by  the  operators 
C  and  TZ,  where  it  is  assumed  that  T,  C,  and  7^  all  have  the  same  rank.  Finally,  |p)  is 
divided  into  |p)  =  |0)  ®  |g).  For  EOM-CCSD  |g)  consists  of  all  determinants  singly  and 
doubly  excited  with  respect  to  |0).  This  is  the  same  notation  used  by  Stanton.^° 

The  ground  state  coupled-cluster  equations  can  be  written  as 

{0\H\0)  =  Ecc,  (3-1) 


.V  :'-- 
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and 

(gW)  =  0.        •  (3-2) 

In  pure  CC  methods  (e.g.  CCSD  or  CCSDT),  H  can  be  written  as 

'.V  '  '  '  '  "  '■'■'.'■> 

.      >.  ,    ,  F  =  e-^//e^  =  (/ffi^)  ,  .^  (3-3) 

where  the  subscript  c  indicates  that  the  terms  are  connected;  they  share  at  least  one  index 


m  common. 


T  =  Y.'t{a^^}  +  \E'^l^{o■Ub^J}  +  ...  (3-1) 


i,a 


is  an  excitation  operator  which  accounts  for  the  ground  state  correlation.  In  the  iterative 
triples  methods  such  as  the  CCSDT-n  methods,"-^  or  in  many-body  perturbation  theory, 
H  does  not  have  such  a  concise  form.  Therefore,  Eqs.  (3-1)  and  (3-2)  will  be  taken 
as  the  defining  equations  of  H. 

It  is  now  possible  to  write  the  EOM-CC  equations^^  as 

(0|£^|p)  =  £;(0|£|p),  (3-5) 

and 

{p\Hn\Q)  ^  E{p\n\0).  (3-6) 

Here,  E  is  the  total  energy  for  the  excited  state, 

C  =  Io  +  Y:  la{i^-}  +  I  E  Ci{':^-i^^}  +  •  •  ■  (3-7) 

a,i  ",6 

'J 

is  a  left  eigenvector  of  H,  and 

n  =  ro  +  E^°{«^4  +  IT.<'{^^'^^^}  +  ■  •  •  (3-8) 
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is  a  right  eigenvector  of  H.  It  is  also  possible  to  include  the  equation-of-motion  in  the 
Tl  equation.    Eq.  (3-6)  then  becomes^^ 

(p|[^,7e]|0)=a;(p|7e|0),  (3-9) 

where 

uj^E-Ecc  (3-10) 

is  the  excitation  energy. 

It  can  be  seen  that  in  EOM-CC  the  excited  states  are  the  eigenvectors  of  H,  with 
the  corresponding  eigenvalues  being  the  energies.  Since  H  is  formed  from  a  similarity 
transformation,  it  is  not  Hermitian,  but  its  eigenvectors  do  form  a  biorthogonal  set. 
Choosing  the  norm  of  the  excited  state  vectors  to  be  one  gives  that 

m,nu\^)^8^k  (3-11) 

for  all  states  j  and  k.  Therefore,  from  either  Eq.  (3-5)  or  (3-6)  we  can  get  an  expression 
for  the  excited  state  energy 

E  =  (S)\CH'R\Q).  .      ^      ■  (3-12) 

When  Stanton  first  derived  the  equations  for  EOM-CC  gradients,^°  he  followed  the 
same  development  as  that  of  Refs.  106,  107,  109  by  taking  the  derivative  of  Eq.  (3-12) 
with  respect  to  a  general  external  perturbation.  In  the  process,  he  had  to  introduce  the 
perturbation  independent  quantity  zeta  (Z)  to  account  for  the  effect  of  the  excited  state 
perturbation  on  the  ground  state  T  amplitudes.  Since  C  and  R  are  eigenvectors  for  the 
excited  state,  they  are  variationally  optimum  and  therefore  stationary  with  respect  to  the 
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perturbation.  On  the  other  hand,  T  is  never  stationary  (even  for  the  ground  state),  and 
another  operator  must  be  included  to  account  for  its  response. 

Szalay'"^  then  used  the  approach  of  a  general  functional  analogous  to  the  A  functional 
in  ground  state  coupled-cluster  theory''^-  ^'^  by  introducing  Z  to  force  T  to  be  stationary 
with  respect  to  the  external  perturbation.  Following  this  approach,  here  is  introduced 
an  EOM-CC  functional  which  has  the  property  that  all  quantities  are  stationary,  and 
therefore  the  functional  will  satisfy  the  generalized  Hellmann-Feynman  theorem."^-  "^ 
Such  a  functional  is 

F  =  (0|£^7^|0)  +  (0|Z^|0)  +  E[l  -  (0|£7e|0)).  (3-13) 

By  construction 

2  =  EcH  +  iECiOV/'}  +  ...  (3-14) 

i,i 

is  a  pure  deexcitation  operator  of  rank  equal  to  T.  Since  the  last  two  terms  of  Eq.  (3-13) 
do  not  contribute  to  the  value  of  the  functional,  the  value  of  the  functional  is  just  the 
energy.  By  taking  the  derivative  of  each  of  the  quantities  on  the  right  hand  side  of 
Eq.  (3-13)  and  setting  them  equal  to  zero,  the  EOM-CC  gradient  equations  will  be 
recovered.    For  example, 

dF 

—  =  0  =  1-  (0|£7^|0)  (3-15) 

is  a  restatement  of  Eq.  (3-11),  the  normalization  condition  for  £  and  R.  Taking  the 
derivative  of  the  functional  with  respect  to  R  (C)  gives  Eq.  (3-5)  (Eq.  (3-6)),  which  is 
the  equation  for  C  (R).  Taking  the  derivative  of  the  functional  with  respect  to  Z  gives 
Eq.  (3-2),  which  is  the  ground  state  T  equation. 


"IV"*    ' 


34 
However,  taking  the  derivative  of  the  functional  with  respect  to  T  is  much  more 
complicated.  Since  the  form  of  H,  and  therefore  the  equations  for  T,  vary  between 
different  CC  methods,  it  should  not  be  surprising  that  the  equations  for  Z  also  vary 
between  different  EOM-CC  methods.  Here  will  be  assumed  a  pure  EOM-CC  method 
(e.g.  EOM-CCSD,  EOM-CCSDT,  ...).  When  the  derivative  with  respect  to  one  of  the 
operators  is  taken,  C,  R,  T,  or  Z,  what  is  really  being  done  is  separately  taking  the 
derivative  with  respect  to  every  coefficient  in  the  operator.  For  example,  ^  really 
means  ^  for  all  t  amplitudes  in  T.  This  still  leaves  the  excitation  or  deexcitation  part  of 
the  operator.  So  when  |f^  is  taken  as  a  part  of  getting  |^,  what  results  is 

:      (3-16) 


—  =  -ftre-^He^  +  e-^He^^T 


Here 

is  the  pure  excitation  operator  without  any  coefficients.  Therefore, 

dF 

^  =  0=(0|£[^,nr]7e|0)  +  (0|Z[F,Or]|0).  (3-18) 

After  significant  simplification,  Eq.  (3-18)  reduces  to 

■      :/  ..        (0|Z|g)  =  (0|£^|q)(q|7e|g)((g|(i?Ccl-^)|g))"\  (3-19) 


where  1  stands  for  a  unit  matrix  of  rank  |g).  This  is  exactly  equivalent  to  the  equation 
for  Z  in  Ref.  90. 

Finally,  since  the  functional  satisfies  the  generalized  Hellmann-Feynman  theorem, 
for  any  perturbation  x 

^  =  ^  =  (0|^^"^|0)  +  (0|Z^>^|0).  (3-20) 
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Here,  H^  represents  the  derivative  of  the  bare  Hamiltonian  elements  with  respect  to  the 
perturbation  x  within  the  similarity  transformed  Hamiltonian  H.  For  the  pure  EOM-CC 
methods  this  means 

m  =  e-^^e^.  (3-21) 

Eq.  (3-20)  is  the  same  as  Eq.  (36)  of  Ref.  90.  Although  the  form  of  the  functional  was 
simply  postulated,  the  fact  that  it  produces  equations  (and  therefore  density  matrices) 
identical  to  those  produced  through  straightforward  derivation  of  the  energy  expression^° 
proves  that  the  functional  is  valid.  For  a  discussion  of  the  properties  of  the  density 
matrices  resulting  from  these  equations,  see  Ref.  119. 

It  is  also  possible  to  derive  an  energy  expression  from  Eq.  (3-9).  Projecting  on  the 
left  by  (0|£|p)  gives  the  expression 

E={0\C[H,n]\0)  +  {0\H\0),  (3-22) 

where  the  second  term,  Eqc,  is  added  in  to  give  a  total  energy  instead  of  just  the 
excitation  energy.  The  functional  which  corresponds  to  this  energy  expression  is 

F  =  {0\ClH.n]\0)  +  {0\H\0)  +  (0|Z^|0)  +a;(l  -  (0|£7^|0)).  (3-23) 

Superficially,  Eqs.  (3-13)  and  (3-23)  appear  similar.  In  fact,  Eq.  (3-22)  can  be  straight- 
forwardly derived  from  Eq.  (3-12),  and  both  functional  have  the  same  value,  the  total 
energy  E.  However,  these  energy  functionals  contain  several  important  differences.  The 
first  is  that  uj,  the  excitation  energy,  appears  on  the  right  instead  of  E.  More  importantly, 
though,  the  equations  for  and  value  of  Z  have  changed.  This  will  seen  later,  but  first  it 
is  important  to  see  that  the  EOM-CC  equations  are  still  contained  in  the  new  functional. 
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As  before,  the  derivative  of  each  of  the  quantities  on  the  right  hand  side  of  the 
functional  will  be  taken  and  set  to  zero.   First, 

^  =  0  =  l-(0|£7^|0)  (3-24) 

again  gives  the  normalization  condition  for  £  and  R. 

OF 

—  =  0  =  (g|if  |0)  (3-25) 

again  gives  the  coupled-cluster  equations. 

Next,  taking  the  derivative  of  the  functional  with  respect  to  C  gives 

dF 

^  =  0=(p|[F,7^]|0)-a;(p|7^|0),     ;  (3-26) 

which  is  equivalent  to  Eq.  (3-9).  Expanding  the  commutator  gives 

a;(p|7l|0)  =  (p|^7e|0)  -  (p|7e^|0).  (3-27) 

Inserting  a  resolution  of  the  identity  |h)(h|  into  the  last  term  gives 

u;(p|7^|0)  =  (p|^7^|0)  -  (p|7^|q)(q|^|0) 

-  (Pl^|g)(g|^|0)  -  (p|7e|0)(0|^|0)  (3-28) 

=  (p|^7e|o)  -  (p|7^|o)Ecc, 

or 

{Ecc  +  a;)(p|7^|0)  =  {p\Hn\Q).  (3-29) 

The  third  term  in  the  first  equality  of  Eq.  (3-28)  disappears  because  of  the  solution  of  the 
coupled-cluster  equations.  The  second  term  in  the  first  equality  of  Eq.  (3-28)  disappears 
because  R,  being  an  excitation  operator,  can  never  decrease  the  excitation  level.  Thus 
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it  cannot  have  any  non-vanishing  terms  between  (p|  and  |q).  Note  that  Eq.  (3-29)  is 
equivalent  to  Eq.  (3-6). 

Taking  the  derivative  of  the  functional  with  respect  to  R  gives 
OF 

-^  =  Q  =  {Q\c[H,nTz]\o)-uj{o\mn\o).     .:  (3-30) 

In  analogy  to  Qt  above,  Q,ti  is  the  excitation  part  of  the  operator  R  without  any 
coefficients.  The  difference  is  that  Qt  acting  on  |0)  gives  |g),  while  Ot^  acting  on  |0) 
gives  |p),  since  Q.-ji  contains  a  constant  part  and  ftj  does  not.  Rewriting  Eq.  (3-30)  gives 

c^(0|£|p)  =  (OlCHnnlO)  -  {0\CnnH\0).  (3-31) 

Inserting  a  resolution  of  the  identity  into  the  last  term  gives 

a;(0|£|p)  =  (0|£^|p)  -  (0|£1^7l|0)(0|^|0>  -  {0\Cnn\p){p\H\0) 

-(0|£l]7e|q)(q|//|0)  (3-32) 

=  {0\CH\p)  -  {0\C\p)Ecc, 


or 


(a;  +  £;cc)<0|/:|p)  =  (0|£^|p).  (3-33) 

This  is  simply  Eq.  (3-5).  The  third  term  in  the  top  equality  of  Eq.  (3-32)  disappears 
because  of  the  solution  of  the  coupled-cluster  equations.  The  last  term  in  top  equality 
of  Eq.  (3-32)  disappears  because  Qji  acting  on  |q)  gives  back  |q),  and  £  cannot  give 
non- vanishing  terms  between  (0|  and  |q). 

For  the  Zeta  equations  a  pure  EOM-CC  method  will  again  be  assumed.  Taking  the 

derivative  of  the  functional  with  respect  to  T  gives 

dF  - 

gj;  =  0  =  {0\C[lH,nT],n]\0)  +  {0\[H,nT]\0)  +  {0\Z[HMT]\0) 

=  {0|£[^.f}^]7^|0)-(0|£7^[^.f]^]|0)  +  (0|[^,O^]|0)  (3-34) 

+  {o\z[H,nT]\o). 
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Notice  that  the  first  and  last  terms  are  the  same  as  Eq.  (3-18).  Expanding  the  commutators 

and  inserting  resolutions  of  the  identity  gives 

0  ={0\CH\q){q\ftTn\0)  +  (0|£F|p)(p|Q^7^|0)  -  (0|£f^T|q)(q|^7^|0) 

-  {o\cnT\p){p\Hn\Q)  -  (o|£7e|q){q|[^,nr]|o) 

-  (0|£7e|g)(g|  [H,  ^t]  |0)  -  (0|£7e|0)(0|  [H.  ^t\  |0)  +  (0|  [iJ.  Qr]  |0)       (3-35) 
+  (0|Z|q)(q|^(^r|0)  +  (0|Z|g)(g|^Or|0)  +  (0|Z|0)(0|fl^f]r|0) 

-  (0|Z.Qr|q)(q|^|0)  -  (0|Znr|g)(g|^|0)  -  (0|Znr|0){0|F|0). 

The  third,  fifth,  and  twelfth  term  disappear  for  the  same  reason  that  the  (0|£fi7;;|q) 

disappeared  above.  Since  the  norm  is  chosen  to  be  one,  the  seventh  and  eighth  terms 
cancel.  The  ninth  and  eleventh  terms  disappear  since  Z  cannot  connect  (0|  to  either  |0) 
or  |q).  The  thirteenth  term  disappears  because  of  the  solution  of  the  coupled-cluster 
equations.   Eqs.  (3-5)  and  (3-6)  can  be  used  to  simplify  the  second  and  fourth  terms. 

Doing  this  gives 

0  =(0|£F|q)(q|7l|g)  +  £;(0|£|p)(p|f2^7^|0)  -  £;(()|£fi^|p>(p|7^|0) 

-  (0|£7^|g)(g|^^^|0)  +  (0|£7e|g)(g|nr|q)(q|^|0) 

(3-36) 
+  (0|£7e|g)(g|f^r|g)(g|^|0)  +  {0|£7l|g)(g|^r|0)(0|^|0) 

+  (o|z|g)(g|^|g)-(o|Z|g)£;cc. 

The  fact  that  R  and  Vij,  both  being  excitation  operators,  commute  has  been  used.  Since 
R  acting  on  |0)  only  connects  with  (p|  and  since  C  acting  on  (0|  only  connects  with  |p), 
the  second  and  third  terms  cancel.  The  fifth  term  disappears  since  fij  cannot  connect 
(g|  to  |q).  The  sixth  term  disappears  because  of  the  solution  of  the  coupled-cluster 
equations.   Finally,  the  Zeta  equation  becomes 

(0|Z|g)(g|(£;ccl  -  ^)|g)  =(0|£^|q)(q|7e|g) 

(3-37) 
-^(0|£7^|g)(g|(i?ccl-^)|g), 
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or 

(0|Z|g)  =  (0|£^|q)(q|7^|g)((g|(£;ccl  -  ^)|g))"'  +  {0\Cn\g).  (3-38) 

The  first  term  in  Eq.  (3-38)  is  the  same  as  Eq.  (3-19).  However,  the  second  term  is  new. 
The  expression  for  the  energy  derivative  corresponding  to  this  functional  is 

—  =  {0\C[m,n]  |0)  +  {0\H^\0)  +  (0|Z^^|0).  (3-39) 

In  fact,  this  equation  and  Eq.  (3-20)  produce  the  same  density  matrices.  In  the  first 
term  of  Eq.  (3-20)  H^  and  7?.  are  not  required  to  be  connected,  as  they  are  in  Eq. 
(3-39).  Those  disconnected  term  can  be  divided  up  into  two  categories.  In  the  first,  C 
is  connected  to  only  7^.  When  all  such  terms  are  added  up,  the  norm  can  be  factored 
out,  leaving  the  second  term  of  Eq.  (3-39).  The  other  category  is  when  C  is  connected 
to  H^.  These  terms  lead  to  the  extra  term  in  Eq.  (3-38). 

It  is  important  to  note  that  the  two  different  functional  have  different  equations  for 
Z,  but  since  they  give  the  same  density  matrices,  either  may  be  used.  Therefore,  the 
operator  Z  is  not  unique,  but  only  has  meaning  within  the  context  of  its  functional. 

To  better  understand  the  differences  between  the  functionals,  let  us  examine  how  the 
ground  state  fits  into  each.  In  EOM-CC  theory,  the  ground  state  is  just  an  eigenstate  of 
H  with  the  special  property  that  7^  =  1,  and  £  =  1  +  A,  where  A,  a  pure  deexcitation 
operator,  is  the  usual  A  from  coupled-cluster  theory. "^"^  Substituting  these  values  into 
the  first  functional  gives 

F=(0|(l  +  A)fl^|0),  (3^0) 
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the  usual  ground  state  functional.  In  this  functional  Z  for  the  ground  state  is  zero.  This 
can  be  seen  from  Eq.  (3-19),  where  the  constant  operator  1  cannot  connect  (q|  and  |g). 
These  give  for  the  energy  derivative,  Eq.  (3-20), 

—  =  {()\{1+A)HX\0),  -■'''/'■  (3^1) 

again  the  usual  coupled-cluster  expression. 

For  the  second  functional  the  picture  is  somewhat  more  complicated.  Substituting  the 
appropriate  values  for  TZ  and  C  into  Eq.  (3-38)  gives  that  (0|Z|g)  =  (0|(1  +  A)l|g)  = 
(0|A|g).  In  other  words,  Z  for  the  ground  state  in  the  second  functional  is  A.  But  in 
this  case,  72.  for  the  ground  state,  a  constant  operator,  gives  zero  for  the  commutators  in 
the  first  terms  of  Eqs.  (3-23)  and  (3-39).  So  the  usual  coupled-cluster  expressions,  Eqs. 
(3-40)  and  (3^1),  are  still  obtained.  .       ,'  "       "       '    . 

The  above  equations  derived  from  the  second  functional  have  a  small  computational 
advantage  over  those  derived  from  the  first  functional.  Specifically,  the  second  Zeta 
equation  suggests  a  simpler  way  to  calculate  some  of  the  terms  in  the  density  matrix 
equations.  Stanton  and  Gauss^^  suggest  combining  Z  and  tqC  into  a  composite  operator, 
and  using  this  new  operator  in  a  ground  state  coupled-cluster  gradient  code  to  calculate 
many  of  the  excited  state  density  matrix  contributions.  The  approach  presented  here 
suggests  including  all  of  (0|£7?.|g)  into  the  composite  operator,  reducing  even  further  the 
number  of  terms  which  would  need  to  be  calculated  separately.  In  practice,  this  should 
make  very  little  difference  in  terms  of  how  long  a  calculation  takes.  However,  it  should 
make  programming  somewhat  easier. 
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It  is  possible  to  derive  a  third  Z  equation.  Starting  with  the  first  line  of  Eq.  (3-34) 
and  not  breaking  the  commutators  gives 

{0\Z\g)  =  {{0\C[[H,n],nT]\0)  +  {0\[HMT]\0)){{g\{Eccl  -  H)\g))-\      (3-42) 

In  the  first  term  H  must  be  connected  to  both  7^  and  flj,  while  in  the  second  term  H 
is  connected  to  just  fir,  but  TZ  and  C  cannot  appear.  In  diagrammatic  language  saying 
that  H  is  connected  to  ^j  is  equivalent  to  saying  that  at  least  one  of  the  lines  from  H 
must  be  dangling  from  the  bottom. 

Since  the  second  and  third  Z  equations  were  derived  starting  from  the  same  place, 
they  must  be  equivalent.  For  EOM-CCSD  the  third  Z  equation  has  two  more  diagrams 
than  the  first  Z  equation.  None  of  the  diagrams  where  the  two  differ  would  change  the 
overall  computational  scaling  of  the  method. 

P-E0M-MBPT(2)  Gradients 

The  defining  equations  for  P-E0M-MBPT(2)  are  identical  to  Eqs.  (3-4)  and  (3-5),  but 
with  a  modified  H.  Detailed  equations  for  H,  TZ,  and  C  have  been  given  previously.""^ 
For  (non-partitioned)  E0M-MBPT(2)  the  functional  is^'' 

F=(o|£[(/^[o]  +  //W+//WTW)7^]    |o) 
+  (o|//toi  +  /f[nf//WTW)    |o\ 
+  (o|z[(//toi-^o)TW  +  //W]  |o) 

The  c  means  that  the  operators  are  connected;  they  share  at  least  one  index.  Here,  the 
usual  many-body  perturbation  theory  partitioning  of  the  Hamiltonian  is  assumed,   /f '"1 


42 


is  the  nth  order  term  of  the  Hamiltionian,  and  Tt"'  is  the  nth  order  T  amplitude  from 
many-body  perturbation  theory.  For  P-E0M-MBPT(2)  the  functional  becomes 

F  =(o  I  (A  +  £2) [(//'°i  +  //'^i  +  //'^'TW)7ei]  1 0) 


+    0 


+ 


(o|£i[(/f[oi  +  //W  +  //WrW)7i2]  |o) 


+    0 


(3^4) 


//[01+//[i]  +  (^^[i]t[i]J    |o 
+  (o|z[(7/[oi-£;M)TW  +  /fW]  |o) 

+  e[i-U)  en  o\y 

This  form  of  the  Zeta  operator  comes  from  the  defining  equation  for  the  T'^'  amplitudes, 
which  is 


g 


^/^[0]_£;[0lJ7^[ll  +  ^[l] 


0)  =0. 


(3^5) 


For  the  rest  of  the  paper  it  will  be  assumed  that  the  reference  determinant  is  the 
Hartree-Fock  solution.    In  this  case. 


H^ 


EM'-^4  +  E^4'^'^'^]- 


(3-*6) 


and 


(3^7) 


pqr.s 


Also, 


TW  =  l^,fl{aW,}, 


(3-48) 


'•J 


where 


'■^  £i  +  £j  -£a-£b' 


(3^9) 


0  )  must  vanish, 
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Since  T'^J  is  now  a  double  excitation  operator,  /o    £i  H^^^T^^^TZ2 
and  the  functional  is  slightly  simpler.  Another  consequence  of  T^^^  being  a  pure  double 
excitation  is  that  Z  consists  only  of  double  deexcitation  operators. 

Taking  derivatives  of  the  (simplified)  Eq.  (3^W)  gives  the  equations  for  the  gradients. 
It  can  be  shown  that  taking  derivatives  with  respect  to  the  various  7^  and  C  operators 
recovers  the  energy  expression.  Again,  by  construction,  taking  the  derivative  with  respect 
to  Z  recovers  the  T'^I  equation.  Finally,  taking  the  derivative  with  respect  to  Tt^l  gives 

(3-50) 


(o|z(/f[oi-E[oi) 


+ 


0>  +  ^0 

d\. 


Here,  fl^  =  ^   Yl   {a^ib^j}  is  the  pure  doubles  excitation  operator.   The  fact  that  TZ 

i,j,a,b 

commutes  with  Qd  has  been  used.   Solving  for  Z  gives 


(o  I  z  I  d)  ={(o  I  £(//W7eind)^  h)  +  (H  ^'''  I  ^)} 

x((d|(^[0]l-i/[0])|d))"'. 


(3-51) 


Previously,  solving  for  Z  involved  inverting  the  full  H,  but  now  only  the  inverse 
of  ^d  ^^[0]  -  fM  j  I  d^  is  needed,  which,  in  the  Hartree-Fock  case,  is  just  the 
denominator  from  MBPT(2).  It  is  now  convenient  to  introduce  a  new  operator,  S,  which 
is  the  inhomogeneous  part  of  the  Z  equation.  For  P-E0M-MBPT(2) 

a,b 

(3-52) 


a,b 


=  (o  I  £(i7W7^lf^d)  J  o)  +  (o  I  jfW  I  d). 

Computationally,  to  solve  for  Z,  H  must  first  be  calculated  and  then  the  denominator 
applied. 
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Now  the  derivative  can  be  rewritten  as 

dE      OF 


dx       dx 


=  (0 


.(0 


\   5x  5x  dx         J     ^ 


°) 


+ 


(o  I  £2(//I"1R2) 
.(0 

Wo 


0 


ai/[oi     5//tii     a//tii  Ml 


(3-53) 


dx  dx 


dx 


0 


z|£^tW  +  ^*'" 


0 


5x   ^       '     5x 

This  expression  only  contains  quantities  which  are  independent  of  the  perturbation  plus 
derivatives  of  the  Hamiltonian  elements.  Thus,  the  derivative  can  be  recast  as  one-  and 
two-particle  effective  density  matrices  times  derivatives  of  integrals; 


dx       ^^  dx       ^''      dx 


(3-54) 


Explicit  spin  orbital  equations  for  the  effective  density  matrix  elements  in  terms  of  Tl, 
C,  Til],  and  Z  are 


e  m,e.J 

n?  —  \^  /"V^    -I-         \^    Amn.ae 
Pb  -  Z^h  ^m  +  2    2^   ^^'^      "1"' 
6  Tn,n,a 


m    ae 

e     im  1 


m,n,e,f 


m,n,e,f 


m,e 


pI  =  0, 


Pki  =  ^P-{iJ)P-{kl)Sa.6ji, 


(3-55) 

(3-56) 

(3-57) 
(3-58) 
(3-59) 
(3-60) 


-Si^-i"- 
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m.e 


Pcd  =  0, 


(3-61) 


^  ni.ej  m,ej 

m,e,f 

pij--^P-iiJ)P-iah)lir^,  (3-62) 


"^  m,n,e  m,n,e 

1  ^  .  (3-63) 

+  -P.{ab)  Y^  Cc^te, 

Tn,n,e 

Phc-^l^kc^m^  (3-64) 

rj). 

"i,e  m,e 


(3-65) 


(3-66) 


ah  n. 

Pij  -  u-  (3_67) 

The  explicit  spin  orbital  equation  for  E  is 

"ijfi  m,e 

+  P-(ii)P_(a6)5^Z>^(mj||e6) 
-P-(v)P_(a6)^Cr^(^j||nfe) 

m,n,e 

+  P-{iJ)P-{<^h)Yil<,{fj\\eh) 

m,e,/ 
m,e,/ 

-P.{rj)  Y  Qr^^iJmWne)  -  ^  J'iym{fm\\ba) 

m,n,e  „,  ,,j 

-  YCb'rtnmen)  +  {tj\\ah). 

m,n,e 


(3-68) 
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In  these  equations  P-{pq)  stands  for  [1  -  P{pq)],  where  P{pq)  means  permute  p  and  q. 

The  equations  for  p  and  for  (  are  subsets  of  the  corresponding  equations  for  EOM- 
MBPT(2)  gradients,^''  and  the  P-E0M-MBPT(2)  gradients  have  been  implemented  in  the 
ACES  II  program  system*  starting  from  Stanton  and  Gauss's  E0M-MBPT(2)  gradient 
code.^^ 

One  final  note  should  be  made  about  the  cost  of  the  P-E0M-MBPT(2)  gradients. 
Previously  it  was  argued  that  the  iterative  Ti^^ocN^irt  steps  should  dominate  the  time  over 
the  non-iterative  nl^^N^-^^  step  for  normal  vertical  excitation  energy  calculations.  This 
was  partially  because  the  cost  of  the  n^^^N^-^f  step  can  be  amortized  over  multiple 
excited  states.  However,  the  gradient  for  only  one  excited  state  can  be  calculated  at  a 
time.  Therefore,  this  amortization  cannot  occur.  More  importantly,  there  are  now  many 
"^ore  nl^^N^.^^  steps,  which,  when  added  together  will  also  contribute  significantly  to 
the  cost  of  the  calculation.  Therefore,  P-E0M-MBPT(2)  gradients  is  most  accurately 
considered  a  non-iterative  nl^^N^^^^  method. 

Application  to  Diatomic  Molecules 

In  order  to  assess  the  performance  of  the  CIS(D)83  method  for  excited  state  potential 
energy  surfaces,  Stanton  et  alP  studied  singlet  excited  states  for  six  diatomic  molecules. 


The  ACES  II  program  is  a  product  of  the  Quantum  Theory  Project,  University  of 
Florida.  Authors:  J.  F.  Stanton,  J.  Gauss,  J.  D.  Watts,  M.  Nooijen,  N.  Oliphant,  S.  A. 
Perera,  P.  G.  Szalay,  W.  J.  Lauderdale,  S.  R.  Gwaltney,  S.  Beck,  A.  Balkova,'  D.  E 
Bernholdt,  K.-K.  Baeck,  H.  Sekino,  P.  Rozyczko,  C.  Huber,  and  R.  J.  Bartlett.  Integrals 
packages  included  are  VMOL  (J.  Almlof  and  P.  R.  Taylor),  VPROPS  (P.  R.  Taylor),  and 
a  modified  version  of  the  ABACUS  integral  derivative  package  (T.  U.  Helgaker,  H.  J. 
Aa.  Jensen,  J.  Olsen,  P.  J0rgensen,  and  P.  R.  Taylor). 
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In  order  to  compare  this  new  work  with  the  methods  presented  there,  those  systems  have 
also  been  studied  here.  The  six  diatomics  are  CO,  N2,  C2,  H2,  BH,  and  BF.  Except  for 
N2,  the  studied  states  are  the  lowest  singlet  excited  state  of  the  molecule.   For  N2  the 
I'llu  state  was  studied. 

The  basis  sets  used  were  6-31G*,'^^  aug-cc-pVDZ,'-^'  and  aug-cc-pVTZ.'^'^  From 
attempting  to  match  the  reported  numbers,  it  would  appear  that  for  the  6-3 IG*  and 
aug-cc-pVDZ  basis  sets  all  six  Cartesian  d  functions  were  used  previously,  and  for  the 
aug-cc-pVTZ  basis  set  only  the  spherical  d  and  /  functions  were  used.  The  same  was 
done  here. 

These  states  have  been  carefully  studied  with  very  accurate  calculations.  However, 
the  point  of  this  work  is  to  compare  the  current  method  with  other  single  reference 
excited  state  methods.  Therefore,  an  attempt  to  review  the  volumous  literature  on  these 
molecules  will  not  be  made.  Instead  the  three  inexpensive  methods  (CIS,  CIS(D),  and  P- 
E0M-MBPT(2))  will  be  compared  to  the  more  expensive  and  more  complete  EOM-CCSD 
and  to  experiment. '-^^  All  of  the  other  methods  can  be  viewed  as  various  approximations 
to  EOM-CCSD.  For  valence  excited  states  there  are  frequently  important  contributions 
from  triples  to  consider.^^^-' 


^  Basis  sets  were  obtained  from  the  Extensible  Computational  Chemistry  Environment 
Basis  Set  Database,  Version  1 .0,  as  developed  and  distributed  by  the  Molecular  Sciences 
Laboratory  which  is  part  of  the  Pacific  Northwest  Laboratory,  P.  O.  Box  999,  Richland, 
Washington  99352,  USA,  and  funded  by  the  U.S.  Department  of  Energy.  The  Pacific 
Northwest  Laboratory  is  a  multi-program  laboratory  operated  by  Battelle  Memorial 
Institute  for  the  U.S.  Department  of  Energy  under  contract  DE-AC06-76RLO  1830. 
Contact  David  Feller,  Karen  Schuchardt,  or  Don  Jones  for  further  information. 
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Results  for  the  diatomic  molecules  are  presented  in  Tables  3-1  to  3-5.   There  are 
some  small  differences  between  the  numbers  reported  here  and  the  values  given  in  Ref. 
93.  Several  values  given  in  the  previous  paper  were  in  error.  Those  errors  have  been 
corrected  here. 

Vertical  and  Adiabatic  Excitation  Energies 

The  vertical  excitation  energies  are  given  for  the  six  molecules,  for  the  four  methods, 
and  for  the  three  basis  sets  in  Table  3-1.  The  adiabatic  excitation  energies  are  given  in 
Table  3-2.  The  ground  state  minimum  geometries  used  are  listed  in  Table  3-3,  with  the 
excited  state  minimum  geometries  in  Table  3-4.  The  appropriate  ground  state  minimum 
for  CIS  is  the  SCF  minimum.  The  appropriate  ground  state  minimum  for  CIS(D)  and  P- 
E0M-MBPT(2)  is  the  MBPT(2)  minimum.  The  appropriate  minimum  for  EOM-CCSD 
is  the  CCSD  minimum. 

As  has  been  noted  before,^-'  the  CIS  answers  are  too  poor  to  be  reliable  CIS  has 
several  errors  of  greater  than  1  eV,  and  for  C2  CIS  even  predicts  the  wrong  sign  for 
the  excitation  energy.  The  mean  absolute  deviation  with  respect  to  experiment  for  the 
adiabatic  excitation  energies  with  the  aug-cc-pVTZ  basis  is  0.724  eV.  These  failings  are 
directly  attributable  to  the  method's  complete  lack  of  dynamic  correlation. 

CIS(D)  improves  on  the  CIS  in  all  cases,  except  those  four  where  the  CIS  is 
accidentally  already  in  good  agreement.  In  all  cases,  the  CIS(D)  energies  are  within 
0.2  eV  of  the  EOM-CCSD  energies,  and  in  all  of  the  cases  with  the  best  basis  set,  the 
error  with  respect  to  experiment  is  less  than  0.3  eV.  The  mean  absolute  deviation  with 
respect  to  experiment  for  the  adiabatic  excitation  energies  with  the  aug-cc-pVTZ  basis  is 
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0.1 12  eV.  This  is  partially  due  to  the  choice  of  excited  states  considered  and  is  actually 
better  accuracy  than  was  initially  reported  for  the  method.^^ 

The  P-E0M-MBPT(2)  method,  on  the  other  hand,  actually  performs  worse  than 
originally  reported.  ""^  The  vertical  excitation  energies  are  slighdy  worse  than  the  CIS(D) 
energies.  The  adiabatic  energies  are  slighUy  better  to  somewhat  worse  than  the  CIS(D) 
energies,  with  their  mean  absolute  deviation  with  respect  to  experiment  for  the  aug-cc- 
pVTZ  basis  being  0.279  eV. 

The  fact  that  the  P-E0M-MBPT(2)  energies  are  worse  than  the  CIS(D)  energies  seems 
to  be  caused  by  the  choice  of  states  and  by  a  curious  feature  of  CIS(D).  In  the  previous 
smdy,"4  it  was  noted  that  P-EOM-MBPT(2)  performed  better  on  Rydberg  states  than 
on  valence  states.  This  can  be  understood,  since  Rydberg  excitations  essentially  involve 
pulling  an  electron  out  of  the  valence  region  and  putting  it  in  a  very  diffuse  orbital,  while 
valence  excitations  put  the  electron  back  into  the  valence  space  in  a  different  arrange- 
ment. Thus  valence  excitations  should  involve  more  orbital  relaxation  and  differential 
correlation  compared  to  Rydberg  excitations.  Therefore,  it  is  reasonable  that  a  method 
as  simple  as  P-E0M-MBPT(2)  could  have  more  trouble  accurately  describing  valence 
excitations  than  Rydberg  excitations.  EOM-CCSD  is  also  more  accurate  for  Rydberg 
states  than  valence  states,  where  it  takes  triples  to  partially  correct  the  problem. '^3 

CIS(D),  on  the  other  hand,  performs  better  for  valence  states  than  for  Rydberg 
states. '14  This  is  curious,  since  the  method  does  not  allow  the  states  to  relax  in  the 
presence  of  electron  correlation.  All  of  the  states  studied  here  are  valence  states,  which 
is  probably  why  the  CIS(D)  energies  are  better  than  the  P-EOM-MBPT(2)  energies. 
Also,  these  states  are  well  separated  from  other  states  of  the  same  symmetry,  meaning 
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that  mixing  of  states,  which  the  zeroth  order  wavefunction,  the  CIS,  would  have  trouble 
handling,  does  not  occur  in  these  problems. 

The  EOM-CCSD  aug-cc-pVTZ  results  agree  quite  well  with  experiment,  which  is  not 
unexpected  for  such  relatively  simple  states.  The  mean  absolute  deviation  with  respect 
to  experiment  for  the  adiabatic  excitation  energies  is  0.106  eV.  Actually,  there  are  not 
significant  differences  between  the  basis  sets,  except  for  H2,  where  6-3 IG*  does  not 
contain  any  polarization  functions. 

Bond  Lengths  and  Vibrational  Frequencies 

Excited  state  equilibrium  bond  lengths  are  presented  in  Table  3^  and  excited  state 
vibrational  frequencies  are  presented  in  Table  3-5.  Certain  trends  hold  true.  CIS 
bond  lengths  are  always  too  short  and  vibrational  frequencies  are  all  too  high.  This 
is  completely  analogous  to  the  situation  for  Hartree-Fock  for  the  ground  state.  However, 
like  Hartree-Fock,  the  bond  lengths  and  vibrational  frequencies  are  reasonable.  The 
CIS(D)  and  P-E0M-MBPT(2)  bond  lengths  and  vibrational  frequencies  are  quite  similar. 
And,  as  should  be  expected,  the  EOM-CCSD  is  superior  to  the  other  methods. 

For  CIS,  the  mean  absolute  deviation  in  the  bond  length  is  0.028  A,  and  the  mean 
absolute  deviation  in  the  frequency  is  188  cm"'.  For  comparison,  the  mean  absolute 
deviation  in  the  ground  state  bond  lengths  for  the  SCF  was  0.015  A.  The  mean  absolute 
deviation  in  the  bond  length  drops  to  0.015  A  for  CIS(D),  with  the  mean  absolute  deviation 
in  the  vibrational  frequencies  being  108  cm"' .  This  compares  to  a  mean  absolute  deviation 
of  0.011  A  for  the  P-EOM-MBPT(2)  bond  lengths,  106  cm"'  for  the  P-E0M-MBPT(2) 
vibrational  frequencies,  and  0.009  A  for  the  MBPT(2)  ground  state  bond  lengths.  The 
EOM-CCSD  bond  lengths  have  a  mean  absolute  error  of  0.008A,  and  the  vibrational 
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frequencies  have  a  mean  absolute  error  of  67  cm"'.  The  ground  state  CCSD  has  a 
mean  absolute  deviation  of  0.004  A  for  the  bond  length.  All  comparisons  were  made 
between  the  aug-cc-pVTZ  results  and  the  experimental  results.  In  general,  the  ground 
state  methods  were  slightly  better  than  the  excited  state  methods,  and  EOM-CCSD  was 
better  than  P-E0M-MBPT(2),  which  was  slightly  better  than  CIS(D). 

Application  to  Polyatomic  Molecules 

The  P-E0M-MBPT(2)  method  has  also  been  applied  to  some  prototypical  polyatomic 
molecules.  Because  of  the  much  more  complicated  chemistry  possible,  polyatomics  can 
provide  a  much  more  demanding  test  of  a  method. 

The  S]  State  of  Ammonia 

The  Si  state  of  ammonia  has  been  studied  extensively. '^'*'  '^^  As  a  methods  test, 
one  of  its  attractive  features  is  its  predissociative  nature.  The  character  of  the  state 
changes  from  Rydberg  like  near  its  Dsh  minimum  to  valence  like  near  the  C2v  transition 
state. '2^  Therefore  for  a  method  to  properly  describe  the  barrier  height,  it  must  equally 
treat  the  valence  and  Rydberg  parts  of  the  potential  energy  surface.  Studying  the  S]  state 
of  ammonia  will  also  provide  an  opportunity  to  compare  the  P-E0M-MBPT(2)  with  the 
(non-partitioned)  E0M-MBPT(2)  method^''  in  order  to  assess  the  effect  of  the  partitioning 
apart  from  the  effect  of  replacing  the  CCSD  ground  state  with  the  MBPT(2)  ground  state. 

The  basis  set  used  for  this  study  was  the  "A"  basis  set  of  Ref.  127.  This  basis  set 
contains  65  contracted  functions  and  includes  multiple  diffuse  functions  at  each  atom.  It 
is  flexible  enough  to  give  a  reasonable  description  of  the  entire  area  of  interest  of  the 
potential  energy  surface.  Table  3-6  gives  the  geometry  and  vibrational  spectrum  for  the 
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D3h  minimum.  It  is  believed  that  the  N-H  bond  distance  is  1.055  ±  0.008  A.'^^  The 
EOM-CCSD  bond  distance  falls  within  these  error  bars.  However,  a  larger  basis  set  would 
tend  to  shorten  it.'^^  The  CIS  bond  length  is  once  again  too  short  and  the  frequencies 
are  too  large.  The  CIS(D)  ,  P-E0M-MBPT(2),  and  E0M-MBPT(2)  bond  lengths  are 
all  similar  and  slightly  shorter  than  the  EOM-CCSD  bond  length.  For  the  vibrational 
frequencies,  the  CIS(D)  results  are  slightly  more  erratic  compared  to  the  EOM-CCSD 
than  the  are  P-E0M-MBPT(2)  and  E0M-MBPT(2),  which  are  similar. 

The  infrared  intensities  show  some  interesting  patterns.  In  every  case  CIS(D) 
substantially  overestimates  the  correlation  correction  to  the  intensity.  The  P-EOM- 
MBPT(2)  is  only  partially  able  to  correct  that  behavior.  Most  of  the  error  in  the  P-EOM- 
MBPT(2)  can  be  traced  to  the  very  approximate  treatment  of  the  doubles-doubles  block, 
since  the  EOM-MBPT(2)  intensities  are  much  closer  to  the  EOM-CCSD  intensities.  The 
intensities  provide  a  very  sensitive  measure  of  the  quality  of  the  wavefunction,  and  these 
large  errors  for  the  P-E0M-MBPT(2)  and  especially  the  CIS(D)  methods  suggest  that 
their  descriptions  of  the  wavefunction  are  not  as  good  as  the  bond  lengths  and  vibrational 
frequencies  suggest. 

The  structure  and  vibrational  spectrum  of  the  transition  state  are  presented  in  Table 
3-7.  The  CIS(D)  bond  length  for  the  hydrogen  being  extracted  is  much  too  long.  This 
bond  length  for  the  P-E0M-MBPT(2)  method  is  also  too  long,  although  shorter  than 
for  CIS(D).  CIS  actually  manages  to  get  this  bond  length  right  compared  to  EOM- 
CCSD.  These  elongated  bonds  for  CIS(D)  and  P-E0M-MBPT(2)  cause  uj(,  to  be  too 
small^-^  with  the  CIS(D)  actually  predicting  that  the  transition  state  is  not  Ciy.  Stanton 
et  alP   suggested  that  the  long  bond  length  was  caused  by  the  tendency  of  CIS(D)  to 
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underestimate  the  energy  of  Rydberg  states,^^  thus  causing  a  "late"  transition  state.  For 
P-E0M-MBPT(2)  the  problem  is  the  opposite.  It  tends  to  overestimate  the  valence  states. 
However,  that  will  also  lead  to  a  "late"  transition  state  and  a  longer  C-H  bond. 

Once  again,  CIS(D)  and  P-E0M-MBPT(2)  have  problems  with  the  intensities,  with  I2, 
I3,  and  I5  having  significant  errors  with  respect  to  EOM-CCSD  and  E0M-MBPT(2).  On 
the  other  hand,  the  CIS(D)  and  P-E0M-MBPT(2)  dipole  moments  are  pretty  reasonable. 
The  barrier  heights,  though,  are  much  too  large,  with  the  CIS(D)  barrier  height  being 
much  worse  than  the  P-E0M-MBPT(2)  barrier  height.  Even  the  EOM-CCSD  barrier 
height  is  well  above  the  experimental  barrier  height,  estimated  to  be  about  2100  cm~'.^28 
Tran^-bent  and  Vinylidenic  Isomers  on  the  S]  Surface  of  C2H2 

The  S]  surface  of  acetylene  provides  a  nice  test  for  excited  state  methods.  The 
C2h  trans-bent  minimum  for  the  Si  state  was  described  in  the  1950's.'^^'32  However, 
in  a  recent  theoretical  paper,  it  was  suggested  that  the  global  minimum  was  a  C2v  vinyli- 
denic structure  with  both  hydrogens  bonded  to  one  of  the  carbons.'^^  The  vinylidenic 
structure  was  calculated  to  be  13  kcal/mol  lower  than  the  trans-hent  acetylenic  structure. 
This  prediction  that  the  vinylidenic  state  was  lower  was  confirmed  with  multireference 
CI,  approximate  coupled-pair  functional, '^4  ^^^  averaged  quadratic  coupled-cluster'^^ 
calculations.  On  the  other  hand,  both  CIS  and  CIS(D)  predict  the  wrong  ordering  of 
the  isomers.^^ 

The  structures,  vibrational  frequencies,  and  energies  of  the  two  isomers  are  presented 
in  Tables  3-8  and  3-9.  The  calculations  were  performed  with  a  TZ2P  basis  set.'^^  Not 
including  zero  point  energy  corrections,  CIS  predicts  the  acetylenic  structure  to  be  more 
stable  by  2.1  kcal/mol.   The  CIS(D)  predicts  the  acetylenic  structure  to  be  more  stable 
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by  0.9  kcal/mol.  The  P-E0M-MBPT(2)  also  has  the  wrong  order,  with  the  acetylenic 
structure  being  2.9  kcal/mol  lower  in  energy  than  the  vinylidenic  structure. 

Simple  Carbonyls 

As  the  simplest  of  the  carbonyls,  formaldehyde's  spectrum  has  been  studied 
extensively.^-'^'  '^^  The  first  excited  state,  n^-K*,  has  two  distinctive  geometrical  fea- 
tures. The  C-0  bond  lengthens  and  the  molecule  becomes  pyramidal.  CIS  severely 
underestimates  the  bond  lengthening, '°^  while  CIS(D)  severely  overestimates  the  bond 
lengthening.^^  Also,  CIS(D)  predicts  the  molecule  to  be  almost  flat.^^ 

Table  3-10  presents  the  P-EOM-MBPT(2)  results.  Other  than  predicting  a  C-0  bond 
that  is  0.02  A  too  long,  the  P-EOM-MBPT(2)  geometry  agrees  well  with  experiment, 
especially  considering  the  relatively  poor  basis  set  used  (6-3 IG*,'^^  the  same  that  was 
used  in  Ref.  93).  The  frequencies  are  also  reasonable,  except  for  uj^.  The  dipole  is 
somewhat  too  large,  but  is  still  in  better  agreement  with  experiment  than  the  CIS(D) 
dipole. 

Table  3-1 1  presents  a  comparison  of  the  formaldehyde,  acetaldehyde,  and  acetone 
geometries,  all  with  the  6-3 IG*  basis  set.'^°  Hydrogen  "a"  is  a  hydrogen  attached  to 
the  carbonyl  carbon.  Hydrogen  "b"  is  the  "in  plane"  hydrogen.  Hydrogen  "c"  is  the 
hydrogen  pointing  into  the  pyramid.  Finally,  hydrogen  "d"  is  the  hydrogen  pointing 
away  from  the  pyramid.  Since  acetaldehyde  and  acetone  are  pyramidal  in  their  first 
excited  state  like  formaldehyde  is,  there  is  no  true  in  plane  hydrogen.  However,  in  both 
cases  there  is  a  hydrogen  with  a  dihedral  angle  only  about  five  degrees  out  of  planarity. 
The  difference  with  the  ground  state,  however,  is  that  the  in  plane  Jiydrogen  points  away 
from  the  oxygen  instead  of  being  hydrogen  bound  to  it.'^^  This  can  be  attributed  to  the 
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removal  of  one  electron  from  the  in  plane  lone  pair  on  oxygen,  letting  the  steric  effects 
dominate  over  the  weakened  hydrogen  bonding. 

Zuckermann  et  aO'^^  carefully  measured  the  spectrum  of  acetone  around  the  origin 
of  the  n^-K*  band.  Their  analysis  led  to  the  positive  analysis  of  three  vibrational 
bands,  with  uncertainty  about  two  others.  Their  assignments,  along  with  the  calculated 
vibrational  frequencies,  are  presented  in  Table  3-12.  Zuckermann  et  al.  could  not 
definitely  distinguish  between  the  first  possibility  and  the  second  possibility;  though 
they  preferred  the  first.  Determining  which  is  correct  depends  upon  the  frequency  of 
uig.  The  computed  frequency  of  368.4  cm"'  is  well  below  the  465.4  cm"'  for  the  second 
possibility.  On  the  other  hand  the  P-E0M-MBPT(2)  is  low  on  all  of  the  other  frequencies 
except  for  i^i2,  considering  that  the  calculated  frequencies  are  harmonic  frequencies  as 
opposed  to  the  measured  fundamental  frequencies.  Also,  in  acetaldehyde  the  P-EOM- 
MBPT(2)  gets  very  good  agreement  between  the  calculated  vibrational  frequency  (368.4 
cm-')  and  the  measured  frequency  (370  cm-')'^^  for  the  corresponding  mode.  The  next 
calculated  mode  in  acetone  is  at  811.1  cm"'.  Hence,  the  calculated  spectrum  supports 
the  second  possibility  over  the  first  possibility.  It  would  take  better  calculations,  for 
example  an  EOM-CCSD  calculation  with  a  better  basis  set  and  possibly  with  triples,  to 
confirm  which  one  it  is. 

Conclusions 

An  alternative  and  simpler  derivation  for  EOM-CCSD  gradients  based  upon  an 
excited  state  EOM-CC  functional  is  presented.  Gradients  for  the  P-E0M-MBPT(2) 
method  have  been  derived  and  their  cost  discussed.     Then,  by  studying  a  series  of 
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diatomics  and  prototypical  polyatomics  the  performance  of  P-E0M-MBPT(2)  for  low- 
lying  valence  potential  energy  surfaces  was  studied  versus  other  single  reference  methods 
for  excited  states  for  which  gradients  have  been  derived. 

The  performance  of  CIS,  CIS(D),  and  EOM-CCSD  has  been  discussed  previously^^, 
and  so  they  will  only  briefly  be  discussed  here.  The  CIS  energies  are  too  poor  to  be 
even  qualitatively  reliable.  However,  the  CIS  geometries  are  normally  quite  reasonable. 
CIS(D)  has  the  opposite  behavior.  Its  energies  for  valence  states  often  has  errors  of  0.3 
eV  or  less,  but  its  geometries  sometime  fail  dramatically.  Even  worse,  there  seems  to 
be  no  indication  of  when  the  CIS(D)  geometries  will  fail.  EOM-CCSD  performed  quite 
well  in  all  of  the  tests  presented  here,  both  for  the  energies  and  for  the  geometries. 

Overall,  the  performance  of  the  P-E0M-MBPT(2)  method  was  mixed.  The  energies 
were  reasonable  but  were  not  as  good  as  had  been  reported  previously.""^  The  geometries 
were  also  reasonable  but  were  noticeably  worse  than  the  EOM-CCSD  geometries.  Over- 
all, in  these  tests  P-E0M-MBPT(2)  performed  very  similarly  to  CIS(D)  when  CIS(D)  did 
not  fail  and  was  qualitatively  correct  for  those  cases  for  which  CIS(D)  did  fail. 

The  discrepancy  between  the  P-E0M-MBPT(2)  and  the  EOM-CCSD  geometries 
is  attributable  primarily  to  the  partitioning  of  the  doubles-doubles  block  in  P-EOM- 
MBPT(2),  since  the  (non-partitioned)  E0M-MBPT(2)  performed  similarly  to  the  EOM- 
CCSD  for  ammonia.  However,  it  is  this  partitioning  that  makes  the  method  attractive, 
since  it  is  the  partitioning  that  reduces  the  cost  from  an  iterative  n^.A^l  .  to  an  iterative 
''^Lc^virf  plus  a  non-iterative  nl^^N^.^^. 
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Table  3-1 :  Vertical  excitation  energies  (in  eV)  for  the  excited  singlet  states 


6-3 IG* 

aug-cc-pVDZ 

aug-cc-pVTZ 

H2 

CIS' 

15.354 

12.639 

12.714 

CIS(D)' 

15.341 

12.574 

12.832 

P-E0M-MBPT(2) 

15.349 

12.574 

12.835 

EOM-CCSD' 

15.257 

12.484 

12.717 

BH 

CIS' 

3.029 

2.849 

2.852 

CIS(D)' 

2.991 

2.864 

2.810 

P-E0M-MBPT(2) 

2.991 

2.875 

2.824 

EOM-CCSD' 

3.117 

2.972 

2.914 

CO 

CIS' 

9.385 

.    9.272 

9.330 

CIS(D)' 

8.911 

8.768 

8.767 

P-E0M-MBPT(2) 

9.196 

9.047 

9.076 

EOM-CCSD' 

8.833 

8.638 

8.666 

N2 

CIS' 

10.378 

10.360 

10.550 

CIS(D)' 

9.444 

9.389 

9.532 

: 

P-E0M-MBPT(2) 

9.883 

9.805 

9.980 

■^ 

EOM-CCSD' 

9.425 

9.339 

9.514 

BF 

CIS' 

6.888 

6.559 

6.601 

i ' 

CIS(D)' 

6.865 

6.481 

6.455 

P-E0M-MBPT(2) 

6.935 

6.560 

6.545     . 

EOM-CCSD' 

6.898 

6.482 

6.454 

C2 

CIS' 

-1.180 

-1.258 

-1.229 

CIS(D)' 

1.310 

1.151 

1.222 

P-E0M-MBPT(2) 

1.423 

1.261 

1.328 

EOM-CCSD' 

1.526 

1.283 

1.308 

'^Ref.  93. 
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Table  3-2:  Adiabatic  excitation  energies  (in  eV)  for  the  excited  singlet  states 


6-3 IG* 

aug-cc- 
pVDZ 

aug-cc- 
pVTZ 

Expt." 

H2 

CIS" 

12.765 

11.269 

11.352     ■ 

CIS(D)'' 

12.980 

11.290 

11.429 

P-E0M-MBPT(2) 

12.943 

11.274 

11.408 

EOM-CCSD'' 

13.110 

11.228 

11.353 

11.3694 

BH 

CIS" 

3.024 

2.845 

2.849 

CIS(D)" 

2.989 

2.861 

2.807 

P-E0M-MBPT(2) 

2.989 

.     2.872 

2.821 

eom-ccsd" 

3.117 

2.972 

2.913 

2.8685 

CX) 

CIS" 

8.784 

8.739 

8.802 

CIS(D)" 

8.269 

8.219 

8.246 

P-E0M-MBPT(2) 

8.566 

8.527 

8.580 

eom-ccsd" 

8.330 

8.229 

8.256 

8.0684 

N2 

CIS" 

9.413 

9.432 

9.582 

CIS(D)" 

8.802 

8.775 

8.864 

P-EOM-MBPT(2) 

9.264 

9.231 

9.346 

eom-ccsd" 

8.760 

8.714 

8.839 

8.5900 

BF 

CIS" 

6.806 

6.519 

6.563 

CIS(D)" 

6.760 

6.428 

6.404 

P-E0M-MBPT(2) 

6.822 

6.507 

6.495 

EOM-CCSD" 

6.811 

6.444 

6.417 

6.3427 

C2 

CIS" 

-1.275 

-1.350 

-1.319 

CIS(D)" 

1.143 

0.997 

1.075 

P-EOM-MBPT(2) 

1.293 

1.138 

1.208 

EOM-CCSD" 

1.301 

1.066 

1.103 

1.0404 

^>Ref 

.   122. 

''^Ref 

.93. 
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Table  3-3:  Equilibrium  distances  (in  A)  for  the  ground  states 

6-31G*  '"S-''-  ^"S-"'^-  Expt« 

pVDZ  pVTZ  P 


Hz 

SCF'' 

,,      0.7300 

0.7481 

0.7345 

MBPT(2)'' 

0.7375 

0.7549 

0.7374 

'      CCSD^' 

0.7462 

0.7617 

0.7431 

0.7414 

BH 

SCF'' 

1.225 

1.233 

1.221 

MBPT(2)'' 

1.233 

1.241 

1.216 

CCSD'' 

1.244 

1.249 

1.220 

1.2324 

CO 

SCF^ 

1.114 

1.110 

1.104 

MBPT(2)'' 

1.150 

1.147 

1.134 

CCSD*' 

1.141 

1.138 

1.124 

1.1283 

N2 

SCF^' 

1.078 

1.078 

1.067 

MBPT(2)*' 

1.130 

1.131 

1.110 

CCSD'' 

1.113 

^      1.113 

1.093 

1.0977 

BF 

SCF^ 

1.260 

1.270 

1.249 

;-  "■      '  ' 

MBPT(2)*' 

1.279 

1.294 

1.264 

CCSD^* 

1.281 

1.296 

1.263 

1 .2626 

C2 

SCF'' 

1.245 

1.253      ■ 

1.241 

MBPT(2)^ 

1.264 

1.276 

1.254 

CCSD^ 

1.252 

1.265 

1.241 

1.2425 

^)Ref. 

122. 

''^Ref. 

93. 

i' 
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Table  3-^:  Equilibrium  distances  (in  A)  for  the  excited  singlet  states 


6-31G* 

aug-cc- 
pVDZ 

aug-cc- 
pVTZ 

Expt.^ 

H2                 CIS'' 

1.544 

1.239 

1.239 

CIS(D)*' 

1.599 

1.256 

1.273    : 

P-EOM-MBPT(2) 

1.626 

1.273 

1.295 

EOM-CCSD*' 

1.616 

1.267 

1.283 

1.2928 

BH                 CIS*' 

1.204 

1.214 

1.204 

CIS(D)'' 

1.219 

1.223 

1.199 

P-E0M-MBPT(2) 

1.216 

1.223 

1.199 

EOM-CCSD*' 

1.241 

1.242 

1.211 

1.2186 

CO                 CIS*' 

1.228 

1.220 

1.213 

CIS(D)'' 

1.295 

1.286 

1.263 

P-E0M-MBPT(2) 

1.293 

1.282 

1.259 

EOM-CCSD*' 

1.252 

1.242 

1.224 

1.2353 

N2                 CIS*' 

1.200 

1.200 

1.192 

CIS(D)*' 

1.250 

1.249 

1.231 

P-E0M-MBPT(2) 

1.245 

1.244 

1.227  '  • 

EOM-CCSD*' 

1.221 

1.220 

1.202 

1.2203 

BF                 CIS*' 

1.316 

1.312 

1.287 

CIS(D)*' 

1.350 

1.349 

1.312 

P-E0M-MBPT(2) 

1.353 

1.349 

1.312 

EOM-CCSD*' 

1.345 

1.342 

1.304 

1.3038 

C2                  CIS*' 

1.293 

1.301 

1.289 

CISCD)" 

1.333 

1.346 

1.320 

P-E0M-MBPT(2) 

1.325 

1.337 

1.312 

EOM-CCSD*' 

1.334 

1.347 

1.318 

1.3184 

^^Ref.  122. 

*')Ref.  93. 

r- 
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Table  3-5:  Harmonic  vibrational  frequencies  (in  cm"')  for  the  excited  singlet  states 


6-3 IG* 

aug-cc- 
pVDZ 

aug-cc- 
pVTZ 

Expt.* 

H2 

CIS" 

1495 

1610 

1589 

CIS(D)" 

1438 

1485 

1422 

P-EOM-MBPT(2) 

1407 

1410 

1337 

EOM-CCSD" 

1428 

1439 

1368 

1358.09 

BH 

CIS'' 

2576 

2536 

2544 

CISCD)" 

2426 

2441 

2517 

P-E0M-MBPT(2) 

2445 

2443 

2515 

EOM-CCSD'' 

2180 

2243 

2372 

2251.0 

CO 

CIS'' 

1646 

1615 

1633 

CIS(D)'' 

1282 

1238 

1323 

P-E0M-MBPT(2) 

1281 

1244 

1330 

EOM-CCSD'' 

1559 

1517 

1592 

1518.2 

N2 

CIS'' 

1939 

1909 

1897 

CIS(D)'' 

1612 

1583 

1615 

P-E0M-MBPT(2) 

1638 

1612 

1635 

EOM-CCSD'' 

1856 

1825 

1854 

1694.21 

BF 

CIS'' 

1339 

1273 

1367 

CISCD)" 

1173 

1097 

1239 

P-E0M-MBPT(2) 

1160 

1095 

1241 

EOM-CCSD'' 

1199 

1130 

1279 

1264.9 

C2 

CIS" 

1830 

1797 

1794 

CISCD)" 

1618 

1577 

1631 

P-E0M-MBPT(2) 

1674 

1635 

1687 

EOM-CCSD'' 

1605 

1563 

1630 

1608.35 

^^Ref.  122. 

*')Ref.  93. 
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Table  3-6:   Geometries  (in  A),  harmonic  vibrational  frequencies  (in  cm"'),  infrared  intensities 
(in  km/mol),  and  energies  (in  Hartrees)  for  the  D3h  equilibrium  geometry  of  the  Si  state  of  NH3 


CIS« 

CIS(D)^ 

P-EOM- 

MBPT(2) 

EOM- 

MBPT(2)'' 

EOM- 
CCSD*' 

fNH 

1.0213 

1.0441 

1.0420 

1.0487 

1.0512 

a;i(ai') 

3180.8 

2814.0 

2971.4 

3180.8 

2993.1 

^^2(32") 

842.4 

736.5 

767.2 

769.6 

741.2 

^3(e') 

3356.7 

3277.2 

3185.5 

3021.0 

2997.5 

a;4(e') 

1517.9 

1378.6 

1370.6 

1331.1 

1335.2 

I2 

0.1 

57.6 

48.3 

14.7 

9.7 

I3 

5803.2 

1538.5 

3620.2 

4622.1 

4447.1 

I4 

19.4 

981.9 

464.4 

293.0 

376.0 

Energy 

-55.968  635 

-56.240  376 

-56.225  152 

-56.234  687 

-56.246  53^ 

^)Ref.  93. 

1' 

''>Ref.  97. 
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Table  3-7:  Geometries  (in  A  and  degrees),  harmonic  vibrational  frequencies  (in 

cm"'),  infrared  intensities  (in  km/mol),  dipole  moments  (in  Debyes),  energies  (in 

Hartrees),  and  barrier  heights  (in  cm"')  for  the  C2v  predissociative  transition  state  of 

the  S|  state  of  NH3.  The  symmetry  unique  hydrogen  is  denoted  by  an  asterisk. 


CIS^ 

CIS(D)^ 

P-EOM- 
MBPT(2) 

EOM- 
MBPT(2)'' 

EOM- 
CCSD'' 

tnh* 

1.3497 

1.4371 

1.4108 

1.3219 

1.3421 

tnh 

1.0106 

1.0421 

1.0370 

1.0423 

1.0441 

^(H*NH) 

123.43 

126.78 

125.82 

124.51 

124.27 

a;i(ai) 

3538.5 

3133.7 

3203.6 

3069.0 

3051.6 

a;2(ai) 

1568.5 

1539.1 

1525.9 

1478.5 

1456.0 

a;3(ai) 

1446.6/ 

1510.3/ 

1543.2/ 

1977.3/ 

1897.3/ 

u;4(bi) 

1025.3 

1276.0 

1191.1 

1042.0 

1023.1 

a;5(b2) 

3737.2 

3313.9 

3384.0 

3274.9 

3255.2 

W6(b2) 

498.8 

104.8/ 

129.9 

396.5 

419.6 

II 

293.7 

640.4 

590.5 

1081.2 

1010.6 

I2 

5.4 

0.3 

0.0 

48.9 

41.0 

I3 

463.7 

5.5 

112.1 

2301.1 

1921.3 

I4 

50.2 

88.4 

68.0 

43.9 

69.9 

I5 

2.0 

6.7 

3.7 

58.1 

56.5 

l6 

73.0 

66.5 

103.4 

158.0 

167.5 

/' 

3.750 

2.951 

2.985 

2.631 

2.727 

Energy 

-55.952  723 

-56.214  426 

-56.205  848 

-56.223  512 

-56.234  405 

Barrier 
Height 

3492 

5695 

4237 

2452 

2663 

^'Ref.  93. 

''^Ref.  97. 

:^'f  ■■■!        ... 
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Table  3-8:  Geometries  (in  A  and  degrees),  harmonic  vibrational 
frequencies  (in  cm"'),  infrared  intensities  (in  km/mol),  and  energies 
(in  Hartrees)  for  the  trans-bent  isomer  of  the  Si  state  of  acetylene 


CIS^ 

CIS(D)^ 

P-EOM- 

MBPT(2) 

EOM-CCSD'' 

rcH 

1.0788 

1 .0889 

1.0890 

1.0907 

rcc 

1.3521 

1.3724 

1.3717 

1.3575 

^(HCC) 

124.55 

122.06 

121.86 

123.64 

a;i(ag) 

3293.5 

3144.7 

3141.7 

3107.7 

^2(ag) 

1545.1 

1423.6 

1436.4 

1471.4 

a;3(ag) 

1162.4 

1108.2 

1114.1 

1106.6 

u;4(au) 

1230.8 

687.2 

907.8 

614.6 

a;5(bu) 

3280.3 

3131.6 

3126.2 

3083.7 

t^eCbu) 

805.7 

882.4 

887.5 

745.8 

I4 

2.8 

111.1 

40.7 

100.1 

I5 

2.5 

3.6 

3.9 

6.9 

h 

342.4 

281.6 

295.6 

405.9 

Energy 

-76.681  20 

-76.969  01 

-76.966  77 

-76.981  76 

^^Ref.  93. 

':<^-  ' 

''^Ref.  133 

v^^it:- 
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Table  3-9:  Geometries  (in  A  and  degrees),  harmonic  vibrational 

frequencies  (in  cm"'),  infrared  intensities  (in  km/mol),  and  energies 

(in  Hartrees)  for  the  Civ  vinylidenic  isomer  of  the  Si  state  of  acetylene 


CIS« 

CIS(D)^ 

P-EOM- 
MBPT(2) 

EOM-CCSD'' 

rcH 

1 .0847 

1.0881 

1.0879 

1.0890 

rcc 

1.3877 

1.4225 

1.4174 

1 .4289 

^(HCC) 

122.94 

122.15 

122.38 

122.28 

wi(ai) 

3173.5 

3086.1 

3089.2 

3060.1 

a;2(ai) 

1618.8 

1464.2 

1476.7 

1486.5 

t^sCai) 

1381.5 

1225.7 

1246.5 

1213.1 

uj  4(h2) 

3235.5 

3172.0 

3171.0 

3138.9 

UJ5{h2) 

1015.5 

894.6 

903.9 

927.1 

'^6(a2) 

977.1 

702.3 

708.6 

763.2 

Ii 

17.1 

0.1 

0.2 

3.2 

h 

1.3 

29.1 

24.5 

14.4 

h 

0.3 

57.0 

41.7 

31.3 

I4 

9.7 

2.2 

2.1 

7.0 

I5 

4.7 

3.9 

3.5 

5.3 

l6 

1.6 

1.1 

0.4 

2.1 

Dipole 
moment 

3.817 

2.926 

3.029 

2.875 

Energy 

-76.677  85 

-76.967  63 

-76.962  08 

-77.003  07 

^)Ref.  93. 

''^Ref.  133 

,•        ;-.... 
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Table  3-10:  Geometries  (in  A  and  degrees),  harmonic  vibrational  frequencies  (in  cm"'),  infrared 

intensities  (in  km/mol),  dipole  moments  (in  Debyes),  and  energies  (in  Hartrees)  for  the  equilibrium 

geometry  of  the  A  'A"  state  of  formaldehyde.  Experimental  frequencies  are  fundamentals. 


CIS^ 

CIS(D)'' 

P-EOM- 
MBPT(2) 

EOM- 
CCSD'' 

Expt.* 

rco 

1.258 

1.384 

1.340 

1.324 

1.321 

rcH 

1.085 

1.088 

1.093 

1.096 

1.097 

^(OCH) 

117.3 

117.8 

117.0 

115.8 

118.0 

r(HOCH) 

148.4 

171.5 

153.0 

145.9 

148.3 

t^i(a') 

3209.5 

3154.1 

3102.3 

3072.9 

2847 

a;2(a') 

1437.9 

1420.5 

1412.7 

1407.9 

1290 

'^sCa') 

1660.8 

895.7 

1187.2 

1248.5 

1173 

a;4(a') 

561.6 

172.6 

564.9 

719.8 

683 

^5(a") 

3299.6 

3295.8 

3223.6 

3189.5 

2968 

ujeia.") 

993.5 

914.0 

932.9 

934.3 

898 

h 

29.6 

2.6 

0.7 

2.8 

h 

52.0 

6.3 

0.8 

3.5 

h 

2.1 

109.6 

66.8 

93.7 

U 

108.2 

38.3 

36.0 

30.8 

I5 

0.0 

1.6 

2.9 

2.8 

l6 

3.4 

3.6 

3.4 

4.0 

Dipole 
moment 

1.507 

2.374 

2.094 

1.776 

1.56(7f> 

Energy 

-113.669  47 

-114.043  84 

-114.028  60 

-114.052  75 

")Ref.  102. 

^■'::-   ', 

''^Ref.  93. 

'^^Ref.  141. 
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Table  3-11:  A  comparison  of  the  geometries  of  the  n^7r*  states  in  formaldehyde, 
acetaldehyde,  and  acetone.  Geometries  are  in  Angstroms  and  degrees. 


Formaldehyde 

Acetaldehyde 

Acetone 

rco 

1.340 

1.359 

1.376 

rcc 

1.494 

1.494 

rcHa 

1.093 

1.095 

rcHb 

1.091 

1.092 

rcHc 

1.099 

1.101 

rcHd 

1.094 

1.094 

^(OCC) 

116.0 

113.4 

^(OCHa) 

117.0 

114.3 

^(CCHb) 

110.2 

110.3 

^(CCHc) 

110.6 

110.2 

^(CCHd) 

110.2 

110.6 

r(HaCOH) 

153.0 

r(HaCOC) 

149.9 

r(CCOC) 

146.9 

r(HbCCO) 

-177.5 

-174.5 

t(HcCCO) 

61.6 

65.0 

T(HdCCO) 

-57.1 

-53.7 

•'  ■*^«f 
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Table  3-1 2:  A  comparison  of  two  possible  assignments  of  the  low  frequency  vibrations  in 
the  Si  state  of  acetone  to  the  calculated  frequencies.  All  frequencies  are  in  cm"'. 


first 

second 

P-EOM- 

possibility^ 

possibility^ 

MBPT(2) 

iy\2  torsion  (antigearing) 

155.5 

155.5 

206.8 

1^24  torsion  (gearing) 

172.5 

172.5 

190.1 

i^s  C-C-C  bend 

373 

373 

353.8 

1/23  C=0  out-of-plane  wagging 

333 

329.7 

//19  C=0  in-plane  wagging 

177.5 

465.4 

368.4 

^)Ref.  140. 


CHAPTER  4 

THE  SIMILARITY  TRANSFORMED  EQUATION- 

OF-MOTION  COUPLED-CLUSTER  METHOD 

Introduction 

The  equation-of-motion  coupled-cluster  (EOM-CC)  method,^^'  ^^  especially  when 
based  on  a  coupled-cluster  singles  and  doubles  (CCSD)^  ground  state,  offers  an  attractive 
and  unified  formalism  to  extend  single  reference  coupled-cluster  methods  to  excitation 
energies  (EE-EOM-CCSD,  sometimes  just  EOM-CCSD),^^-  ^^  ionization  potentials  (IP- 
EOM-CCSD),^^  and  electron  affinities  (EA-EOM-CCSD).^^  The  method  is  conceptually 
single  reference,  in  that  the  entire  calculation  is  built  upon  one  reference  determinant. 
Therefore,  the  only  choices  necessary  in  the  method  are  the  choice  of  basis  set,  the  choice 
of  reference  determinant,  usually  the  Hartree-Fock  solution,  and  the  choice  of  excitation 
level  to  include  in  the  ground  and  excited  state  operators. 

Several  applications^^'  '^^-us,  123  ^f  EE-EOM-CCSD  and  its  twin,  coupled-cluster 
singles  and  doubles  linear  response,^^'  ^'^  have  been  made.  The  method  has  also  been 
extended  to  include  various  effects  of  triple  excitations.^*^^^  But  EOM-CCSD,  and 
especially  EE-EOM-CCSD,  has  some  drawbacks.  The  first  is  the  cost  of  the  method. 
Every  EE-EOM-CCSD  excited  state  is  approximately  as  expensive  as  a  ground  state 
CCSD  calculation.  Approximations  can  be  made  to  reduce  the  cost,'''^  with  varying 
degrees  of  success.  The  second  drawback  is  the  accuracy  of  the  method.  For  ethylene, 
butadiene,  and  cyclopentadiene,  the  average  error  for  the  Rydberg  states  was  0.17  eV, 
but  the  average  error  for  the  valence  states  was  0.26  eV.'^^      ' 
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The  key  to  EOM-CC  theory  is  the  similarity  transformation  of  the  Hamiltonian.  After 
the  ground  state  coupled-cluster  equations  have  been  solved,  the  cluster  amplitudes  are 
used  to  transform  the  Hamiltonian.  A  similarity  transformation  changes  the  eigenvectors 
of  an  operator  without  changing  its  eigenvalues.  Hence,  the  energies  of  the  excited  states, 
electron  attached  states,  and  ionized  states,  which  are  all  eigenvalues  of  the  Hamiltonian 
matrix  in  terms  of  the  appropriate  bases,  have  not  changed.  ■/ 

So  far  nothing  has  been  accomplished.  To  calculate  the  energies  of  these  states,  it  is 
still  necessary  to  diagonalize  a  matrix  that  is,  in  principal,  infinite.  After  a  finite  basis  set 
has  been  introduced,  the  matrix  is  now  finite,  but  its  size  grows  as  the  size  of  the  basis  set 
factorial.^^  Therefore,  approximations  must  be  introduced.  The  typical  approximations 
involve  truncating  the  Hamiltonian  matrix  to  include  only  a  subset  of  possible  excitations. 
In  EE-EOM-CCSD  the  matrix  is  limited  to  only  single  and  double  excitations,  i.e. 
determinants  where  one  or  two  virtual  orbitals  replace  one  or  two  occupied  orbitals.  For 
EA-EOM-CCSD  the  space  is  limited  to  determinants  where  an  electron  has  been  added 
to  a  single  virtual  orbital  and  determinants  where  along  with  the  electron  being  added 
to  a  virtual  orbital,  an  occupied  orbital  is  replaced  with  a  virtual  orbital.  These  are  the 
so  called  \p  and  2p\h  states.  For  IP-EOM-CCSD  the  space  consists  of  the  l/i  and  2h\p 
determinants,  which  are  those  determinants  with  one  electron  removed  from  an  occupied 
orbital  and  those  determinants  with  one  electron  removed  and  one  electron  excited. 

After  the  truncation  of  the  matrix  is  when  the  similarity  transformation  becomes 
important.  The  truncation  will  change  the  eigenvalue  spectrum.  Many  eigenvalues 
will  no  longer  appear,  and  those  that  remain  will  be  shifted  by  an  amount  which  will 
depend  on  how  much  of  the  true  eigenvector  lies  in  the  omitted  space.  What  gives  the 
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similarity  transformation  its  power  is  that  by  choosing  the  transformation  so  that  the 
important  eigenvectors  are  contained,  to  the  greatest  extent  possible,  in  the  space  kept 
after  the  truncation,  the  effect  of  the  truncation  on  the  energies  of  the  states  of  interest 
can  be  minimized.  Note  that  an  similarity  transformation  is  needed,  since  a  unitary 
transformation  leaves  the  eigenvectors  unchanged.  An  equivalent  picture  is  to  choose 
the  transformation  such  that  the  parts  of  the  matrix  that  couple  the  kept  space  with  the 
excluded  space  are  minimized,  since  these  parts  of  the  matrix  determine  the  extent  to 
which  eigenvectors  span  both  spaces. 

For  singly  excited  states,  the  most  important  terms  deleted  in  the  truncation  to  singles 
and  doubles  are  the  triple  excitations.  In  the  untransformed  Hamiltonian  matrix,  the  terms 
that  couple  single  excitations  to  triple  excitations  are  the  two  electron  integrals  {ah\\ij), 
where  a,  b,  c,  and  d  are  unoccupied  in  the  reference  determinant  and  /,  ;,  k,  and  /  are 
occupied  in  the  reference  determinant.   These  are  precisely  the  terms  which  are  set  to 
zero  by  the  similarity  transformation.  The  transformation  does  add  new  coupling  terms, 
but  they  are  of  the  form  of  a  t2  amplitude  times  a  different  two  electron  integral.  Since 
for  a  normal  single  reference  problem,  the  largest  t  amplitudes  are  of  the  order  of  0.1,  the 
new  coupling  terms  are  at  least  an  order  of  magnitude  smaller  than  the  former  coupling 
terms.    Therefore  the  "true"  eigenvector  for  a  singly  excited  state  is  contained  more 
fully  within  the  space  of  single  and  double  excitations,  and  the  eigenvalue  for  the  state 
with  the  truncated  similarity  transformed  Hamiltonian  is  much  closer  to  the  eigenvalue 
of  the  untruncated  Hamiltonian  than  is  the  eigenvalue  for  that  state  with  the  truncated 
untransformed  Hamiltonian. 
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In  the  similarity  transformed  equation-of-motion  coupled-cluster  (STEOM-CC) 
theory^'*"^^  the  blocking  of  the  matrix  by  similarity  transformations  is  carried  one 
step  farther.  Starting  with  the  EOM-CC  transformed  Hamiltonian,  a  second  similar- 
ity transformation^^  is  performed  in  order  to  minimize  the  coupling  between  the  singles 
space  and  the  doubles  space.  The  dominant  terms  that  couple  these  two  blocks  in  the 
EOM-CC  Hamiltonian  are  the  (afe||c?;)  and  {ia\\jk)  two  electron  integrals.  The  most 
important  of  these  coupling  terms  are  set  to  zero  in  the  second  similarity  transformation, 
leaving  as  the  main  coupling  term  between  the  singles  and  doubles  blocks  again  a  ?2 
amplitude  times  a  two  electron  integral.  Once  again,  the  coupling  is  reduced  by  at  least 
an  order  of  magnitude. 

Since  the  second  similarity  transformation  is  done  in  such  a  way  as  to  preserve 
the  reduction  of  the  singles-triples  coupling  from  the  first  similarity  transformation,  the 
eigenvectors  for  states  dominated  by  single  excitations  are  almost  exclusively  composed 
of  only  single  excitations,  meaning  that  the  truncation  needs  to  leave  only  the  space  of 
singles.  Now  the  final  diagonalization  to  get  the  eigenvalues  and  eigenvectors  only  scales 
as  the  fourth  power  of  the  basis  set,  as  opposed  to  the  sixth  power  scaling  for  EOM- 
CCSD.  The  trade-off  is  having  to  perform  the  second  similarity  transformation,  but  the 
second  transformation  only  scales  as  the  fifth  power  of  the  size  of  the  basis  set. 

If  the  final  diagonalization  is  over  the  space  where  an  electron  has  been  removed 
from  an  occupied  orbital  and  placed  in  an  unoccupied  orbital,  then  excited  states  are 
calculated,  and  the  method  is  known  as  excitation  energy  STEOM-CC  (EE-STEOM- 
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the  double  ionization  potential  STEOM-CC  (DIP-STEOM-CC).^^  Finally,  if  the  space  is 
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over  determinants  with  two  electrons  added,  the  method  is  the  double  electron  attachment 
STEOM-CC  (DEA-STEOM-CC).66  ^^^y^  of  these  methods  provides  an  important  tool  for 
describing  interesting  chemistry,  and  although  the  rest  of  the  dissertation  will  focus  on 
the  EE-STEOM-CC  variant,  everything  said  can  be  applied  to  the  other  two  with  simple 
modifications. 

The  STEOM-CC  Method 

In  STEOM-CC  theory,  the  orbital  space  is  divided  into  occupied  and  virtual  orbitals, 
sometimes  referred  to  as  holes  and  particles,  based  on  the  orbital 's  occupancy  in  the 
reference  determinant.  The  orbital  space  is  then  subdivided  into  active  and  inactive 
occupied  orbitals  and  active  and  inactive  virtual  orbitals.  To  distinguish  the  sets,  the 
labels  /, ;,  k,  and  /  will  refer  to  generic  occupied  orbitals,  and  a,  b,  c,  and  d  will  refer 
to  generic  virtual  orbitals.  Active  occupied  orbitals  will  be  m  and  n  while  active  virtuals 
will  be  e  and/.  Explicitly  inactive  orbitals  will  be  indicated  with  a  prime,  such  as  /'  and 
a'.  The  labels/?,  q,  r,  and  s  can  refer  to  any  orbital.  All  such  orbitals  are  spacial  orbitals, 
as  opposed  to  the  previous  chapters  where  spin  orbitals  have  been  used. 

Although  a  STEOM-CC  calculation  could  be  based  on  any  coupled-cluster  or  many- 
body  perturbation  theory  treatment  of  the  ground  state,  so  far  only  CCSD  and  MBPT(2) 
have  been  used.^  This  discussion  will  focus  on  STEOM-CC  with  a  CCSD^  ground 
state,  giving  the  STEOM-CCSD  method.^^  in  coupled-cluster  theory,  the  ground  state 
of  the  molecule  is  represented  as 

l^o)=e^|0),  (4_1) 
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where  |0)  is  a  closed  shell  single  determinant,  often  the  Restricted  Hartree-Fock  (RHF) 
determinant.  For  CCSD  the  operator  T  consists  of  pure  excitation  operators  of  the  form 

T  =  T,+T2  =  Y.t'^{o^^]  +  \Y.'^^{^^ih^j].  (4-2) 

i,a  '■] 

a.b 

The  curly  brackets  mean  that  the  operator  is  normal  ordered  with  respect  to  |0). 
The  untransformed  Hamiltonian,  in  second  quantization,  is 

/f  =  /,o  +  ^/p,{pt5}  +  -^l/,^„|pV5j.  (4_3) 

p,q  P.I 

The  constant  term,  Iiq,  is  the  energy  of  the  reference  determinant.  The  fpq  refer  to  the 
appropriate  parts  of  the  Fock  operator.  The  V^prs  are  non-antisymmetrized  two  electron 
integrals.  . 

After  the  first  transformation  the  Hamiltonian  becomes  , 

=  ho  +  Yl  hpg  jp^g }  +  4  Z]  ^^Pi'" {p^^(l^^^ } 

p.q  p"  ,  ■     ..  .,         (4-4) 

+  -^^  hpqrstu  \p^sqU.r^u\  +  .... 

With  a  CCSD  reference  state,  H  will  contain  up  to  six-body  terms,  but  the  four-  and  higher 
body  terms  do  not  appear  in  the  EOM-CCSD  or  STEOM-CCSD  equations.  Explicit 
equations  for  these  H  elements,  as  used  here,  can  be  found  elsewhere. ''^^  In  terms  of 
H,  the  CCSD  equations  can  be  written  as 

(4-5) 
The  CCSD  energy  is 

EccSD  =  {0\H\0)  =  7^0.  (4-6) 
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Thus,  solving  the  CCSD  equations  is  equivalent  to  setting  the  one-  and  two-body  pure 
excitation  parts  of  H  to  zero. 

In  STEOM-CC,  a  second  many-body  similarity  transformation^^  is  performed  in  order 
to  also  set  to  zero  selected  remaining  terms  which  increase  the  net  excitation  level  by 
one.  But  since,  in  diagrammatic  language,  the  necessary  operators  have  a  line  on  the 
bottom,  they  do  not  commute  and  can  connect  to  each  other.  Thus  instead  of  using 
a  transformation  of  the  form  e^ ,  as  originally  done  by  Stolarczyk  and  Monkhorst,^^  a 
normal  ordered  exponential  operator  [e^]  is  used.  The  normal  ordered  exponential, 
introduced  by  Lindgren,^^  excludes  terms  inside  the  normal  ordering  from  connecting  to 
each  other,  and  this  simplifies  the  equations. 

The  S  operator  used  in  the  transformation  consists  of  two  parts  S  =  S^  +  S~ ,  where 

5+  =  St  +  St  =  Y:  4{a''e}  +  \  Y:  <j{atefct,|  (4-7) 


a',e 


and 


S-  =  5r  +  S^  =  Y1  s^{mh'}  +  i  Y^  s^'[nM^jy  (4-8) 


It  is  the  presence  of  the  active  ^-annihilation  operators  which  cause  the  components  of 
S  to,  in  general,  not  commute.^ 

The  double  similarity  transformed  Hamiltonian  takes  the  form 

-1 


^  (4-9) 


p,q  P.I 

r.s 
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The  similarity  transformation  preserves  the  zeros  for  the  one-  and  two-body  pure  exci- 
tation parts,  so  that^^ 

gai  =  Ki  =  0, 

(4-10) 

9  alii  j  —  '^abij  —  "• 

Also, 

go  =  ho  =  EccsD-  (4-1 1) 

In  principle,  like  the  CCSD  equations  for  T  above,  the  S  equations  could  be  derived 
by  setting  the  appropriate  parts  of  G  to  zero.   These  are^^ 

Qmi'  =  ($i'|G|$m)  =  0, 

(4-12) 

gmbij  =  {^'ji\G\^m)=0, 
g,6,j  =  ($;''iG|r)=0. 

Note  that  the  second  two  equations  set  to  zero  the  primary  coupling  terms  between  singly 
and  doubly  excited  determinants.  In  practice,  the  5  coefficients  are  not  calculated  this 
way.  Doing  so  would  lead  to  systems  of  nonlinear  equations  which  in  some  situations  can 
be  numerically  unstable.^  Instead,  the  equivalence  of  IP-EOM-CC  and  EA-EOM-CC  to 
Pock  space  coupled-cluster  theory ''^^  is  exploited  to  rewrite  S  in  terms  of  IP-EOM-CCSD 
and  EA-EOM-CCSD  eigenvectors.  These  equations  are^^ 

42  =  E'-5"W'-Ae'- 
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The  A's  stand  for  the  set  of  IP-  and  EA-EOM-CCSD  eigenvectors.     There  is  one 
eigenvector  per  active  orbital.  The  terms  r'^^  stand  for  the  coefficients  in  the  inverse  of 
matrix  constructed  from  the  active  components  of  the  eigenvectors  A.^ 

After  the  S  coefficients  are  determined,  the  problem  remains  of  how  to  determine 
the  G  coefficients.  Since  the  inverse  of  the  normal  ordered  exponential  is  not  known,^^ 
Eq.  (4-9)  cannot  be  used  directly.   Instead  the  equation 

[e']G  =  H[e']  (4-14) 

is  used.  It  has  been  shown^^  that  this  equation  is  equivalent  to 

or 

«=('^pa-(p-'H,         (4-.6) 

where  the  c  means  that  the  equations  must  be  connected. 

In  practice,  all  terms  involving  Si  are  dropped.^^  Since  Si  does  not  change  the  net 
excitation  level,  it  is  not  involved  in  the  decoupling  of  excitation  blocks.  Ignoring  all 
of  the  Si  terms  is  equivalent  to  a  simple  rotation  of  the  eigenvector,  and  since  the  G 
matrix  is  diagonalized  over  all  singly  excited  determinants  this  "rotation"  is  automatically 
compensated  for  in  the  diagonalization.  Therefore,  including  the  Si  terms  does  not  change 
the  final  answer,  and  for  convenience  they  are  left  out.  Leaving  out  Si  also  simplifies 
Eq.  (4-16).  Since  only  the  singles-singles  block  of  G  is  needed,  the  second  term  does 
not  contributed^  leaving  straightforward,  linear  equations  for  G  in  terms  of  H  and  S. 
The  G  equations  in  orbital  form  can  be  found  in  Ref.  66. 
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Finally,  the  G  matrix  is  diagonalized  over  the  space  of  single  excitations  to  get  the 
excited  states.  This  is  actually  the  fastest  step  in  the  process,  since  the  final  matrix  to 
be  diagonalized  is  so  small.  Afterwards,  excited  state  and  transition  properties  can  be 
calculated,  subject  to  certain  approximations.^^ 

Discussion  About  STEOM-CC 

Formally,  the  introduction  of  active  and  inactive  orbitals  is  not  needed  in  STEOM- 
CC.  All  orbitals  could  be  considered  active.  It  would  simply  mean  that  S  would  have 
no  one-body  component  and  that  an  IP-EOM-CCSD  eigenvector  would  need  to  be  found 
for  each  occupied  orbital  and  an  EA-EOM-CCSD  eigenvector  would  be  needed  for  each 
virtual  orbital.  In  practice,  this  would  never  work,  except  for  trivial  cases.  For  example, 
there  is  normally  a  one-to-one  relation  between  IP-EOM-CCSD  states  and  high  lying 
occupied  orbitals,  but  as  the  spectrum  moves  into  the  inner  valence  region,  shake-up 
states  appear.'^  At  that  point,  many  states  appear  per  orbital,  and  even  the  concept 
of  principal  ionization  potential  becomes  unclear.  The  same  will  happen  for  EA-EOM- 
CCSD.  Eventually,  the  iterative  diagonalization  procedure^^  will  have  trouble  converging. 

Even  if  it  were  possible  to  find  an  IP-EOM-CCSD  or  EA-EOM-CCSD  eigenvector 
for  every  orbital,  it  would  not  be  desirable.  The  cost  of  each  IP-EOM-CCSD  state  is 
'"^cc^virt  times  the  number  of  iterations  needed  for  the  diagonalization.  The  cost  of  each 
EA-EOM-CCSD  state  is  r?4^Ar3.^^  per  iteration.  If  the  number  of  states  equals  the  number 
of  orbitals,  then  the  cost  of  calculating  S  becomes  as  large  as  the  cost  of  calculating  T. 
It  is  not  necessary  to  include  all  of  the  orbitals  in  the  active  space,  anyway.  The  second 
similarity  transformation  serves  to  describe  the  differential  dynamic  correlation  between 
the  ground  and  excited  states.    But  relatively  few  orbitals  actually  are  involved  in  the 
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excitation,  and  it  is  only  those  orbitals  involved  in  the  excitation  that  need  to  be  in  the 
active  space. 

The  relationship  of  STEOM-CCSD  to  both  Fock  space  coupled-cluster  (FSCCSD) 
theory^^^^  and  EOM-CCSD^^-  ^^  have  been  discussed  at  some  length.^^  Only  the 
highlights  will  be  reviewed  here.  If  the  active  space  includes  all  orbitals,  then  STEOM- 
CCSD  and  FSCCSD  are  numerically  identical  for  singly  excited  states.^^  For  doubly 
excited  or  higher  states,  the  two  remain  quite  different.  If  the  active  space  does  not 
include  all  orbitals,  then  the  two  are  different  for  singly  excited  states.  The  first  difference 
is  that  STEOM-CCSD  can  be  viewed  as  a  series  of  successive  similarity  transformations, 
whereas  FSCCSD  has  just  one  transformation.  The  second  difference  is  that  in  STEOM- 
CCSD  the  final  diagonalization  is  over  all  singly  excited  determinants,  while  the  final 
diagonalization  in  FSCC  is  over  those  determinants  where  only  the  active  orbitals  change 
occupation. 

To  consider  the  difference  between  STEOM-CCSD  and  EOM-CCSD,  it  is  useful  to 
fully  expand  out  the  STEOM-CCSD  wavefunction.   It  is 


6^(1  +  5  +  ^52  )^1 


oUe^ 


\  (4-17) 


where  Ri  is  the  single  excitation  part  of  the  STEOM-CCSD  eigenvector  and  rg  is  the 
contribution  of  the  reference  determinant.  This  should  be  compared  to  the  EOM-CCSD 
wavefunction,  which  is 

|1'.)  =  e^(7ei  +  7l2)|0)+e^|0)ro.  (4-18) 
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Since  Si  does  not  change  the  excitation  level,  and  is  not  included  in  the  calculation 
anyway,  it  can  be  ignored  in  this  discussion.   It  is  possible  to  directly  equate  terms  in 
Eqs.  (4-17)  and  (4-18).   They  match  up  as  follows: 

e^|0)ro  <^  e^|0)ro 

(4-19) 
e^ S2Ri\0)  ^  e^n2\0) 

The  STEOM-CCSD  term  is  on  the  left  and  the  EOM-CCSD  term  is  on  the  right.  The 
primary  differences  between  the  two  methods  are  therefore  that  the  doubles  term  in 
EOM-CCSD  has  the  flexibility  to  determine  its  optimum  weight,  while  the  doubles  term 
in  STEOM-CCSD  is  a  product  of  a  fixed  Ri  and  a  fixed  52.  On  the  other  hand,  the 
only  triple  excitations  included  in  the  EOM-CCSD  wavefunction  are  of  the  form  e^Tl, 
while  the  STEOM-CCSD  also  includes  a  ^S^Ri  term.  Thus  STEOM-CCSD  has  a  more 
complete  description  of  triple  excitations  but  gives  up  some  flexibility  in  describing  the 
doubles.^6 

The  cost  of  a  STEOM-CC  calculation  in  almost  dominated  by  the  cost  of  calculating 
the  ground  state  CCSD  solution  and  the  ground  state  left  hand  state,  lambda. '°^  Lambda 
is  needed  for  transition  properties.^^  Both  of  these  are  iterative  "o^c-^mr^  ^^^P^-  The  next 
most  expensive  step  is  the  formation  of  H,  which  also  scales  as  'n^occ^tirt  but  is  done  only 
once.  The  costs  of  the  IP-EOM-CCSD  and  EA-EOM-CCSD  calculations  were  discussed 
above.  Calculating  the  5  coefficients  will  scale  as  r?-occA^^,>f A^^,:rf,<icf'  where  Nyirt,act  is 
the  number  of  active  virtual  orbitals.  Forming  G  will  scale  as  n„c.cN^-^fN„irt,act.  The 
final  diagonalization  is  only  an  iterative  Nf,-^^  step. 
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What  all  of  these  scalings  imply  is  that  the  cost  of  the  calculation  is  effectively 
independent  of  the  number  of  excited  states  desired.  In  fact,  the  only  practical  effect  that 
the  number  of  states  has  on  the  cost  is  that  as  higher  lying  excited  states  are  calculated, 
more  orbitals  will  need  to  be  included  in  the  active  space.  This  independence  of  the  cost 
of  the  calculation  to  the  number  of  excited  states  is  in  contrast  to  almost  every  other 
excited  state  method  in  use  today,  and  it  makes  STEOM-CC  well  suited  to  calculating 
large  numbers  of  excited  states. 

However,  STEOM-CC  does  have  some  drawbacks.  The  primary  one  is  that  the 
method  is  only  suitable  for  states  dominated  by  single  excitations.  It  would  be  possible 
to  derive  a  STEOM-CC  variant  whereby  the  final  diagonalization  takes  place  over  the 
space  of  both  singles  and  doubles.  This  would  allow  the  description  of  both  doubly 
excited  states  and  states  which  are  mixed  singles  and  doubles.  The  problem  is  that  the 
equations  become  much  more  complicated.  For  example,  the  second  term  in  Eq.  (4-16) 
would  contribute.  It  would  also  be  at  least  as  expensive  as  an  EOM-CCSD  calculation. 

STEOM-CC  also  has  strong  and  poorly  understood  dependencies  on  the  size  and  the 
nature  of  the  active  space.  It  is  known  that  the  component  of  the  excited  state  vector 
which  lies  within  the  active  space  should  be  at  least  95-99%.  But  it  is  not  clear  how  ,. 
the  error  in  the  energy  and  properties  depends  on  that  percentage.  This  is  an  area  that 
will  require  much  more  practical  experience  before  it  becomes  clear  how  best  to  choose 
the  active  space. 

Finally,  the  calculation  of  excited  state  and  transition  properties  in  STEOM-CC  is 
unsatisfactory.  Several  approximations^^  are  currently  made  in  the  properties  calculations. 
There  is  some  hope,  at  least  for  the  excited  state  properties,  in  that  the  advent  of 
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gradients  for  STEOM-CC  will  give  a  way  to  calculate  one  electron  excited  state  properties 
rigorously. 


CHAPTER  5 
THE  SPECTRUM  OF  FREE  BASE  PORPHIN 

Introduction 

Because  of  their  importance  in  such  biological  processes  as  photosynthesis,  electron 
transfer,  and  oxygen  absorption  and  transport,  the  porphyrins  have  been  extensively 
studied. ''^^  As  the  base  molecule  for  the  porphyrins,  the  electronic  spectrum  of  free 
base  porphin  has  received  much  attention,  both  with  semiemperical  (see,  for  example, 
Ref.  149)  and  ab  initio  methods. i^^- ''5'  151-153,  65  jj^g  interesting  part  of  the  spectrum 
consists  of,  in  order,  two  visible  peaks  known  as  the  Q  bands,  a  very  intense  peak  known 
as  the  B  (or  Soret)  band,  a  shoulder  on  the  B  band,  called  the  N  band,  and  two  other 
small  peaks,  the  L  and  M  bands.  The  spectrum  can  be  found  in  Ref.  154  and  Ref.  153. 

The  traditional  interpretation  of  the  spectrum  is  that,  with  the  molecule  in  the  xy  plane 
and  with  the  two  internal  hydrogens  along  the  jc  axis,  the  lowest  energy  band,  the  Qx  band, 
comes  from  exciting  to  the  I'Bsu  state.  The  Qy  band  then  comes  from  the  1'B2u,  with 
the  B  band  assigned  to  the  2'B3u  and  2'B2u  states.*^°  Nakatsuji  et  al.}^^  based  on  their 
SAC-CI  (symmetry  adapted  cluster-configuration  interaction)  calculations,'^^  reassigned 
the  spectrum.  They  agreed  with  the  assignments  of  the  Q  bands,  but  they  claimed  that  the 
B  band  should  be  assigned  to  only  the  2'B3u  state,  with  the  N  band  being  the  2'B2u  state. 

However,  their  results  have  three  significant  weaknesses.  The  first  is  the  method  used. 
In  principle  SAC-CI  energies  (but  not  oscillator  strengths)  could  be  equivalent  to  EOM- 
CCSD  (equation-of-motion  coupled-cluster  singles  and  doubles)  ^^  and  coupled-cluster 
singles  and  doubles  linear  response^^  energies.  In  practice,  though,  some  non-linear  terms 
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in  the  underlying  ground  state  coupled-cluster  result  are  always  omitted.  Also,  in  these 
calculations  many  double  excitations  were  omitted  based  on  a  perturbation  selection. '^^ 
For  complex  organic  molecules  even  the  untruncated  EOM-CCSD  may  not  be  sufficient. 
In  a  study  of  benzene  and  the  azabenzenes,  EOM-CCSD  had  an  average  error  of  0.32 
eV  for  the  tt^tt*  states. '^^  SAC-CI  should  do  no  better  than  this,  unless  it  has  some 
fortuitous  error  cancellations. 

The  second  criticism  is  the  basis  set  used.  Nakatsuji  et  aO^^  discuss  the  importance 
of  a  rearrangement  to  the  excitation  energies,  yet  their  basis  set  only  had  2s  type  and 
2p  type  functions  on  the  carbons  and  nitrogens,  giving  it  no  flexibility  to  describe  the 
2s  orbitals.  They  also  dropped  some  of  the  occupied  and  virtual  orbitals  corresponding 
to  the  2s  orbitals,  along  with  all  of  the  \s  orbitals  on  the  carbons  and  nitrogens.  The 
result  is  that  the  basis  set  had  limited  flexibility,  no  polarization  functions,  and  no  diffuse 
functions.  The  lack  of  diffuse  and  polarization  functions  corresponds  to  the  conventional 
viewpoint  that,  for  at  least  the  lowest  states  of  free  base  porphin,  polarization  and  diffuse 
functions  are  not  needed. '^^ 

Finally,  there  is  a  problem  with  the  oscillator  strengths.  The  N  band  appears  as  a 
shoulder  to  the  B  band,  but  Nakatsuji  et  a/.  '^^  calculate  the  excitation  to  the  2'B2u  state  to 
have  an  oscillator  strength  sixty-eight  percent  larger  than  the  excitation  to  the  2'B3u  state. 
They  argue  that  the  N  band  is  actually  quite  broad,  with  the  B  band  being  a  narrow  peak 
on  top  of  it.  But  if  that  is  true,  then  the  splitting  between  the  vertical  excitation  energies 
for  the  B  and  N  bands  would  be  less  than  the  0.32  eV  reported. '^"^  To  answer  some  of 
these  questions,  a  series  of  STEOM-CCSD  (similarity  transformed  equation-of-motion 
coupled-cluster  singles  and  doubles)^  excited  state  calculations  have  been  performed. 
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Computational  Details 

STEOM-CC 

In  a  STEOM-CCSD^'  ^  calculation  (i.e.  a  STEOM-CC  calculation^'^  based  on  a 
CCSD^  ground  state),  two  similarity  transformations  are  applied  to  the  second  quantized 
Hamiltonian,  such  that  the  one-  and  two-body  terms  in  the  Hamiltonian  that  increase 
the  excitation  level  are  set  to  zero.  This  effectively  blocks  the  Hamiltonian  matrix  by 
excitation  level,  so  that  the  single  excitations  can  be  accurately  calculated  with  a  diago- 
nalization  over  just  the  singles-singles  block.  The  higher  excitation  parts  of  the  excited 
state  wavefunction  are  then  implicitly  included  through  the  similarity  transformations.  In 
this  way,  STEOM-CC  has  the  same  conceptual  appeal  as  CI  singles,  but  now  in  a  fully 
correlated  structure.  The  STEOM-CCSD  method  has  been  implemented  within  ACES  II. 

The  first  similarity  transformation,  the  same  as  in  EOM-CCSD,^^  involves  the  T 
amplitudes  from  a  ground  state  CCSD^  calculation.  The  new  Hamiltonian  is  then 

H  =  e-^He:^.  (5-1) 

The  next  step  is  to  solve  for  a  set  of  states  with  one  electron  removed  and  one  electron 
added  by  means  of  the  IP-EOM-CCSD  (ionization  potential  EOM-CCSD)^''  and  the  EA- 
EOM-CCSD  (electron  attachment  EOM-CCSD)^^  methods.  These  involve  diagonalizing 
H  over  the  space  of  \h-2h\p  and  \p-2p\h  determinants,  respectively.  These  eigenvectors 
are  used  to  determine  the  S  coefficients  in  the  second  similarity  transformation 

G=p}-//p}.  (5_2) 

The  final  double  similarity  transformed  Hamiltonian  G  is  then  diagonalized  over  the  space 
of  single  excitations.  For  a  detailed  description  of  the  method,  see  Ref.  66. 
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In  the  current  implementation  the  final  wavefunction  in  terms  of  G  has  only  single 
excitations,  limiting  the  method  to  only  singly  excited  states,  but  in  terms  of  the  normal 
Hamiltonian,  the  wavefunction  actually  includes  all  possible  excited  determinants  (with 
only  the  singles  having  optimized  coefficients).  Specifically,  the  wavefunction  contains 
double  excitations  in  the  form  of  S2R1  and  TiRi,  and  triple  excitations  in  the  form  of 
S2R1,  T2R1,  and  T1S2R1,  where  /?i  is  the  STEOM-CC  single  excitation  eigenvector. 

Since  the  higher  excitations  and  the  differential  correlation  between  the  ground 
and  excited  states  are  described  via  the  S  operator,  it  is  critical  that  the  excitation  be 
described  within  the  set  of  active  orbitals  chosen  for  the  IP-EOM-CC  and  EA-EOM-CC 
calculations  and  included  in  S.  In  the  calculations  presented  here,  the  active  components 
of  the  excitations  are  almost  always  above  99%  and  in  all  singlets  are  above  98%.  In 
test  calculations  the  energy  seems  to  be  converged  to  within  0.05  eV  when  the  active 
component  is  98%.^ 

Basis  Set  and  Geometry 

The  basis  set  used  for  the  calculations  presented  here  come  from  the  large  ANO 
basis  set  of  Widmark,  Malmqvist,  and  Roos.^^^  This  is  a  very  large,  generally  contracted 
basis  set.  It  has  14.^9/74^3/ primitives  for  carbon  and  nitrogen  and  Ss4p3d  primitives  for 
hydrogen.  From  this  set,  the  first  3*  and  2p  contractions  were  used  for  C  and  N  and  the 
first  s  contraction  was  selected  for  H.  This  basis  set,  consisting  of  230  functions,  was  the 
same  as  used  by  Merchan  et  al}^'^  It  is  the  same  number  of  contracted  functions  as  that 
used  by  Nooijen  and  Bartlett,^^  but  their  basis  set  had  a  much  smaller  set  of  primitives. 
This  set  is  also  significantly  larger  than  the  one  used  by  Nakatsuji  et  aO^^ 
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This  basis  set  was  then  extended  in  two  different  ways.  First,  to  gauge  the  effect  of 
polarization  functions  on  the  excitation  energies,  the  first  d  contraction  from  the  ANO  set 
was  added  on  the  carbon  and  nitrogen  atoms,  while  the  second  s  function  was  added  to 
the  hydrogens.  This  gave  a  [is2p\d\  set  on  C  and  N  and  a  [2^]  set  on  H,  with  a  total  of 
364  contracted  functions.  In  other  calculations,  to  gauge  the  effect  of  diffuse  functions 
on  the  excited  states,  a  set  of  2s  and  2p  uncontracted  functions  were  added  to  the  center 
of  the  molecule  and  to  the  geometrical  center  of  each  ring.  The  diffuse  exponents  were 
taken  from  Ref.  158,  and  have  been  used  for  naphthalene '^^  and  biphenyl,'^^  where  they 
were  placed  at  the  center  of  the  molecule.  Since  two  functions  of  each  type  were  used, 
the  exact  exponent  chosen  should  not  matter  significantly.  In  all  calculations  the  first 
24  occupied  orbitals,  corresponding  to  the  \s  orbitals  on  carbon  and  nitrogen,  were  left 
uncorrelated.  , 

Two  different  geometries  were  used  in  this  study.  The  first  is  an  idealized  x-ray 
structure,'^  where  the  molecule,  without  the  two  internal  hydrogens,  would  be  of  D4h 
symmetry.  The  hydrogens  make  the  structure  D2h.  It  was  used  in  several  of  the  previous 
studies'^3'  ^^  and  will  be  referred  to  as  "D4h".  With  the  {is2pl\s]  basis  set,  this  geometry 
had  a  SCF  energy  of  -982.955034  Hartrees  and  a  CCSD  energy  of  -985. 154443  Hartrees. 
The  [7>s2p\dl2s\  basis  set  gave  a  SCF  energy  of  -983.416854  Hartrees  and  a  CCSD 
energy  of  -986.676253  Hartrees.  The  basis  with  the  diffuse  functions  gave  a  SCF  energy 
of  -982.967070  Hartrees  and  a  CCSD  energy  of  -985.182042  Hartrees. 

The  other  geometry  is  a  B3PW91/6-31G*  optimized  geometry,'^'  which  will  be  re- 
ferred to  as  "Opt".  It  also  has  D2h  symmetry.  This  geometry  gave  a  SCF  energy  of 
-982.964858  Hartrees  and  a  CCSD  energy  of  -985.163441  Hartrees  with  the  [?>s2pl\s] 
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basis  set.  The  diffuse  basis  gave  a  SCF  energy  of  -982.977085  Hartrees  and  a  CCSD  en- 
ergy of -985. 19 1231  Hartrees.  The  [3^2^1^/25]  basis  gave  a  SCF  energy  of -983.430681 
Hartrees  and  a  CCSD  energy  of -986.688179  Hartrees,  putting  this  geometry  7.5  kcal/mol 
below  the  D4h  geometry. 

Results 
Ionized  and  Electron  Attached  States 

In  all  of  the  current  calculations,  the  occupied  part  of  the  active  space  consisted  of 
eleven  orbitals.  Therefore,  eleven  IP-EOM-CCSD  states  were  calculated.  In  Table  5-1 
the  current  results,  with  the  three  basis  sets,  are  compared  with  the  SAC-CI  results, '^^ 
with  the  previous  STEOM-CC  results,^^  and  with  experiment. '^^  Because  of  the  low 
resolution  of  the  experimental  spectrum,  it  is  very  difficult  to  relate  the  measured  peaks 
to  the  calculated  states.  Therefore,  the  assignments  given  should  be  considered  tentative. 

While  adding  the  diffuse  functions  does  little  to  change  the  ionization  potentials, 
every  other  time  that  the  basis  set  was  enlarged,  from  the  previous  STEOM-CC  results,^^ 
to  the  DZ  (the  [3s2p/\s]  basis)  results,  to  the  polarized  (the  [3s2p\d/2s]  basis)  results,  the 
electron  became  more  bound.  Adding  the  polarization  functions  had  the  most  dramatic 
effect.  It  caused  the  first  two  states  to  switch  and  the  sixth  and  seventh  states  to  switch. 
However,  the  differences  are  close  to  an  order  of  magnitude  smaller  than  the  errors  in 
the  calculations,  so  nothing  can  be  said  definitively.  Going  from  the  "D4h"  geometry 
to  the  "Opt"  geometry  had  effects  of  less  than  0.2  eV.  The  "Opt"  geometry  typically 
had  larger  IP's. 

Twenty-three  virtual  orbitals  (4  ag,  2  bjg,  3  bjg,  3  bjg,  3  au,  2  bj^,  3  bau,  and  3  bgu) 
were  included  in  the  active  space  for  the  DZ  and  the  polarized  basis  sets.  For  the  diffuse 
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basis  set  several  more  orbitals  had  to  be  included.  The  diffuse  basis  added  forty  Rydberg 
orbitals,  and  most  of  them  had  orbital  eigenvalues  very  close  to  zero.  Therefore,  to 
include  all  of  the  equivalent  orbitals  that  were  included  for  the  DZ  basis,  the  active  space 
needed  to  consist  of  fifty-nine  virtual  orbitals.  But  that  many  orbitals  caused  convergence 
problems.  The  final  set  included  in  the  active  space  consisted  of  fifty-one  virtual  orbitals 
(10  ag,  4  big,  5  b2g,  5  bsg,  3  au,  8  bju,  5  biu,  and  5  bsu).  If  an  IP-EOM-CC  or  an 
EA-EOM-CC  eigenvector  is  less  than  70%  singles,  the  code  automatically  excludes  it 
from  the  second  similarity  transformation.  For  the  diffuse  basis,  one  of  the  Bsg  states 
was  excluded,  leaving  fifty  states  for  the  similarity  transformation. 

In  each  of  the  calculations  two  positive  electron  affinities  were  predicted.  They  are 
listed  in  Table  5-2.  The  other  electron  attached  states,  even  though  they  do  not  correspond 
to  stable  states,  are  still  essential  for  the  calculation  of  the  excited  states.  They  go  into 
the  second  similarity  transformation  and  help  describe  the  differential  correlation  between 
the  ground  and  excited  states. 

Excited  States 

The  singlet  excited  states  of  free  base  porphin  are  listed  in  Tables  5-3  to  5-7,  and  the 
triplet  excited  states  are  in  Tables  5-8  to  5-12.  The  important  states  will  be  discussed 
later,  but  some  general  comments  can  be  made.  The  first  is  that  the  Rydberg  states 
start  about  4.5  eV,  which  is  right  in  the  region  of  most  interest.  The  second  is  that  the 
addition  of  the  diffuse  functions  never  changes  the  energy  of  the  valence  states  by  more 
than  0.02  eV.  This  has  several  consequences.  It  implies  that  the  DZ  basis  set  has  diffuse 
enough  tails  that  the  diffuse  functions  are  not  needed  to  describe  the  valence  region. 
It  also  means  that  there  is  essentially  no  mixing  between  Rydberg  and  valence  excited 


90 
states.  Finally,  it  suggests  that  adding  diffuse  functions  to  the  polarized  basis  would  not 
significantly  change  the  energetics.   Adding  diffuse  functions  might,  however,  have  an 
effect  on  the  oscillator  strengths.  Those  tend  to  change  a  little  with  the  addition  of  the 
diffuse  functions. 

In  general,  going  from  the  "D4h"  geometry  to  the  "Opt"  geometry  has  little  effect  on 
the  spectrum.  The  excitation  energy  typically  increases  by  about  0. 1  eV.  The  important 
exceptions  are  the  'B2u  state  which  drops  from  5.22  eV  to  5.07  eV  and  the  'Bsu  state 
which  drops  from  5.33  eV  to  5.23  eV. 

A  note  should  be  made  about  the  oscillator  strengths.  The  oscillator  strengths  are 
calculated  as 

/  =  ^a;,(*,|//|^,)(^,|//|1',),  (5-3) 

where  oj  is  the  excitation  energy.  The  right-hand  ground  state  is  the  coupled-cluster 
wavefunction,  the  left-hand  ground  state  is  the  lambda  solution  from  coupled-cluster 
theory, '^^  and  the  right-hand  excited  state  is  the  STEOM-CC  state.  Currently,  several 
approximations  are  introduced  when  calculating  the  left-hand  excited  state.^  These 
sometimes  can  cause  numerical  problems  with  the  properties  calculated.  When  the 
oscillator  strength  is  listed  as  ( — ),  it  means  that  the  approximations  are  too  severe, 
and  the  calculated  oscillator  strength  is  unreliable.  Since  the  energies  are  calculated 
with  the  right  hand  excited  state  wavefunction,  they  are  still  correct.  For  the  polarized 
basis  results,  another  approximation  is  made.  Because  the  cost  of  calculating  the  ground 
state  lambda  vector  is  prohibitive  for  that  large  a  basis  set,  the  left-hand  ground  state  is 
estimated  with  the  same  approximations  as  the  excited  state  left-hand  wavefunction.  In 
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a  test  calculation  with  the  DZ  basis,  this  makes  a  difference  in  the  oscillator  strengths 
of  less  than  20%  in  every  case. 

The  "D4h"  geometry,  polarized  basis  results  for  the  singlet  valence  states  from  Tables 
5-3  to  5-7  are  listed  in  Table  5-13,  along  with  CASPTl,^^^  SAC-CI,  '^^  previous 
STE0M-CC,^5  and  EOM-CCSD  results.  The  EOM-CCSD  calculation  is  at  the  "D4h" 
geometry  with  the  DZ  basis.  The  valence  states  are  numbered  for  convenience,  but  these 
numbers  are  only  accurate  through  3^Big.  After  that  the  Rydberg  states  should  enter 
into  the  numbering. 

The  energy  of  the  first  triplet  state  has  been  measured  in  a  solvent/ethyl  iodide  mixture 
at  77K.'^^  The  phosphorescence  peak  is  at  1.58  eV,  well  above  the  STEOM-CC  polarized 
basis  result  of  1.26  eV.  The  CASPT2  result  is  between  the  two  at  1.37  eV.'^^ 

Discussion 

The  recent  controversy  over  the  assignment  of  the  spectrum  centers  around  how  to 
assign  the  N  band.  The  argument  of  Nakatsuji  etal.}^^  essentially,  is  that  they  calculated 
no  other  optically  allowed  states  in  the  3^  eV  range,  and  therefore,  by  default,  the  N 
band  must  be  2'B2u.  Table  5-13  shows  the  poor  quality  of  their  results.  Had  they  used  a 
decent  basis  set  and  had  they  not  made  the  approximations  that  are  always  used  in  their 
SAC-CI  calculations,'^^  their  calculated  energies  would  have  approached  the  EOM-CCSD 
excitation  energies,  meaning  that  they  all  would  have  increased.  At  that  point,  it  becomes 
impossible  to  draw  any  meaningful  conclusions  about  the  assignment  of  the  spectrum. 

The  CASPT2'^2  results  for  the  B  band  are  too  low.  This  could  have  been  caused  by 
the  size  of  their  active  space.  It  has  been  shown  several  times'^'*-  '^''  ^^^-  ^^  that  for  the 
2'B3u  and  2'B2u  states,  Gouterman's  four  orbital  model'^^-'^"^  is  not  sufficient;  the  4biu 
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orbital  must  also  be  included.  Merchan  et  al}^^  did  not  include  it  in  their  active  space. 
More  CASPT2  calculations,  with  a  larger  active  space  and  with  d  functions,  would  be 
very  informative. 

The  Qx  and  Qy  bands  belong  to  the  l^Bsu  and  1'B2u  states.  The  current  results  show 
fortuitously  good  agreement  with  the  Qy  band.  For  the  Qx  band  all  of  the  methods  predict 
too  low  an  excitation  energy.  A  fluorescence  spectrum  of  free  base  porphin  taken  in  a 
supersonic  jet  expansion '^^  placed  the  0-0  transition  for  the  Qx  band  at  2.0234  eV  and 
the  0-0  transition  for  the  Qy  band  at  2.4653  eV,  slightly  higher  than  the  vertical  excitation 
energies  reported  by  Edwards  et  aO^'^  The  addition  of  polarization  functions  increased 
the  Qx  excitation  energy  by  0.05  eV,  so  the  use  of  much  larger  basis  sets  may  further 
increase  the  calculated  excitation  energy,  moving  it  towards  the  experimental  number. 

The  current  results  strongly  support  the  original  assignment  for  the  B  band  being 
both  the  2'B3u  and  2'B2u  states.  The  N  band  is  then  assigned  to  the  3'B3u  state.  If 
the  assignments  of  Nakatsuji  et  al.  were  correct,  it  would  mean  that  the  polarized  basis 
STEOM-CC  energies  would  have  to  be  0.19  eV  too  low  for  the  first  L  peak  and  0.32 
eV  too  low  for  the  second  L  peak.  It  is  highly  unlikely  that  the  STEOM-CC  would 
consistently  underestimate  the  excitation  energy  like  that.  Instead,  this  assignment  puts 
2'B3u  0.14  eV  above  the  B  peak  and  2'B2u  0.29  eV  above  the  B  peak.  The  3'B3u  state 
is  0.41  eV  above  the  N  peak.  This  error  is  uncomfortably  large.  However,  the  energies 
of  all  of  these  states  dropped  substantially  when  the  polarizations  functions  were  added. 
Larger  basis  set  calculations  should  further  decrease  these  gaps. 

The  intensities  still  present  a  problem.  The  experimental  oscillator  strength  of  the 
N  band  is  less  than  0.1,  and  the  calculated  oscillator  strength  is  0.93.  Edwards  et  aO^'^ 
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mention  that  the  shape  of  the  B  and  N  bands  would  be  consistent  with  the  two  states  of 
the  B  bands  being  split  by  1500  cm"'  and  an  intense  N  band  donating  intensity  into  the 
B  band.  These  calculations  give  a  splitting  of  1200  cm~'  between  the  2'B3u  and  2'B2u 
states.  Thus  these  calculations  agree  with  the  possible  interpretation  given  by  Edwards 
et  al}^^  A  B  band  splitting  of  240  cm"'  was  measured  in  a  low  temperature  crystal 
spectrum  and  assigned  to  the  energy  difference  between  the  two  electronic  states, '^^  but 
this  assignment  has  been  disputed. '^^ 

These  assignments  then  leave  the  two  peaks  of  the  L  band  assigned  to  the  3'B2u 
and  4'B2u  states.  The  diffuse  basis  energy  for  the  3'B2u  state  agrees  well  with  the  first 
L  peak,  but  the  4'B2u  state  is  0.33  eV  higher  than  the  second  L  peak.  This  difference 
should  also  be  reduced  with  larger  basis  sets.  The  experimental  oscillator  strength  for  the 
two  states  combined  is  about  0.1.  The  calculated  oscillator  strengths  are  well  above  that. 
The  I'Biu  state  sits  under  the  L  peak,  but  its  intensity  is  so  small  that  it  is  not  visible. 

The  one  remaining  problem  with  this  assignment  is  the  M  peak.  The  calculated 
excitation  energy  for  the  4'B3u  state  is  5. 17  eV,  0.33  eV  lower  than  the  M  band  maximum. 
The  intensities  match  well,  though,  and  there  is  no  other  assignment  that  makes  sense. 
To  say  that  4'B3u  is  part  of  the  L  band  would  require  that  the  excitation  energy  to  drop 
by  0.5  eV,  and  it  would  cause  problems  assigning  the  other  states.  Still,  it  is  not  clear 
why  the  calculated  excitation  energy  would  be  so  low. 

Finally,  a  point  should  be  made  about  the  cost  of  these  calculations.  For  the  polarized 
basis  calculation  there  are  57  occupied  and  283  virtual  functions.  This  gave  2479  T^  and 
34  170  895  7^2  amplitudes.  The  time  taken  for  each  step  on  a  Cray  C90  is  listed  in  Table 
5-14.  It  took  only  36  seconds  per  excited  state  to  calculate  the  energies  and  properties. 
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Most  of  that  time  was  used  calculating  the  properties.  The  energies  took  1 .4  seconds 
each.  This  is  a  very  important  advantage  of  the  STEOM-CC  method.  The  total  time  for 
the  calculation  is  effectively  independent  of  the  number  of  excited  states  calculated.  That 
makes  it  possible  to  study  entire  spectra  instead  of  just  a  few  selected  states. 

Conclusions 

STEOM-CCSD^  calculation  on  free  base  porphin  using  a  [3s2p\d/2s]  basis  set  are 
presented.  These  are  the  first  reported  calculations  on  excited  states  above  the  lowest 
that  use  polarization  functions.  For  the  important  optically  allowed  excited  states,  the 
polarization  functions  have  a  significant  effect.  These  calculations  strongly  support  the 
traditional  interpretation  of  the  spectrum,  in  that  the  intense  B  band  is  assigned  to  both 
the  2'B3u  and  2'B2u  states. 

This  study  is  possible  only  because  the  second  similarity  transformation  in  STEOM- 
CCSD  means  that  the  excited  states  can  be  represented  in  terms  of  single  excitations 
with  respect  to  a  highly  modified  Hamiltonian.  Thus  it  is  possible  to  efficiently  calculate 
many  excited  states  at  once. 
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Table  5-1 :  Ionization  potentials  (in  eV)  for  free  base  porphin. 


SAC- 
CP 

STEOM- 

L'4h 

"Opt" 

DZ 

diffuse 

polarized 

DZ 

diffuse 

Expt.'^ 

l'Au(7r) 

6.35 

6.34 

6.65 

6.63 

6.66 

6.71 

6.69 

6.9 

l'B„(7r) 

6.09 

6.16 

6.46 

6.43 

6.68 

6.48 

6.46 

7.2 

l^B3g(7r) 

7.98 

7.84 

8.11 

8.09 

8.19 

8.29 

8.27 

8.4 

22B,u(7r) 

8.21 

8.11 

8.38 

8.36 

8.43 

8.53 

8.51 

8.8 

l2Ag(«) 

8.65 

8.50 

8.78 

8.76 

9.18 

8.96 

8.94 

9.1 

l^B2g(7r) 

9.10 

8.79 

9.08 

9.06 

9.20 

9.01 

8.99 

l2B2u(n) 

8.79 

8.57 

8.84 

8.82 

9.23 

8.99 

8.96 

32Biu(7r) 

9.34 

9.10 

9.37 

9.36 

9.47 

9.22 

9.20 

22B2g(7r) 

9.42 

9.22 

9.50 

9.47 

9.55 

9.43 

9.41 

22B3g(7r) 

9.54 

9.31 

9.59 

9.57 

9.64 

9.66 

9.63 

22Au(7r) 

10.85 

10.43 

10.71 

10.69 

10.72 

10.76 

10.74 

«)Ref. 

153. 

^>Ref. 

65. 

,  ■"' ,_  ■-•"    ' 

'=>Ref. 

162. 
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Table  5-2: 

Electron  affinities  (in 

;V)  for  free  base  porphin. 

IP-EOM- 

"D4h" 

"Opt" 

DZ 

diffuse 

polarized 

DZ 

diffuse 

l'B2g 

l'B3g 

0.50 
0.37 

0.83 
0.70 

0.81 
0.67 

1.07 
0.99 

0.80 
0.78 

0.78 
0.76 

^)Ref.  65. 
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Table  5- 
ineV. 

-3:  Singlet  excited  states  of  free  base  porphin  through  4.89  eV.  Energies  are 
Oscillator  strengths  are  given  in  parentheses  below  the  excitation  energies. 

L'4h 

"Opt" 

DZ 

diffuse 

polarized 

DZ 

diffuse 

'B3u  (V) 

1.70 

1.70 

1.75 

1.71 

1.71 

(-) 

(— ) 

(0.0007) 

(0.0003) 

(0.0004) 

'B2u  (V) 

2.59 
(0.018) 

2.59 
(0.017) 

2.40 
(0.013) 

2.65 
(0.007) 

2.65 
(0.006) 

'B,g  (V) 

3.57 

3.56 

3.44 

3.77 

3.78 

'B3u  (V) 

3.64 
(0.982) 

3.63 
(0.981) 

3.47 
(0.693) 

3.71 
(1.21) 

3.71 
(1.23) 

'B2u  (V) 

3.74 
(1.36) 

3.74 
(1.37) 

3.62 
(1.20) 

3.78 
(1.46) 

3.77 
(1.48) 

'Ag  (V) 

4.04 

4.04 

3.95 

4.09 

4.10 

'B2g  (V) 

4.05 

4.05 

4.21 

4.25 

4.25 

'Au  (V) 

4.09 

4.09 

4.24 

4.26 

4.26      - 

^B3u  (V) 

4.22 
(0.721) 

4.22 
(0.741) 

4.06 
(0.931) 

4.34 
(0.494) 

4.35 
(0.501) 

'Bsg  (V) 

4.43 

4.41 

4.51 

4.55 

4.54 

'B,g  (V) 

4.44 

4.44 

4.53 

4.41 

4.41 

^Biu(V) 

4.50 
(0.004) 

4.48 
(0.003) 

4.56 
(0.002) 

4.59 
(0.004) 

4.58 
(0.005) 

'B,u  (R) 

7.18 
(0.008) 

4.51 
(0.006) 

5.48 
(0.013) 

7.23 
(0.006) 

4.54 
(0.004) 

'B2u  (V) 

4.63 
(0.441) 

4.63 
(0.444) 

4.35 
(0.422) 

4.65 
(0.239) 

4.66 
(0.241) 

'Ag  (V) 

4.67 

4.66 

4.46 

4.72 

4.72 

'Au  (R) 

7.33 

4.72 

5.52 

7.40 

4.78 

'B2g  (R) 

6.63 

4.74 

5.89 

6.63 

4.78 

'B,g  (V) 

4.79 

4.78 

4.66 

4.67 

4.67 

'Bsg  (R) 

6.78 

4.82 

6.04 

6.80 

4.84 

'Ag  (R) 

4.89 

4.94 

98 


Table  5-4:  Singlet  excited  states  of  free  base  porphin  from  4.94  eV  to  5.50  eV.  Energies  are 
in  eV.  Oscillator  strengths  are  given  in  parentheses  below  the  excitation  energies. 


"D4h" 


"Opt" 


DZ 


diffuse 


polarized 


DZ 


diffuse 


'B3g  (R) 

6.83 

4.94 

5.91 

6.87 

5.01 

'Ag  (V) 

5.00 

4.99 

4.87 

5.01 

5.02 

'B2g  (R) 

6.97 

5.01 

6.08 

7.03 

5.06 

'Au  (R) 

7.14 

5.07 

6.62 

7.14 

5.09 

'Big  (R) 

5.13 

5.19 

'Biu  (R) 

5.14 

6.89 

7.93 

5.17 

(0.0006) 

(0.000003) 

(0.002) 

(0.0005) 

'Biu  (R) 

5.20 
(0.001) 

5.22 
(0.0008) 

'B2u  (V) 

5.22 

5.22 

5.00 

5.06 

5.07 

(0.172) 

(0.182) 

(0.153) 

(0.290) 

(0.312) 

'B3u  (R) 

5.23 
(0.007) 

5.28 
(0.010) 

'Ag  (V) 

5.28 

5.25 

5.12 

5.26 

5.27 

'Big  (V) 

5.26 

5.26 

5.04 

5.25 

5.25 

'Biu  (R) 

5.28 
(0.000008) 

7.37 
(0.00002) 

5.33 
(0.00002) 

'B2u  (R) 

5.31 
(-) 

6.65 
(0.075) 

5.32 
(0.002) 

'Au  (R) 

8.19 

5.33 

6.78 

8.22 

5.40 

'B3u  (V) 

5.33 

5.33 

5.17 

5.22 

5.23 

(0.400) 

(0.387) 

(0.272) 

(0.471) 

(0.476) 

'Au  (R) 

5.39 

5.45 

'B2u  (R) 

5.45 
(0.002) 

5.52 
(0.001) 

'B2g  (R) 

7.84 

5.46 

7.43 

7.90 

5.48 

'Bsg  (R) 

5.47 

7.91 

5.50 

'Bbu  (R) 

5.50 
(— ) 

5.55 
(— ) 
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Table  5-5:  Singlet  excited  states  of  free  base  porphin  from  5.51  eV  to  5.80  eV.  Energies  are 
in  eV.  Oscillator  strengths  are  given  in  parentheses  below  the  excitation  energies. 


DZ 


"D4h" 


"Opt" 


diffuse 


polarized 


DZ 


diffuse 


'Biu  (R) 

'B2g  (R) 
'B3g  (R) 
'Ag  (R) 
'Ag  (R) 
'B3g  (R) 
'B2g  (R) 
'B2g  (R) 
'B3g  (R) 
'Ag  (R) 
'Au  (R) 
'B3g  (R) 

'Biu  (R) 

'Ag  (R) 
'Au  (R) 

'B,u  (R) 

'B2g  (R) 

'B3u  (R) 

'B,g  (R) 

'2u  (R) 


5.71 


'B 


5.51 
(0.001) 

5.53 

5.56 

5.57 

5.61 

5.65 

5.66 

5.66 

5.68 

5.69 

5.71 

5.72 

5.72 
(0.00004) 

5.72 

5.73 

5.75 
(0.001) 

5.75 

5.76 
(0.002) 

5.78 

5.80 
(0.002) 


7.50 


5.51 


5.68 


5.54 
(0.001) 

5.55 

5.58 

5.60 

5.64 

5.71 

5.68 

5.72 

5.70 

5.72 

5.77     :. 

5.78 

5.74 
(0.00007) 

5.68 

5.75 

5.78 
(0.0007) 

5.81 

5.80 
(0.0005) 

5.84 

5.82 
(0.002) 
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Table  5-6: 

Singlet  excited  states  of  free  base  porphin  from  5.80  eV  to  6.01  eV.  Energies  are 

ineV. 

Oscillator  strengths 

are  given  in 

parentheses  below  the  excitation  energies. 

"D4h" 

"Opt" 

DZ 

diffuse 

polarized 

DZ                diffuse 

'Big  (R) 

5.80 

5.85 

'Au  (R) 

5.82 

5.84 

'B3g  (R) 

5.85 

5.92 

'B2g  (R) 

5.87 

5.93 

'Big  (R) 

5.88 

5.88 

5.59 

5.88                 5.90 

'B2g  (R) 

5.90 

5.92 

'Big  (R) 

5.90 

5.96 

'B3g  (R) 

5.91 

5.93 

'B3g  (R) 

5.91 

5.95 

'B2g  (R) 

5.92 

5.94 

'Au  (R) 

5.92 

5.99 

'Biu  (R) 

5.93 
(0.00005) 

5.99 
(0.000002) 

'Au  (R) 

5.94 

6.00 

'Ag  (R) 

5.95 

5.98 

'B2u  (R) 

5.97 
(0.001) 

6.03 
':                (0.0008) 

'Bsg  (R) 

7.89 

5.99 

6.73 

■  V              6.10 

'Bsu  (R) 

5.99 
(— ) 

6.05 

(— ) 

'Biu  (R) 

7.94 

6.00 

6.70 

6.07 

(0.091) 

(0.006) 

(0.021) 

(0.001) 

'Ag  (R) 

6.01 

6.03 

'B,u  (R) 

7.81 

6.01 

6.72 

6.39 

(0.003) 

(0.006) 

(0.020) 

(0.006) 
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Table  5-7:  Singlet  excited  states  of  free  base  porphin  from  6.04  eV  to  6.23  eV.  Energies  are 
in  eV.  Oscillator  strengths  are  given  in  parentheses  below  the  excitation  energies. 


"D4h" 

"Opt" 

DZ 

diffuse 

polarized 

DZ                diffuse 

'Blu  (R) 

6.04 
(0.001) 

6.06 
(0.005) 

'B2g  (R) 

6.05 

6.08 

'B3g  (R) 

6.08 

'B3g  (R) 

6.11 

'Biu  (R) 

6.15 
(0.018) 

'Bag  (R) 

6.09 

6.16 

'B2g  (R) 

6.10 

6.16 

'B,u  (R) 

6.10 
(0.0003) 

6.12 
(0.005) 

'B2g  (R) 

6.13 

6.19 

'Au  (R) 

6.14 

6.17 

'Big  (R) 

6.17 

6.24 

'B,u  (R) 

6.19 
(0.022) 

6.21 
(0.023) 

'Ag  (R) 

6.19 

6.28 

'Big  (R) 

6.20 

6.25 

'B3u  (V) 

6.20 

6.20 

6.07 

6.33 

(0.166) 

(0.160) 

(0.182) 

(0.110) 

'Au  (R) 

6.23 

6.29 
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Table  5-8: 

Triplet  excited  states  of  free  base  porphin  through 

1  4.73  eV.  Energies 

are  in  eV. 

"D4h" 

"Opt" 

DZ 

diffuse 

polarized 

DZ 

diffuse 

^B2u  (V) 

1.19 

1.19 

1.26 

1.15 

1.15 

^Bsu  (V) 

1.82 

1.83 

1.80 

1.87 

1.88 

^B2u  (V) 

2.02 

2.02 

Lg5 

2.07 

2.07 

^B3u  (V) 

2.29 

2.28 

1.98 

2.31 

2.31 

^Big  (V) 

2.96 

2.97 

2.87 

3.09 

3.10 

'B3u  (V) 

3.10 

3.11 

2.98 

3.21 

3.21 

3Ag  (V) 

3.49 

3.48 

3.33 

3.61 

3.61 

^B,g  (V) 

3.50 

3.50 

3.45 

3.47 

3.47 

'Ag  (V) 

3.81 

3.81 

3.64 

3.82 

3.82 

'Big  (V) 

4.02 

4.02 

3.93 

3.83 

3.84 

'B2g  (V) 

4.04 

4.04 

4.19 

4.24 

4.24 

'B2u  (V) 

4.08 

4.08 

3.85 

3.90 

3.90 

'Au  (V) 

4.10 

4.10 

4.23 

4.27 

4.27 

'B3g  (V) 

4.15 

4.14 

4.29 

4.27 

4.26 

'B,u  (V) 

4.20 

4.18 

4.33 

4.28 

4.27 

'B2U  (V) 

4.23 

4.24 

4.07 

4.29 

4.31 

'Ag  (V) 

4.48 

4.47 

4.31 

4.47 

4.47 

'Biu  (R) 

7.17 

4.50 

5.43 

7.22 

4.52 

-'Au  (R) 

7.29 

4.71 

5.49 

7.36 

4.77 

3B2g  (R) 

6.57 

4.73 

5.85 

6.58 

4.77 
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Table  5-9:  Triplet  excited  states  of  free  base  porphin  from  4.74  eV  to  5.47  eV.  Energies  are  in  eV. 


"D4h" 

"Opt" 

DZ 

diffuse 

polarized 

DZ 

diffuse 

^Big  (V) 

4.76 

4.74 

4.57 

4.82 

4.82 

^B3g  (R) 

6.71 

4.81 

5.98 

6.72 

4.83 

3Ag  (R) 

4.90 

4.94 

^B3u  (V) 

4.92 

4.92 

4.74 

4.83 

4.83 

^B3g  (R) 

6.80 

4.93 

5.89 

6.84 

..      5.01 

3Ag  (V) 

4.98 

4.97 

4.76 

4.97 

4.97 

3B2g  (R) 

6.50 

5.01 

6.05 

7.00 

5.06 

^Au  (R) 

7.02 

5.06 

6.53 

7.02 

5.08 

3Big  (R) 

5.10 

5.18 

^B,u  (R) 

5.13 

6.87 

7.79 

5.17 

^B,u  (R) 

5.20 

5.22 

^Big  (V) 

5.24 

5.24 

5.23 

^Bsu  (R) 

5.24 

5.27 

5.28 

^B,u  (R) 

7.34 

5.28 

6.70 

7.37 

5.33 

^B2u  (R) 

5.30 

5.32 

^Au  (R) 

5.32 

6.76 

8.19 

5.39 

^B2u  (R) 

5.39 

5.28 

5.50 

^Au  (R) 

5.39 

5.45 

3B2g  (R) 

7.74 

5.45 

7.80 

5.48 

^B3g  (R) 

5.47 

5.49 

._  J. . 
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Table  5-10:  Triplet  excited  states  of  free  base  porphin 
from  5.49  eV  to  5.75  eV.  Energies  are  in  eV. 


DZ 


"D4h" 


"Opt" 


diffuse 


polarized 


DZ 


diffuse 


^B3u  (R) 

^B3u  (R) 

^B,u  (R) 
^B2g  (R) 
3B3g  (R) 
^B2u  (R) 

^Ag  (RA^) 
3Ag  (R) 
3Ag  (R) 
^Bsg  (R) 
^B2g  (R) 
'B2g  (R) 
3B3g  (R) 
3Ag  (R) 
^Au  (R) 
^B3g  (R) 
^B,u  (R) 
^Au  (R) 
^B2g  (R) 
^Biu  (R) 


5.51 


5.52 


5.49 
5.51 
5.51 
5.52 
5.55 
5.56 
5.56 
5.57 
5.64 
5.65 
5.66 
5.66 
5.68 
5.70 
5.71 
5.71 
5.71 
5.72 
5.75 
5.75 


5.26 


5.55 

5.57 

5.63 

5.53 

5.55 

5.57 

5.59 

5.66 

5.59 

5.57 

5.60 

5.66 

8.09 

5.71 

5.68 

5.72 

5.70 

5.72 

5.77 

5.78 

5.74 

5.75 

5.80 

5.78 
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Table  5-1 1 :  Triplet  excited  states  of  free  base  porphin 
from  5.77  eV  to  5.97  eV.  Energies  are  in  eV. 


^B3u  (R) 

^B,g  (R) 

^B2u  (R) 
^B,g  (R) 
^A,  (R) 
^B3g  (R) 
^B2g  (R) 
^B,g  (R) 
^B3g  (R) 
^B2g  (R) 
3B3g  (R) 
^B2g  (R) 
^Au  (R) 
^B,u  (R) 
^Au  (R) 
'B,u  (R) 
^Big  (R) 
^B2u  (R) 
^B3g  (R) 
^Ag  (R) 


DZ 


7.75 


lip.    II 

diffuse 
5.77 
5.77 
5.80 
5.80 
5.82 
5.85 
5.87 
5.88 
5.89 
5.89 
5.91 
5.92 
5.92 
5.93 
5.94 
5.96 
5.96 
5.96 
5.96 
5.97 


polarized 


DZ 


6.64 


6.63 


'Opt" 


diffuse 
5.80 
5.84 
5.82 
5.87 
5.84 
5.92 
5.93 
5.95 
6.07 
5.92 
5.93 
5.94 
5.98 
5.99 
6.00 
6.09 
6.06 
6.03 
5.94 
5.99 
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Table  5-12:  Triplet  excited  states  of  free  base  porphin 

from  5.97  eV  to  6.11  eV.  Energies  are  in  eV. 

"D4h"                                 I  "Opt" 

DZ                diffuse            polarized  DZ                diffuse 

3B3u  (R)             5.97                 5.97                  5.77  6.05                 6.07 

3B3u  (R)                                    6.00  6.04 

^Ag  (R)                                     6.00  6.03 

^Biu  (R)                                    6.01  6.06 

^Biu  (R)                                    6.04  6.06 

3B2g  (R)                                   6.05  6.08 

3B3g  (R)                                    6.08  6.11 

^Bsg  (R)                                    6.09  6.15 

3B2g  (R)                                    6.10  6.16 

^Biu  (R)                                    6.10  6.13 

^B3g  (R) 6AI 6.31 


.»,v ' 
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Table  5-13:  Singlet  valence  excited  states  for  free  base  porphin.  Energies  are  in 
eV.  Oscillator  strengths  given  in  parentheses  below  the  excitation  energies. 


CASP12^ 

SAC-CI'' 

EOM- 
CCSD 

STEOM- 

"D4h" 

polarized 

Expt.'' 

I'Bsu 

1.70 

1.75 

2.15 

1.72 

1.75 

1.98  Qx 

(0.001) 

(0.0001) 

(0.0007) 

(— ) 

(0.0007) 

(0.01) 

l'B2u 

2.26 

2.23 

2.76 

2.61 

2.40 

2.42  Qy 

(0.016) 

(0.0006) 

(0.007) 

(0.016) 

(0.013) 

(0.06) 

I'Big 

3.55 

3.63 

3.44 

2'B3u 

2.91 

3.56 

3.89 

3.66 

3.47 

(1.66) 

(1.03) 

(1.00) 

(1.03) 

(0.693) 

3.33  B 

2'B2u 

3.04 

3.75 

4.05 

3.77 

3.62 

(1.15) 

(1.54) 

(1.73) 

(1.56) 

(1.42) 

(1.20) 

2'Ag 

4.25 

4.08 

3.95 

3'B3u 

4.24 

4.51 

4.28 

4.06 

3.65  N 

(0.976) 

(0.812) 

(0.71) 

(0.931) 

(<o.i) 

l'B2g 

4.05 

4.08 

4.21 

I'Au 

4.18 

4.14 

4.24 

3*B2u 

4.52 

4.83 

4.67 

4.35 

4.25  L 

(0.350) 

(0.313) 

(0.44) 

(0.422) 

(~0.1) 

l'B3g 

4.37 

4.45 

4.51      . 

2'B,g 

4.62 

4.50 

4.53 

I'Biu 

4.51 

4.82 

4.54 

4.56 

(0.0005) 

(0.004) 

(0.004) 

(0.002) 

,      3'Ag 

4.74 

4.72 

4.46 

3iB,g 

5.13 

4.82 

4.66 

4iAg 

5.28 

5.06 

4.87 

•    4'B2u 

5.31 

5.26 

5.00 

4.67  L 

(   \    X;  V? 

(0.280) 

(0.20) 

(0.153) 

CO.l) 

5'Ag 

5.32 

5.12 

4'B,g 

5.31 

5.04 
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Table  5-13:  (Continued)      Singlet  valence  excited  states  for  free  base  porphin.  Energies 
are  in  eV.  Oscillator  strengths  given  in  parentheses  below  the  excitation  energies. 


CASPT2^ 

SAC-CI'' 

EOM- 
CCSD 

STEOM- 

"D4h" 

polarized 

Expt.'' 

4'B3u 

5.45 

5.38 

5.17 

5.50  M 

(0.351) 

(0.41) 

(0.272) 

(-0.3) 

5'B3u 

6.26 

6.07 

(0.19) 

(0.182) 

5'B2u 

6.57 

6.21 

(0.002) 

(0.004) 

'"Ref.  152. 

''^Ref.  153. 

'^^Ref.  65. 

'^^Ref.  154 

...     s;  :/'    . 
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Table  5-14:  User  CPU  time  needed  for  the  steps 
in  the  STEOM-CCSD  calculation  with  the  polarized  basis. 


Step 


User  CPU  Time  in  Seconds 


Integrals 

Property  Integrals 

SCF 

Integral  Transformation 

Integral  Reordering 

CCSD 

Forming  Hbar 

1 1  IP-EOM-CCSD  States 

23  EA-EOM-CCSD  States 

Forming  G 

84  STEOM-CC  States 


29229 

136 

166 

2491 

'  2744 

56405 

3937 

243 

9873 

789 

3016 


CHAPTER  6 
GRADIENTS  FOR  THE  SIMILARITY  TRANSFORMED 
EQUATION-OF-MOTION  COUPLED-CLUSTER  METHOD  ,.; : 

Theory 

As  described  in  Chapter  4,  the  similarity  transformed  equation-of-motion  coupled- 
cluster  singles  and  doubles  (STEOM-CCSD)^"^^^  energy  for  an  excited  state  is  an 
eigenvector  of  the  double  similarity  transformed  Hamiltonian 

_i  (6-1) 

Since  G  is  not  Hermitian,  it  has  different  right  and  left  hand  eigenvectors.  Given  a  right 
and  left  hand  eigenvector  for  a  specific  state,  the  energy  of  the  state  can  be  represented  as 

E={0\LGR\0).     ■  (6-2) 

where  the  defining  equations  for  L  and  R  are  .      , 

{p\GR\0)  =  E{p\R\0),  (6-3) 

and 

(0|LG|p)  =  £;(0|L|p).  (6^) 

The  notation  here  is  that  the  space  of  all  possible  determinants  is  divided  into  |h)  = 
|p)  ®  |q),  as  was  done  by  Stanton  for  equation-of-motion  coupled-cluster  singles  and 
doubles  (EOM-CCSD)  gradients.^^^^  yj,g  difference  is  that  since  the  STEOM-CCSD 
eigenvector  only  consists  of  the  reference  determinant  and  single  excitations,  |p)  only 
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Ill 

includes  the  reference  determinant  (|0))  and  single  excitations  (|s)).  The  double  (|d)) 
excitations  are  a  part  of  |q). 

It  is  possible  to  write  down  an  energy  functional^'*  and  take  the  derivatives  of  the 
appropriate  quantities  in  order  to  derive  expressions  for  the  derivative  of  the  STEOM- 
CCSD  energy  with  respect  to  a  general  external  perturbation.  However,  care  must  be 
taken  since  the  STEOM-CCSD  functional  is  much  more  complicated  than  the  EOM- 
CCSD  functional.^'*  For  EOM-CCSD  gradients^°  the  quantity  Zeta  (Z)  was  introduced 
in  order  to  guarantee  the  solution  of  the  ground  state  coupled-cluster  equations  and  to 
make  the  ground  state  T  amplitudes  stationary  with  respect  to  the  perturbation,  i.e.  the 
role  of  Lambda  (A)  for  ground  state  coupled-cluster  theory. ^^^  This  Z  is  still  needed  for 
STEOM-CCSD,  but  now  two  more  Z  like  quantities  are  needed  in  order  to  make  the 
S^  coefficients  stationary  with  respect  to  the  perturbation.  These  will  be  called  Z^  and 
Z",  and  they  have  the  form 

(6-5) 


Z-  =  Zr  +  ZJ  =  Y:  d{':"",}  +  \  Y,  G{.*mit6}. 


m,b 
'J 


Even  though  5j  are  not  used  in  the  final  diagonalization  over  singly  excited  states,  they 
appear  in  the  equations  for  S}^^  and  therefore  must  appear  in  the  functional.  Since  5^ 
appear,  Z  must  also  contain  one-body  components.  The  forms  of  Zr  can  be  compared 
to  5"*"  and  S~,  which  are 


(6-6) 


S-  =  5r  +  5^-  =  J]  .15'  {mt.'}  +  i  Y:  s7;'{mUb^j}. 
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Of  course,  T  is  the  same  as  previously,  namely 


a,b 

'      In  order  to  be  consistent  with  ground  state  coupled-cluster  theory, '^'^  what  was  called  Z 
for  EOM-CCSD  gradients^°  will  be  called  A.  It  has  the  form 

A=E^^H  +  iE^:^{^Vi}-  (6-8) 


a,b 


While  this  notation  can  be  confusing,  it  should  be  clear  from  context  what  is  meant. 
The  STEOM-CCSD  functional  is 

F  =(0|LGi?|0)  +  E(r  |Z+G|r)  +  ^(^^JZ-GI^^) 

<■  m  (6-9) 

+  (0|A^|0)  +  £;(1-(0|L/?|0)). 
The  l^*")  and  !$„,)  are  the  determinants  where  an  electron  has  been  added  to  an  active 
virtual  orbital  and  where  an  electron  has  been  removed  from  an  active  occupied  orbital, 
respectively. 

The  procedure  for  deriving  gradient  equations  from  the  functional  is  the  same  as 
that  used  in  Chapter  3;  the  derivative  of  the  functional  with  respect  to  each  operator  in 
the  formula  will  be  taken  and  set  to  zero.  What  taking  the  derivative  with  respect  to  an 
operator  means  is  that  the  derivative  with  respect  to  each  coefficient  within  the  operator  is 
taken,  and  the  annihilation  and  creation  part  of  the  operator  remains.  By  inspection  it  can 
be  seen  that  taking  the  derivative  of  F  with  respect  to  L  and  R  and  setting  them  equal  to 
zero  give  Eqs.  (6-3)  and  (6-^)  respectively.  It  also  can  be  seen  by  inspection  that  taking 
the  derivative  with  respect  to  A  yields  the  ground  state  coupled-cluster  equations.  The 
derivative  with  respect  to  E  gives  the  normalization  condition  for  L  and  R.  Therefore, 
L,R,  A,  and  E  are  all  stationary  with  respect  to  the  perturbation. 


or  equivalently 
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The  result  of  taking  the  derivative  with  respect  to  Z"^  and  Z~  is  more  difficult  to  see. 
For  this  two  new  spaces  must  be  introduced.  They  are  \Q'^)  and  \Q~),  and  they  are 
defined  as  the  spaces  into  which  5^  project  when  acting  on  IP"*")  and  |P~)  respectively, 
and  where  |P^)  are  the  space  of  all  determinants  with  one  electron  added  to  an  active 
virtual  orbital  or  with  one  electron  removed  from  an  active  occupied  orbital. ^^ 

Using  this  notation,  the  result  of  taking  the  derivative  of  the  functional  with  respect 
to  Z+  is 

dF 
0  =  Q^  =  {Q^\Gm,  (6-10) 

(6-11) 

o  =  ($f|G'|r).         •;  ;  ^ 

The  derivative  with  respect  to  Z~  is,  for  all  m 

dF  ■■    ■ 

0=g^  =  {Q-\G\^m),  (^12) 

or 

(6-13) 

0   =    ($;,;|G'|$„,)- 

These  are  the  equations  for  the  parts  of  G  set  to  zero  in  the  second  similarity 
transformation,^^  and  these  equations  can  be  used  as  the  defining  equations  for  S.  The 
fact  that,  in  practice,  5  is  calculated  in  a  different  manner  does  not  matter.  It  is  sufficient 
that  these  equations  are  equivalent  to  the  ones  used. 

Before  continuing,  it  is  necessary  to  introduce  another  piece  of  notation.  The  term 
connected  has  been  introduced  previously,  along  with  its  symbol  c.  Algebraically  it 
means  that  the  two  operators  share  at  least  one  index.  Operators  which  are  not  connected 
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are  disconnected;  they  share  no  indices.  The  symbol  disc  will  be  used  to  indicate  when 
operators  that  potentially  could  be  connected  are  explicitly  disconnected.  Disconnected 
operators  containing  an  even  number  of  annihilation  and  creation  operators  commute.'^' 
Therefore,  for  the  generic  operators  A  and  B  with  an  even  number  of  annihilation  and 
creation  operators, 

(^^)disc  =  (5^W.-  ■  '      ■  (6-14) 

By  definition, 

AB  =  {AB\  +  (AB),,,^  =  {AB)^  +  (fi^),,,^.  (6-15) 

For  a  discussion  about  connectedness,  see  Ref.  116. 

In  order  to  continue,  it  is  necessary  to  expand  out  G  in  Eq.  (6-9).    It  has  been 
shown^^  that 

«=(*P}X-(P-'R-  <^'6) 

For  the  first  term  in  Eq.  (6-9),  {e^' }  is  not  used,  and  therefore  the  second  term  of  Eq. 
(6-16)  does  not  contribute.  Breaking  the  commutator  in  the  first  term  gives 

where  G2  is  used  to  show  that  this  Hamiltonian  was  built  using  {e-^^  j  Q^jy^  ^^^  jj^g  (-1)  in 
the  second  term  is  to  show  that  the  constant  part  of  the  exponential  is  implicitly  connected 
and  is  in  Eq.  (6-16).  Substituting  Eq.  (6-17)  into  the  first  term  of  Eq.  (6-9)  gives 

(0|LGi?|0)  =  (0|L^{fi^^}i?|0)  -  (0|L({e^^  -  i]h)^^^R\Q).  (6-18) 
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The  second  term  is  zero,  since  (e-^^  -  l}  excites  into  the  |q)  space,  and  L  can  only 
connect  (0|  to  |p). 

It  has  also  been  shown^^  that 

(Q±|G|P±)=(Q±|^+(^5±)JP±) 

(6-19) 
-  {Q^\S^\P^){P^\H  +  [HS^)^  -  Ecc\P^), 

where  Ecc  stands  for  the  ground  state  coupled-cluster  energy.  The  connectedness  in  the 
second  term  on  the  right  is  redundant;  S^  must  be  connected  to  H  to  give  an  nonzero 
contribution. 

Substituting  Eqs.  (6-18)  and  (6-19)  into  Eq.  (6-9)  gives 

-  ^(r|Z+5+|P+)(P+|(^(l  +  5+)  -  Ecc)m 

(6-20) 

m 

-Y.{^.r\Z-S-\p-){p-\{H{\  +  S-)-Ecc)\^^) 

ni 

'■        :  +(0|A^|0)  +  E(l-(0|Li?|0)).  ^- 

It  is  now  possible  to  continue  taking  derivatives  of  the  functional  with  respect  to  the 
operators  it  contains.  Since,  L,  R,  A,  E,  and  7r  have  been  considered,  only  5+,  S~ ,  and 
T  remain.  The  annihilation  and  creation  portions  of  these  operators  will  be  called  Q5+, 
Vts-,  and  9.T,  respectively.  Remember,  these  are  the  parts  of  the  operators  remaining 
after  the  derivative  with  respect  to  the  operator  has  been  taken. 
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The  first  term  to  consider  is  5"*".  Taking  the  derivative  of  the  functional  with  respect 

to  it  gives  for  all  /  and  k,  where  ac  is  a  generic  \Q)  space  determinant, 

dF        dF 


0  = 


95+       dS'j 


+  5^($nZ+(^f^5;)j$^)<^e/  (6-21) 

e 


or 


0=(0|LF|e^^f^5+}7?|0) 

+  ($^|Z+(^l]5;)j$^)  -  ($^|Z+5+|P+)(P+|^|$«)  (6-22) 

-  ^($^|Z+|$«)($/|(^(1  +  5+)  -  £;cc)|$'=). 

e 

Inserting  a  resolution  of  the  \Q)  space  into  the  second  and  third  terms  yields  for  all/and  k 

+  ^($^|Z+|$")($"|(^n5.)|$/)  ;;-^ 

-  ^($^|Z+|$'')($''|5+|P+)(P+|^|$'^)  ^^^^^ 

-  J]($^|(^(l  +  5+)  -  Ecc)W)mz-^\^''), 

e 

where  r]  is  another  generic  |Q)  space  determinant.  In  matrix  form,  Eq.  (6-23)  is 

+  (P+|Z+|Q+)(Q+|5+|P+)(P+|^|g+)  (6-24) 

-(P+|Z+|Q+)(Q+|(^Q5,)  |P+). 
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Taking  the  derivative  with  respect  to  S~  and  following  the  same  procedure  as  above 
gives 

{0\LH{e'^ns-]R\0)  ={P-\{H{1  +  S-)  -  Ecc)\P-){P-\'^-\Q-) 

+  {P-\Z-\Q-){Q-\S-\P-){P-\H\Q-)  (6-25) 

-{p-\z-\Q-){Q-\{ms-)jp-). 

Equations  (6-24)  and  (6-25)  can  be  solved  iteratively  to  calculate  7r. 

Before  taking  the  derivative  of  the  functional  with  respect  to  T,  it  will  be  useful  to 
remember  that 

^  =  ^  ('"^^'^)  =  "^'^^^^^T  -  nre-^He:^  =  [H,  ^t]  .  (6-26) 

The  derivative  of  Eq.  (6-20)  with  respect  to  T  is  then 
Q=Qf=mlH.^T]{e'^}R\0) 

e 

(6-27) 

m 

-Y,{^m\z-s-\p-){p-\[H,nT]{i  +  s-)\^m) 

m 

+  (0|A[^,r2r]|0). 
Focusing  on  the  last  term, 

(0|A[^,Or]|0)  =  (0|A^J]r|0)  -  (0|Al]r^|0) 

=  (0|A^|s  +  d)  -  (0|A|s  +  d)(0|^|0)  (6-28) 

=  (0|A|s  +  d)(s  +  d|(^-l£;cc)|s  +  d), 


(6-29) 
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where  1  is  a  unit  matrix  of  the  appropriate  rank.  Inserting  Eq.  (6-28)  into  Eq.  (6-27)  gives 

0={0\L[H,nT]{e^'}R\0) 

e 

-^($''|Z+5+|P+)(P+|[^,l^r](l  +  5+)|r) 

e 

m 

-J2{^m\Z-S-\P-){p-\[HMT]{l  +  S-)\^,n) 
m 

+  (0|A|s  +  d)(s  +  d| [H  -  lEcc) |s  +  d), 
and  solving  for  A  yields 

(0|A|s  +  d)  =  -[(0|L  [H,  Or]  {e^^  }/?|0) 

e 

-  5]($''|Z+5+|P+)(P+|  [^,  Or]  (1  +  5+)  l^'^) 

e 

+  ^($^  |Z-  ( [H,  Or]  +  ( [^.  ^r]  5-)  J  |$™ ) 
-^($,„|Z-5-|p-)(p-|[F.0r](l  +  5-)|$„,)] 

With  these  definitions  of  Z*  and  A,  the  derivative  the  functional  with  respect  to 
all  of  the  terms  which  appear  in  it  is  zero,  and  the  generalized  Hellmann-Feynmann 
theorem"'^'  '^^  can  be  applied.  Therefore,  the  derivative  of  the  STEOM-CCSD  energy 


(6-30) 


(6-31) 


119 
with  respect  to  any  perturbation  x  is 

|.|H0|.H{e-H0) 

e 

e 

+  5]($„,|Z-(F^'+(^^5-)J|$-) 

m 

-J2{^m\Z-S-\p-){p-\H^{l  +  S-)\^m) 
m 

+  (0|A^^|0). 
with 

H^  =  e-^^e^.  (^32) 

OX 

So  far,  the  derivation  has  been  for  excitation  energy  STEOM-CCSD  (EE-STEOM- 
CCSD)^  gradients,  but  their  extension  to  double  ionization  potential  STEOM-CCSD 
(DIP-STEOM-CCSD)^^  or  double  electron  attachment  STEOM-CCSD  (DEA-STEOM- 
CCSD)^^  is  straightforward.  The  only  difference  in  the  algebraic  equations  is  that  for  DIP- 
STEOM-CC,  5+  does  not  appear,^^  and  therefore  Z*  will  not  appear  in  the  functional. 
Likewise,  for  DEA-STEOM-CC,  S~  does  not  appear,^^  which  leads  to  Z'  not  being  in 
the  functional. 

One  extra  consideration  arises  with  STEOM-CC  gradients.  In  normal  coupled-cluster 
gradients,  the  equations  for  the  response  of  the  orbitals  to  the  perturbation  are  simplified 
by  using  the  fact  that  the  method  is  invariant  to  rotation  among  the  occupied  orbitals 
and  among  the  virtual  orbitals. '°^  In  STEOM-CC,  because  of  the  splitting  of  the  orbitals 
into  an  active  set  and  an  inactive  set,  this  invariance  to  orbital  rotations  does  not  hold. 
The  solution  is  simple,  though.    In  the  derivation  of  gradients  for  the  case  that  some 
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orbitals  are  not  correlated, ''^^  the  extra  terms  that  arise  have  already  been  included,  and 
the  same  code  can  be  used  here. 

Applications 

In  order  to  test  how  well  STEOM-CCSD  describes  excited  state  potential  energy 
surfaces,  the  diatomic  molecule  examples^^  from  Chapter  3  were  tested.  The  STEOM- 
CCSD  gradients  are  an  iterative  nl^^N^-^^  method,  both  from  the  ground  state  CCSD 
and  from  the  calculation  of  the  excited  state  A.  Therefore,  the  only  fair  comparison 
for  STEOM-CCSD  is  EOM-CCSD.  The  STEOM-CCSD  gradients  will  be  about  half  as 
expensive  as  EOM-CCSD  gradients.  Gradients  for  the  STE0M-PT2  method^  should 
have  roughly  the  same  cost  as  P-EOM-MBPT(2)'"^  gradients.  '  ■ 

The  active  space  used  for  H2  consisted  of  the  l^a+,  1-3^-,  Itt,,,  Itt,,  and  \8g 
orbitals,  giving  a  percent  active  space  of  99.7%  for  the  l^S"  state  at  the  excited  state 
minimum  geometry.  The  active  space  for  BH  was  \-la  and  l-37r,  giving  a  98.2%  active 
space  component  for  the  1  ^n  state  at  its  minimum.  For  CO  the  active  space  was  4-10(7, 
l-47r,  and  \6,  giving  a  98.7%  active  space  component  for  the  I'n  state  at  its  minimum. 
For  N2  the  active  space  was  3-6a+,  2-5^-,  l-27r„,  l-27r,,  and  Ug,  giving  a  99.2% 
active  space  component  for  the  l^n,  state  at  its  minimum.  For  BF  the  active  space  was 
4-lOa,  2-47r  and  \8,  giving  a  98.9%  active  space  component  for  the  I'n  state  at  its 
minimum.  Finally,  the  active  state  for  C2  was  3-5a+,  2-3^-,  l-27r„,  l-2x,,  and  \6„ 
giving  a  99.8%  active  space  component  for  the  l^n„state  at  its  minimum. 

The  active  spaces  were  chosen  by  including  all  orbitals  for  which  the  corresponding 
IP-EOM-CCSD67  or  EA-EOM-CCSD^^  state  could  be  converged  and  for  which  the  state 
consisted  of  greater  than  seventy  percent  single  excitation  character.  Since  the  decoupling 
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caused  by  the  second  similarity  transformation  depends  on  the  52  coefficients  being 
small,  any  IP-EOM-CCSD  or  EA-EOM-CCSD  states  with  less  than  seventy  percent 
single  excitation  character  are  automatically  excluded  from  the  transformation. 

The  vertical  excitation  energies  for  each  molecule  are  listed  in  Table  6-1.  The 
aug-cc-pVTZ  basis  set^^'  was  used  for  all  of  the  calculations,  and  only  the  spherical 
components  of  the  d  and  /  functions  were  kept.  The  ground  state  geometries  used  were 
the  aug-cc-pVTZ  CCSD  geometries  from  Table  3-3.  Except  for  N2.  the  EOM-CCSD 
and  STEOM-CCSD  vertical  excitation  energies  agree  pretty  well. 

The  adiabatic  excitation  energies  are  given  in  Table  6-2.  The  EOM-CCSD  and 
STEOM-CCSD  are  very  close  for  BH,  BE,  and  C2.  For  H2  the  EOM-CCSD  and  STEOM- 
CCSD  straddle  the  experiment,  with  the  EOM-CCSD  being  a  little  closer.  For  CO  the 
STEOM-CCSD  is  0.15  eV  lower  than  the  EOM-CCSD,  and  for  N2  the  STEOM-CCSD 
is  0.28  eV  lower.  In  both  cases  the  STEOM-CCSD  agrees  better  with  the  experimental 
adiabatic  excitation  energy. '^^  Presumably,  for  excitations  out  of  the  triple  bond,  the 
better  description  of  higher  excitations  given  by  STEOM-CCSD  is  important. 

The  equilibrium  distances  for  the  excited  states  are  given  in  Table  6-3.  Except  for 
H2,  the  EOM-CCSD  and  STEOM-CCSD  geometries  are  quite  close.  It  is  not  clear  why 
H2  shows  such  a  big  discrepancy.  EOM-CCSD  for  H2  is  exact  within  a  basis  set,  so 
the  difference  must  be  attributed  to  a  shortcoming  in  the  STEOM-CCSD.  The  fact  that 
the  rest  of  the  equilibrium  bond  lengths  agree  well  with  the  EOM-CCSD  and  with  the 
experiment  gives  confidence  that  the  problem  with  H2  is  an  anomaly. 

That  the  H2  excited  state  vibrational  frequency,  as  given  in  Table  6-4,  does  not  agree 
with  the  EOM-CCSD  frequency  should  be  attributed  to  the  problem  with  the  bond  length. 
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Even  still,  its  agreement  with  the  experimental  vibrational  frequency  is  satisfactory.  For 
all  of  the  molecules,  the  error  in  the  vibrational  frequency  is  about  the  same  for  the 
two  methods. 

One  final  note  should  be  made  about  the  dependence  of  the  STEOM-CCSD  results  on 
the  choice  of  active  space.  Conspicuous  by  their  absence  in  the  active  space  for  BF  are 
the  (occupied)  Itt  orbitals.  They  are  not  included  because  at  the  ground  state  minimum 
geometry,  the  EA-EOM-CCSD  state  which  corresponds  to  them  is  only  about  fifty  percent 
single  excitations.  However,  at  the  excited  state  minimum,  the  percentage  is  just  over 
seventy  percent.  Some  test  calculations  were  run  to  determine  the  effect  of  including 
them  in  the  active  space.  When  that  was  done,  the  adiabatic  excitation  energy  dropped 
to  6.331  eV,  in  better  agreement  with  the  experimental  adiabatic  excitation  energy.  On 
the  other  hand  the  excited  state  bond  length  decreased  to  1.291  A,  and  the  vibrational 
frequency  decreased  to  1 1 83  cm"' ,  both  in  noticeably  worse  agreement  with  experiment. 
The  strangest  part,  though,  was  that  adding  the  functions  to  the  active  space  decreased 
the  percent  active  space  component  to  97.8%  from  98.9%.  It  is  not  clear  why  adding 
functions  to  the  active  space  would  decrease  the  percent  active  space  component.  But 
this  does  show  that  the  issue  of  active  spaces  in  STEOM-CCSD  needs  further  study. 
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Table  6-1 :  Vertical  excitation  energies  (in  eV)  for  the  excited  singlet  states 


H2 

BH 
CO 

N2 
BF 

C2 


EOM-CCSD^ 

STEOM-CCSD 

12.717 

12.723 

2.914 

2.934 

8.666 

8.552 

9.514 

9.302 

6.454 

6.469 

1.308 

1.285 

«)Ref.  93. 
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Table  6-2;  Adiabatic  excitation  energies  (in  eV)  for  the  excited  singlet  states 


EOM-CCSD^ 

STEOM-CCSD 

Expt.'' 

H2 

11.353 

11.422 

11.3694 

BH 

2.913 

2.934 

2.8685 

CO 

8.256 

8.107 

8.0684 

N2 

8.839 

8.560 

8.5900 

BF 

6.417 

6.425 

6.3427 

C2 

1.103 

1.098 

1 .0404 

''>Ref.  93. 

■  •:  .? 

''^Ref.  122. 
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Table  6-3:  Equilibrium  distances  (in  A)  for  the  excited  singlet  states 


EOM-CCSD^ 

STEOM-CCSD 

Expt.''        : 

H2 

1.283 

1.258 

1.2928 

BH 

1.211 

1.214 

1.2186 

CO 

1.224 

1.232 

1.2353 

N2 

1.202 

1.208 

1.2203 

BF 

1.304 

1.308 

1.3038 

C2 

1.318 

1.314 

1.3184 

')Ref.  93. 

'^'''•'■-..' 

''^Ref.  122. 
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Table  6-4:  Harmonic  vibrational  frequencies  (in  cm"')  for  the  excited  singlet  states 


EOM-CCSD^ 

STEOM-CCSD 

Expt.^ 

H2 

1368 

1344 

1358.09 

BH 

2372 

2371 

2251.0 

CO 

1592 

1586 

1518.2 

N2 

1854 

1837 

1694.21 

BF 

1279 

1248 

1264.9 

C2 

1630 

1652 

1608.35 

^)Ref.  93. 
"^Ref.  122. 


CHAPTER  7 
A  FINAL  WORD 

In  this  dissertation  two  new  methods  were  presented  for  calculating  excitation  ener- 
gies of  atoms  and  molecules.  The  new  methods  are  the  partitioned  equation-of-motion 
coupled-cluster  singles  and  doubles  (P-EOM-CCSD)  method  and  the  partitioned  equation- 
of-motion  second  order  many-body  perturbation  theory  (P-E0M-MBPT(2))  method."'* 
These  methods  were  shown  to  have  great  promise  for  calculating  vertical  excitation  ener- 
gies, but  they  consistently  were  a  little  worse  for  valence  excited  states  than  for  Rydberg 
excited  states. 

Gradients  for  the  P-E0M-MBPT(2)  method  were  developed.  Being  able  to  ana- 
lytically calculate  the  derivative  of  the  energy  with  respect  to  an  external  perturbation 
(such  introducing  an  electric  field  or  moving  an  atom)  not  only  allows  for  the  efficient 
and  accurate  calculation  of  one  electron  properties  for  the  state  of  interest,  but  also  it 
makes  possible  efficient  searches  of  potential  energy  surfaces  for  extrema,  such  as  min- 
ima or  transition  states.  Ultimately,  if  a  computational  method  is  going  to  be  useful  for 
describing  a  variety  of  chemistry,  it  must  have  gradients  available. 

Being  able  to  search  potential  energy  surfaces  also  provides  more  demanding  tests 
for  a  computational  method.  If  structures  are  to  be  determined,  it  is  not  sufficient  to  be 
able  to  accurately  describe  just  one  region  of  the  potential  energy  surface.  Any  method 
that  can  only  accurately  describe  a  small  region  of  the  surface  is  of  limited  use,  since 
it  may  be  difficult  to  know  a  priori  whether  or  not  the  structures  of  interest  lie  within 
the  well-described  region. 
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In  the  end,  the  performance  of  the  P-E0M-MBPT(2)  method  for  excited  state  surfaces 
was  mixed.  The  method  was  almost  always  qualitatively  correct,  but  only  sometimes 
quantitatively  correct.  Because  of  its  reduced  cost  relative  to  EOM-CCSD  gradients,^^ 
the  method  is  still  useful  as  a  relatively  simple  and  inexpensive  way  to  account  for 
electron  correlation  in  excited  state  energies  and  structures  and  as  a  way  to  generate  a 
good  initial  guess  for  an  EOM-CCSD  structure  search. 

Since  an  inexpensive  method  that  gives  consistently  reliable  results  for  excited  state 
structures  is  still  needed,  focus  shifts  to  the  new  similarity  transformed  equation-of- 
motion  coupled-cluster  singles  and  doubles  (STEOM-CCSD)  method.^  As  an  example 
of  the  power  of  the  STEOM-CCSD  method,  the  vertical  excited  state  spectrum  of  free 
base  porphin  was  studied.  This  study  shows  how  accurate  calculations,  even  without 
including  the  effects  of  geometry  relaxation,  can  be  used  to  interpret  complicated  spectra. 
The  ability  of  STEOM-CCSD  to  calculate  many  excited  states  quickly  was  important  in 
that  it  allowed  the  calculation  of  essentially  the  entire  spectrum  in  the  energy  range  of 
interest.  These  calculations  also  showed  the  importance  of  care  in  choosing  the  details 
of  the  calculation,  such  as  the  basis  set  and  the  active  space. 

Finally,  gradients  for  STEOM-CCSD  were  developed.  Because  of  the  complicated 
structure  of  the  STEOM-CC  method,  the  gradient  equations  are  tedious.  But  eventually 
the  equations  take  on  a  very  nice  structure.  Initial  tests  on  very  simple  molecules  look 
promising,  but  it  is  clear  that  many  issues  involving  the  role  of  the  active  space  in 
STEOM-CC  remain  to  be  explored. 
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